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ABSTRACT 


Some new results in the scheduling of dynamic multi- 
commodity flows in data communication networks are presented. 

A new performance measure for effective delivery of back- 
logged data to their destinations is defined and the solution 
to the resulting delivery problem is obtained through a se- 
quential linear optimization methodology. Properties of an 
optimal dynamic multicommodity flow schedule are studied in 
detail, taking advantage where possible of the linear pro- 
gramming formulation. The special case of the delivery 
problem in a single destination network also is analyzed. 

Application of the results to stochastic delivery prob- 
lems in which the data inputs to the network are modelled as 
Poisson processes is addressed, and a new dynamic data com- 
munication network analysis is presented. 

Finally, the delivery problem on networks with capacitated 
links and with traversal delays is considered and some new 


results obtained. 
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I. INTRODUCTION 


A. MOTIVATION 

The design of data communication networks has received 
much attention during the past decade and the interest in 
pastel is constantly growing. This fact is explained by 
the rapidly expanding role being played by data processing 
in today's society and the apparent advantage in sharing power- 
ful computational resources. 

The overall subject of data communication networks can be 
looked upon from many different points of view, each repre- 
senting an intellectually challenging problem. Just to name 
two: there is the problem of cost effective topological net- 
work design and there is the problem of determining, in a given 
network, the routes which data should follow from their source 
to their destination. The problem of assigning routes has 
been one of the most intensively studied areas in the field 
of data communication networks. 

It is possible to view the complete routing problem as com- 
Besecsor two parts: the first, calculation of routing assign- 
ments and second, their implementation in an actual data 
cA Control strategy in a real, basically stochastic, net- 
work environment. It is the first part that we are mostly con- 
cerned with in this research, although we touch briefly on the 
other. 

A network model for calculating routing assignments can 


be roughly classified as static or dynamic. These reflect the 


15 








long and the short term stochastic behavior of a network, 
respectively. A finer taxonomy is possible, regarding the 
frequency with which a routing assignment has to be recomputed 
to support efficient communication. A desirable complete rout- 
ing methodology must provide fast adaptation to a changing net- 
work environment and be stable at the same time. It should 

be dead-lock free as well as distributed to minimize computa- 
Plenauzesmplexity and control information flows. Clearly, there 
are substantial difficulties in achieving all these attributes 
Simultaneously or even in qualifying trade-off relations among 
them. 

In spite of the progress made in some areas, and in par- 
ticular in the computation of static routing assignments, many 
basic questions remain unanswered. For example, definition of 
a meaningful performance measure for dynamic routing assignments 
which will give rise to mathematically tractable problem formu- 
lation, seems to have been misSing. Thus, it appears that new 
insight into the relationships that make up a dynamic network 
model should be helpful in coping with these questions. The 
new conceptual ideas and analytic tools for studying them 


introduced in this research are a step in this direction. 


B. SURVEY OF PREVIOUS WORK 

In {15]) Kleinrock introduced an analytical model of a data 
communication network. Many algorithms use this model to com- 
pute routing assignments, for example, the algorithms due to 


Cantor and Gerla [17], Fratta, Gerla and Kleinrock [16] and 
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Defenderfer [18]. In all those algorithms the steady-state 


delay of messages in each link is calculated explicitly as 


where: 
57 o the data flow rate in link [i J] (messages 
J second); 
Sis is the capacity of link [i,j] (messages/second); 


is the expected delay/message experienced by 


CaN 
A all messages using that link. 


A routing assignment is selected to minimize the expected weighted 


delay T, 


m = 


pos E 

Bo 1J 1] 

of messages traversing the network. This analysis is based on 
the assumption that message arrivals to each link can be con- 
sidered as Poisson process with independent, exponentially 
distributed, message length, which requires Kleinrock's famous 
"independence assumption" that messages "lose their identity" 
at each node and are assigned new independent lengths. 

All the algorithms referenced thus far exploit the convexity 
of the objective function. A different approach, leading to 
a linear problem formulation, was suggested by Wozencraft in 
[29]. Instead of minimizing the average expected delay, the 
objective is to minimize the maximal saturation ratio a 


j 


E7 








in the network. Subsequently the next maximal saturation ratio 
is minimized, etc. This Min-Max policy was extensively studied 
by Ros [8] and as a result considerable insight into the problem 
of static routing was gained. One of the objectives of this 
research is to study possible generalization of the Min-Max 
approach and in particular its sequential character to deal with 
dynamic network models. 

Kleinrock's steady-state and independence assumptions make 
his model appropriate to describe long term stochastic behavior 
of a data communication model, but less applicable for computa- 
tion of dynamic routing assignments. 

One of the most important efforts in the direction of find- 
ing a more appropriate model for describing network dynamics 
was proposed by Segall [1]. His approach is to view the con- 
tents of data queues at the network nodes in a continuous state 
space setting rather than as an integer number of messages or 
bits. The continuous nature of the state is justified by 
recognizing that any individual message contributes very little 
to the overall behavior of a network. Having established the 
model, Segall expresses the dynamic routing assignment problem 
as a linear-optimal control problem for which a close-loop solu- 
tion is sought. This formulation has been used to investigate 


how to minimize 
k 
D = rent 


where: 
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te is some time at which all queues are empty; 


as (t) is the queue size at node i of data destined to 
node k, at time t. 

Here D is the total delay experienced by the queued messages 
during the period [torte], during which all the queued data 
aeRO be delivered to their destinations. 

The solution to this problem is approached by means of 
the minimum principle of Pontryagin [30], since this principle 
provides not only necessary but also sufficient condition for 
optimality in that case. A comprehensive set of principles 
describing the closed loop solution has been obtained by Moss 
(2] for the case in which all backlogged data have the same 
destination and the input to the network is continuous and con- 
stant. Additional study of the single destination case is pro- 
vided in [ll] and [31]. Unfortunately the optimal control 
approach runs into difficulties when the general, multicommodity 


case is considered and no general solution has yet been obtained. 


erEYNOESIS OF THESIS 

To cope with the difficulties posed by the multicommodity 
dynamic routing assignment problem, a different approach is 
taken in this research. It involves a change in the perform- 
Miecesmeasure, Rather than trying to minimize the total delay 
experienced by the backlogged data, the concept is to deliver 
Mo: that data to their destinations in the shortest time 
possible. Since eventually there are many dynamic routing 


assignments that can do so, one which maximizes the total amount 
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of delivered data over time, is selected. We refer to this 
problem formulation as the "delivery problem" and to the corres- 
ponding dynamic routing assignment as the "flow schedule." 

There are several advantages associated with the above state- 
ment of the delivery problem. Most important is the ability 

to find the desired flow schedule by solving a sequence of 
hierarchically related linear programs. Also, it turns out 

that the new performance measure reveals some structural Proper- 
ties that bring new insight to the problem of dynamic flow 
scheduling. The main purpose of this research is to study those 
properties in detail. 

The new performance measure (we call it the "delivery func- 
tion") and related optimality criterion are explained in Ch. II. 
BIlso, by introducing a concept of "global optimality" we are 
able to relate the optimal delivery function to the "total de- 
Em ceriterlon. In Ch. If we also introduce the basic network 
model and the notational convention to be used throughout the 
ehesis, 

In Ch. III a solution algorithm to the delivery problem 
is presented. It consists of solving a sequence of hierarchi- 
cally related linear programs. In principle, each linear 
program that is solved contracts the space of remaining feasible 
flow schedules, until finally an optimal piecewise constant 
flow schedule is obtained. The optimality of the resulting 
flow schedule is formally derived in Ch. III using some basic 
results regarding the properties of feasible flow schedules. 

In Ch. IV we study in detail the structural properties of 


an optimal flow schedule. By exploiting various properties of 
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linear programming we are able to derive several results that 
characterize an optimal flow schedule. One key result is the 
description of critical sets of commodities and the capacity 
resources (links) they must saturate for various periods of 
time. Also, the properties of those saturating flows, and in 
particular their total rates, are determined. In this chapter 
the important idea of optimal solution "stability" is intro- 
duced, which allows one to express most of the above properties 
in terms of the optimal dual variables associated with the 
linear programs of the solution algorithm. 

We devote Ch. V to discussion of the single destination 
Poser The solution algorithm of Ch. III is specialized to 
handle this case. Considerable additional simplifications are 
obtained by observing that the single destination delivery 
problem may be interpreted as a single commodity flow problem 
so that advantage may be taken of many well established results. 
The concept of stability 1s revised and exploited as part of 
the solution algorithm. 

Continuing the discussion from Ch.II we show that the opti- 
mal single destination flow schedule is also globally optimal 
and thus also solves the single destination "minimal total de- 
lay problem" [11]. The computational advantages of the new 
algorithm are addressed briefly. 

In Ch. VI we analyze the stochastic delivery problem. Here, 
in addition to the backlogged data considered so far, we are 
concerned with Poisson arrivals of messages to the network. 


Following Yee [12] the expected minimal time to empty an M/M/1 
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queueing system is taken to be the new performance measure. 
It is shown that the theory of dynamic flow scheduling, derived 
deerit the deterministic case, and in particular the se- 
quential linear optimization methodology, can be applied (at 
least, in principle) to solve the stochastic delivery problem. 

In Ch. VII we consider a more general setting for a de- 
livery problem. Here we associate with each link a traversal 
delay, in addition to a capacity constraint. Although it would 
be possible to continue our discussion within the context of 
data communication networks, it is more natural to choose the 
transportation problem as the framework for our investigation. 
This allows a generalization of link capacity constraints to 
include loading and unloading constraints. The addition of 
traversal delays greatly complicates the delivery problem. 
It is possible, however, to exhibit a (conceptually) simple 
solution procedure for the case of bi-partite transportation 
networks. A discrete time approximation for general network 
models also is discussed and a particular example of military 
application is presented. 

In Ch. VIII we summarize the most important results of 
this research and indicate areas for future study. 

Finally, we defer to the Appendix a number of proofs, 
examples and short discussions which would tend to blur the 


Maınsıdeas if left within the body of the thesis. 
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II. PROBLEM FORMULATION 


A. COMMUNICATION NETWORK MODEL 
l. Topological Representations 
A data communication network may be modelled as a set 
of nodes interconnected by a number of links. The nodes repre- 
sent physical locations at which data may enter or exit the 
network and the links represent unidirectional channels over 
which data is transmitted from node to node. A typical data 


communication model is shown in Fig. II.l. 


<—l— 
oe 


data entry 


Er 


data exit 


Face... Ccomnmuntceacton Network 


With each link we associate a channel capacity which indicates 
the upper bound on the data flow rate for that channel. With 
each nođe we associate, at every instant of time, the amounts 

of data awaiting transmission to each destination at the corres= 


ponding location. The collection of node descriptors for all 
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the nodes in the network constitutes the state of the system, 
or equivalently, system congestion at any given time t. 

We will say that any data in the network is commodity 
E its origin (entry node) was node i and its final 
destination (exit node) is node k. We shall also say that any 
data in the network is commodity k if its final destination is 
node k. 

Consider a data communication network G(V,Lp) » where 


Be 2, ...,n} is a set of n nodes and L. = [[i,jl} is a set 


0 
Or links. By the notation [i,j] we mean the link that con- 
Heeeoepode 2 tO node j, in that direction. We will also use 
mie notation Ly = ich 2p Co and ae note an element of Ly by 
e. If link e corresponds to link [1,3] then h(e) = i and 
t(e) = 5.7 We say that node i can communicate to node k iff 
there is at least one directed chain of links going form node 


i to node k. We also define: 


Ny = set {(i,k)} of node pairs such that node i can 
communicate to node k, 1 # k. 
EE) amount of commodity k, stored in node llat time t, 
7 (i,k) € Na. 
0 
Een Mo rate oF commodity k on link [1,;]] at timeret, 
2) oa = Ee and Na 
5 Sapa rot Line lla e Lo: 
a(t) Destroy rate of commodity k arriving at node 1 from 
3 outside the network, at time t, *(i,k) e€ No: 
h=-head, t=-tail hie) > tle) 
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We reserve the use of respective capital letters for sets and 
vector notation, interchangeably. For example, 
P(e) È (£f (t), f£? (t) f% (t)...) denot 

12 EE BERR i TE enotes a vector (set) 


of flows. 
2. Dynamic System Equations and Constraints 
The flows and the queues just defined must satisfy three 
basic constraints: non-negativity, conservation and capacity. 
Memon negativity Constraint states that 


ee 0, “[i,j] € La. “(i,k) € Nn and “t. GN 


OF 0 


The conservation constraint may be written as 


k i K 
o a or 


I 


meee) e 2 1 Ea) 
jt#i) j(#i) E 
= a,(t), otherwise, 
(i,k) € No andes (AZ) 


eaS traint (11.2) accounts for the fact that at all times the 
amount of any commodity stored at any node is a non-negative 
quantity. Moreover, the fact that the net delivery rate of 
commodity k from node i is non-negative, ee UL No 


and ¥t, implies that data is not stored at intermediate nodes 


o 


! 


en route from its entry node to its exit node in the network. 


"Intermediate storage of commodities is discussed in 
Appendix C. 
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Finally, the capacity constraint is 


Mh 
No 


k sr | 
ij (È e < Sizr *li,j] e Ly and «t (11.3) 


where Ba denotes the aggregate flow rate on link (A aE 


time t. 


Detintcion (11.1). 
Peset vor: flows P(t) , to SS ti is a feasible multicommodity 
flow schedule if it satisfies constraints (Il) ALE for arm 


E € [tyrt,]- P 


— 


We will assume, for mathematical convenience, that data input 
flow rates are identically zero during the time interval under 


consideration, i.e. 


k 
a, (t) 


tht 
O 
~- 
+ 
H- 
A 
m 
Z 


Se (DR, 


In Ch. VI the behavior of a communication network in a stochas- 
tic environment is considered and as (t) = ar, (el N will 
be interpreted as a rate of a Poisson process. 

We follow the model proposed by Segall [1], where the 
contents of the queues at the nodes are viewed as continuous 
quantities, rather than as integer number of messages (in Ch. 
VI we recognize the existence of separate messages but model 
their size as a continuous quantity). This macroscopic point 


of view not only provides a model that is analytically simpler 


than others, but also is justified by recognizing that any 
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individual message (bit) contributes very little to the over- 
all behavior of the network. 


Let Q(t) and Q(t,) be the system states at times to 


and ti respectively. We say that the state Q(t,) is reachable 
from the state Q (to) if there exists a feasible multicommodi ty 


flow schedule F(t), t Sees ty for which 


0 


t 
k ER l k | 
A E . | ritt)dt, “(i,k) e Ny (Dies) 


0 


Or in shorthand notation, F(t): Q (to? = Q(t). To every feasi- 
ble multicommodity flow schedule we adjoin a delivery function 


Die), Es ER t which represents the total amount of data 


delivered to their destination by time t. 


k k 
D(t) = } lese. (t)i, amt Se (II.6a) 


or equivalently 


D(t) = J f ríta)da, ty <t<t, (II.6b) 


From the nature of the model that we constructed and 
particular from the fact that qs (ty) Sa LK Nor we 
conclude that the zero state is reachable from any other finite 


state within some finite time Es y i.e. there exist both ty and 


e 
'To avoid excessive notation complexity we will assume in 


the sequel that o) > 0 < (i,k) <e Na. The case of 
AN EROS (1,K) € No is included in the examples. 
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Et), Eg £t < t] such that F(t): Q(t) 7 0(t,) = 0. Consequently, 


for every initial state there is some minimal value of ty - We 


define the minimal total delivery time E? 


1 3S 


= 0} a) 


ar 


= min O Olt,) + A(t 


) 
ee} 5 


We conclude this section with a basic result regarding the 


nature of feasible multicommodity flow schedules. 


Theorem Br 


Bert Q (to? and Q(t) be any two states such that Q(t) is 


reachable from Q(to) - Then there exists a feasible multi- 
commodity constant flow schedule F(t) = F, to SEAS ty such 


enart P(t): Q(t) > Q(t). 


E 
nd 


This result implies that in order to transfer a system into a 
reachable state it is sufficient to look for an appropriate 
flow schedule within the subset of constant flow schedules. 
The benefits of this property will become evident in following 
chapters. 

3. Piecewise Constant Flow Schedules 

Prom this point on we will narrow our interest to the 

subset of feasible multicommodity flow schedules which are 
piecewise constant. By an M-part constant-flow schedule 


Ey (Y) y E e E we mean 


0 1l 


Tsee pewendix A, for proof Of this theorem. 
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Fylt) = Fu m, E et (II.3) 


Consulting (II.6) we immediately conclude that the correspon- 


wog delivery function Dy(t) y to SEL E, is piecewise linear. 


We write 
Mi AOS Ce 
Dy (+) = Du (Em? 70M (0D (tat)? Em+1 Se 8 En (12292) 
Do U DE Or oe 
where o,,(m) is the total delivery rate in the m-th interval, 


M 


k 
er 


o., (m) a ; 


M (m) 117.90) 


(i,K)eN, 


The rm), (i,k) e N, is the net delivery rate of the k-th 


0 
commodity from node i (see (II.2)) in the m-th interval, 
corresponding to the flow schedule segment Fy, (m) ay ek ee). 


Clearly, (m) is the slope of the delivery function in that 


PM 
interval. An example of a piecewise linear delivery function 


is shown in Fig. II.2 
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N 








Br I Ez t2 ty t 


Fig. 11.2. Piecewise Linear Delivery Function 


4. Example 
To fix ideas, in this subsection we provide a simple 
example of a communication network, and use it to illustrate 
the various notation and definitions introduced previously. 


Consider the network G(V,L,) shown in Rige LL. 


gy (E) = 10 
C19 = S53 = Sq, = Le 
1 
le) =5 q (Ep) = 5 


Fig. 11.3. Communication Network with Queued Data 
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For this netowrk 


VS AS) 

A IO 

No aa az 
O00) = LO, 10,5,0,075) 


Consider a constant flow schedule Pitt), US 


given by its components: 


ean au eee 1.41 a A 
a a ara Se 
2 wh 
Ze: 


It is easy to check that F, (t) satisfies constraints (II.1)-(II.3) 
cr all t, t e [0,15] and that F(t): A O = One 
find it useful to decompose the flow pattern into commodity 


evan £lows as shown in Fig. II.4. 


a (0) = 10 


2/3 192 
q) (0) = 5 q3 (0) 5 


Fig. 11.4. Chain Flow Decomposition 
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Using (II.2) we find the net delivery rates: 


2 _ Su: x i l 3 ji 
A E A E A UA 
nt 
AS 
The total delivery rate is thus (their sum) 97 (1) - $ and the 
Zeswueing delivery function is given in Fig. Il.>. 
D, (t) 
20 
15 
_ 4 
10 
5 
0 5 10 15 e 


Hoer S. Delivery Function 


B. OPTIMAL DELIVERY FUNCTION CONCEPT 

Recall that our objective is to find a flow schedule that 
will efficiently deliver a given set of data backlogs to their 
destinations. Thus, we need to establish what are the desired 
properties of an optimal delivery Eunice ton tanda na eS 
generating flow schedule. In defining the optimality criteria 


we were led by the following goals: 


52 





2) physically meaningful criterion; 

Ti) computational tractability of the solution procedure; 

(111) gaining new insight into the problem. 
We will show that the criterion we have chosen actually attains 
mese goals. In this section we discuss property (i) above. 
Properties (ii) and (iii) will be considered in the next 
chapters. 

l. Convex Delivery Functions 
Consider some non-convex piecewise linear delivery 


Function Dy, (t) - A typical example of such a function is shown 


mir rg. 11.0. 


D ie) 





Fig. 11.6. Non-Convex Delivery Function 


MOE.) and Q(t.) be the system states at times E; and Es: 
1 


respectively. Since Sen) is reachable from Q(t.) (by assump- 


tion), then according to Thm. 11.1, it is possible to construct 


a constant flow schedule F(t), e: 


<< fr sucn- that 
I == 
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F(t): Q(t) > Q(t,). To this flow schedule corresponds a 
linear delivery function, indicated by the broken line in 


ma. 11.6. 


We define a new flow schedule Py-jyit), t st <t 


0 E 


as 


Bel Eo E AE An EE 


ct 
| 


rr] 
Pa 
ct 
ch 
N 
ct 
IA 
ct 


where Ey (t) is the original flow schedule corresponding to 


ols 
Dy (t) - The new delivery function Dy-=1 18) » to 


convex. It is not difficult to see that a similar procedure 


< t< t] is 


may be applied, repeatedly if necessary, to any non-convex flow 


schedule. We summarize this fact in Lemma II.l. 


Lemma II.I 

Let Dy(t), to SMER ty be a non-convex piecewise linear de- 
livery function. Then there exists a convex piecewise linear 
delivery function D(t), tọ < t < tį and K < M such that 


0 
Dy (t) = Dy ft). wt e [tyrt,)- 


There should be no doubt that in the context of our problem the 


delivery function Dz (t) is preferable to D, (t). As a result 


M 


“By "delivery function” we mean, unless otherwise indicated, 
a feasible, convex, piecewise linear delivery function. Simi- 
larly, we use loosely "flow schedule" to mean a feasible, 
multicommodity piecewise constant flow schedule. 
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we conclude that the search for "optimal" delivery functions 
may be confined to the subset of convex delivery functions. 


2. Preference Relation 


Here we answer the following question: If Be 
i 


A B 


o and Dy we, to E $ tt are two delivery func- 


0 


tions, which is preferable? From our previous discussion of 
piecewise linear delivery functions (see (II.9)) we know that, 
for a given initial system state, a delivery function is 
completely specified by its corner points and delivery rates. 
We define T to be a descriptor vector for a delivery function 
Du't). 

TÉ (e 


pile ) (Gir eo, 


Pi eg Parcs: M'PM 


Deezimitıon IL.2 


A 
Given two convex, piecewise linear delivery functions Deel 
Ms 


A B B A : 
to Bee ty and Dg (E) > E, SS th we say that Dy (€) dominates 
B A B Deren BE , 
D,, (t) E e < BS and T = Tis ea ee pe ay a OC eS MC.) 
J < min(M,K), where T, denotes the j-th component of T. 


i 
The implication of Def. II.2 will become evident in the follow- 


ing examples. In Fig. II.7, DS (t) dominates Ds since 


09 (2) < a). imei, Lies, D5 (t) dominates DS (t) since 
_B A ; A B A B 
ti <t]. TOE Igea TT 9, D, (t) dominates D3 (t) since t, <t,. 


O SE ma lit Criterion 


The definition of optimal flow schedule now follows 


directly from the preference relation that we established. 
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Higg 


ek. 


Su 


A 
— Dz (t) 


==. Dale) 


N U 


té A B 
2 ty tz 


Comparison of Delivery Functions 


A 
/ o3(1) 
TA 
a 
JA 
y 05 (1) 
/ A 
EZ 3 
VA ae B 
P 2 
Zr B A 
DE al al 


Comparison of Delivery Functions 
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D(t) 


Be A 
Br A 
| m 
y D, (t) 
/ a 
Y 3 
B A E 
to $ ty t, eat, = 
Fig. 11.9. Comparison of Delivery Functions. 


Derınıtion 11.3 


We say that a delivery function Deco to me E a is optimal 


for an initial system state Q(t) TEE 
E ale A = 2,...,2M (11.11) 
i w 3 = 


El 


In words, Dy (t) issanwontimal delivery function if it is not 


dominated by any other delivery function. We shall call a flow 


schedule PU (t), to T EN which generates DAE) an optimal 


flow schedule. 
There is a technical question concerning the existence 
of Dy (t) : Can we be sure that M, the number of corner points 


Snte, because if not the delivery function will not be piecewise 


oy 








Became We defer the discussion of this point to Appendix E, 
since we need additional results before we can COBENWIEH SE: 

Our definition of optimality is based on a particular 
preference relation, which we believe to be a natural one within 
the framework of the optimal delivery problem. Obviously, other 
Bir zerenee relations (and thus optimality criteria), based on 
the delivery function concept, can be designed. For example, 
the minimal total delay objective (see [1], [2]) may be speci- 


med, In terms of a delivery function, as 


E 


l 
min w S q; (t ) - D(t)ldt (11.12) 
{pD(t)} + (i,j)eN,) 0 
0 0 
where Q(t,) = 0, and the optimal delivery function is that one 


meh Satisfies (11.12). 

A legitimate question to be raised is: Since the optimal 
delivery function approach as well as the minimal total delay 
objective seem to have merits in the same workframe, are they 


related? To answer this question we need to introduce the con- 


e prsorzesdobal Optimality. 


Definition II.4 


* * 
We say that a delivery function Due), to sts ti is globally 


Optimal iff 


Die) me), se elle E ana der. Grae) 


U] 
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Theorem II.2 
* ® ® . x 
Fmppose there exists a globally optimal delivery fuaction Dy (E) - 
* * 
Then Dy (+) and its generating flow schedule Ey (E) solve the 


optimal celivery as well as the minimal total delay problems. 


Proof: 
Tf we denote the total amount of data stored in the network 


at time t, by Go then (11.12) may be written 


0 


nn 


D(t) = arg {min [(t,-t,)q, - f D(t)dt]} (II.14) 
K r 5 | i 0 0 


The graphic representation of (11.14) is shown in Fig. 11.10 


D(C) 
m 
0 
=a D (C) 
N x | 
* 
o Dy (t) 
un a 
to Er t = 


BETT. 10. Graphic Representation of Eg. (11.14). 
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The objective is to minimize the shaded area. Assume that 


DP (t) solves the problem and Dye (t) A Dy (t) . Then, since by 
= . * . * * 
Aer Tinition By (+) = DL (t), YE x6 [ty rt] the double shaded area 


is smaller than the shaded area. This contradicts the assump- 
pon that De (t) = Dy (t) and completes the proof with regard 
to the minimal total delay problem. 

Next, assume that De (t) solves the optimal delivery 


problem. Locking at Dae (t) from t, backwards, and arguing that 


l 

eeen it satisfies (II.11), we find that it is identi- 
x 

Gad to Dy (E) - This completes the proof with regard to the 


optimal delivery problem. 
J 


Dewzomelude from Thm. [1.2 that 1f either criterion 
produces a globally optimal D(t), then they are equivalent. 
Unfortunately, the conjecture that a globally optimal delivery 
function always exists for any multicommodity delivery problem 
is disproved by counterexample as shown in Appendix B. We 
Shall prove later, however, that a globally optimal delivery 
function does always exist when all flows have a single 
destination. 

4. Example 

We conclude this section with a simple example of opti- 
mal delivery function and its generating flow schedule. We 
use the same delivery problem as in Sec. A.4 of this chapter. 

An optimal flow schedule is shown in Fig. 11.11, were 


we use the chain flow representation. 
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10<t<15 


EAT. il. Optimal Flow Schedule 


Miera» ermal delivery function 1s plotted 1n Fig. 11.12, 


The broken line depicts the delivery function that we obtained 


D(C) 


20 


iS 





Figs, 11.12. Optimal Delivery Function 


41 





for this problem in the previous section. A methodology for 
obtaining an optimal delivery function is treated in the next 


chapter. 
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III. SOLUTION ALGORITHM 


The definition of optimal delivery function (see Def. II.3) 
suggests in itself a sequential structure of the solution al- 


gorithm. We start by solving for E the minimal time in which 


1 


all the queues can be delivered to their destinations. Then, 
keeping 5 fixed we search for a flow schedule which generates 
a delivery function with a minimal possible delivery rate in 


its first interval. We denote this rate by 02. The next step 


sito hold EN and o, fixed and search for a 


point, etc. We will refer to the problems in which the mini- 


the next corner 


mal time (or corner location) are being found as Minimal Time 


Problems (MTP, for short). We will call the other problems, 
Minimal Rate Problems (or MRP). Since there may be more than 


one corner point to consider, we will usually add an index to 
indicate which corner we are dealing with. We will show that 
both types of problems can be formulated as linear programming 
problems (LP), and as such enjoy tremendous computational 
benefits. The mathematical properties of the optimal solutions 


as well as their meaning will be the subject of Ch. IV. 


PRE FIRST CORNER POINT 
l. The First Minimal Time Problem 
Let Fy (t) be an optimal flow schedule and Dy (t) its 
Peiveryerunction. Thm. II.l assures us of the existence of 


an F(t), 0 < t < t? such that FU(t): Q(0) > Q(t3) =0. The 
relation between DY (t) and Dy (t) is pr eclured (wen ta bie, 
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BKE) 


0 0 
O De 





paige III. Ll. DI (t) versus Da 


An important consequence of the above is that our search for 


0 


tí can be confined to the subset of constant flow schedules 


al 
Bruch one corner only) {F, (t)}. 


Be inition IIl.l 
The First Minimal Time Problem is given by 


MTP (1): MENEL 


l 
Sos 
k k 
A E ESO 
lji) 1 o E 
k 
De 
k(ži + 
ty, £;,(1) > 0, vlij] e Lgr 
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J A 


K ur 

į (0), TK) < No 

Ir Vlad) € Lo 
¥(1,k) N (ELT alka) 





u 


The minimal value of the objective function of MTP (1) 
represents the minimal time by which the set of commodity 
queues Q(0) can be completely delivered to their destination 
by some constant flow schedule F(t). We summarize this fact 
with the previous observation that we need to look only within 


the class of one corner flow schedules by 


Theorem III.l 
The minimal value of the objective function of MTP (1) 


equals e 


Problem (III.la) has some quadratic constraints. To 


overcome this difficulty we use a simple transformation of 


variables. Define 
K A K -e 
en = eee le Klee Nos [i,j] € Lo (IIT) 


The transformation relies on the assumption that E 100 


viously this is the case in any problem of interest. Intro- 


ioemnae(itl.2) into (I11.1a) results in 


MTP (1): min ty 
Su ta 
k k k ae 
A Ra DE ri) N 
(gži 3 3(#i) J? > E 
k we. 
oe ee) lu. . (1) Zee, ne, Giri) 
A) a 77 
Bu wy (1) Oy) ile ey (eee 
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which is an LP. The new variable represents the total 
amount of data destined to node k that traverses link De 


during a period of length et). 


2. The First Minimal Rate Problem 


Now that we have found EN we want to find a flow 


schedule which satisfies both oy and Be 





D(t) Q(t} -e) 
0 >” 
9,4) 


F, (t) = D, (t) ee 


0 0 
Ey (4) > D, (t) = 





RO 


ine LIT. 2. Dy (t) versus D,(t) 


Let Dy (t) be an optimal delivery function and let 
Q(ty-e) be the state of the system at some time tyne, where 
U < ees From Thm. II.1 we know that there exists a con- 
stant flow schedule which can transfer the system from its 
iieeial state O(0) to Q(ty-e). If we combine this flow schedule 
With part of the optimal flow schedule in the interval (tee eae 


we obtain a two part flow schedule F(t), ERS De The 
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Besultine delivery function D,(t), Or Gis ey is shown in 


Fig. III.2. We observe that the delivery function D, (t) has 
a delivery rate, in the interval Bere, o7 (1) which is 


identical to 07. 


Br rımıtıon III.2 


The First Minimal Rate Problem is given by 


MRP(1): min J ae 
(i EN 
0 
Set. 
ee er an 
eo i eee 0 


k — 
eee CU) ecto lela Gere (IE LS) 
k(Zi) E] 1) 0 
See O E ee 
ki A ee, E 
k k k k Tr 
Í Ve 3 O 2 eee, Ke) eet 


oTo 
Co any e Such that 0 < gms t¡7t>1 


where 


CJ 


As in the case of MTP(1) also here we can Stace thar 


Theorem IITI:;2 
The minimal value of the objective function of MRP (1) equals 
o0 


I: - 
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In geonditien one an (TII 3) requires knowledge of 
$ which has not yet been determined. This dit tiem) Ss only 
of a theoretical nature, however, since in practice we can 
choose e as small as we wish. If the very small e we picked 
is too large, i.e., £ > Sa then we will miss one corner 
point of the optimal delivery function and the solution we 
obtain will be suboptimal in this sense. On the other hand, 


euere za not much to lose by overlooking corners in the optimal 


delivery function which are an infinitesimal distance apart. 


B. SUBSEQUENT CORNER POINTS 
By now it should be clear that the procedure of Section A 
in principle can be applied M-times in sequence to obtain the 


optimal delivery function. Thus, for example, the optimal 


fetmery unction Solution to MTP (m), which we denote by D? (t) 


0 
M 


2m-1 elements of the descriptor vectors (see (II.10) and Def. 


0 0 0 0 , 
II.2). These elements are Eroare rom] Em The optimal 


eermetdes with the optimal delivery function D. (t) in the first 


delivery function solution to MRP(m), which we denote by Dt) 
coincides with the optimal delivery function Du (t) in 2m des- 
criptor elements, namely Er Since the term 
"optimal delivery function" may be confused with partial solu- 
tions we will usually use "optimal delivery function of order m" 
whenever we mean a solution to MTP(m) or MRP(m) for some m, 
Im = M. 

denon. 211,3 and Fig.) bil.4 to help idenrify the 


variables we can write linear programs for solving MTP(m) and 


MRP(m) as follows: 
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l. The m-th Minimal Time Problem 
Befinition III.3 


The m-th Minimal Time Problem is given by 


MTP (m): min t 
m 

s.t. 

T k k k 

EI. We) ul) q¡(0), w(i,k) e Ny 

p=1 4(#i) 7 Gan) 

) uf. (p)- uz. (p) O (i,k) e Ng’ ¡A ‚Mm 
3(4i) 7? a) 


K rs: 
-t =) . . (m) A LE 
m ly k(#i) 15) 0 
K te ae 
C.. + ) a lo (LLE AN 
A) K (#4) T) 0 
k 0 am 
ei CD) = Ne a! e Lg» 
k (Hi) 7 = 
p= 1,2; ‚m-2 
0 k E Sr 
-t__.P__ a) y O ) u,,(m IL) 
EN, RA aa 
= Gea 
DO CT ate) g utp = Atop p =1,2,-+ 2m? 
i,keN, j(#i) *? j (#4) 
0 
a mel + “m-1 
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k u 
een 2 9, v(1,k) e N, *(1,3] e Lp, p= 1,2, ‚Mm 
for given 4° Be = A 
ele m m= |e - 
We used the following change of variables: 
USE = 
i A, ¡O 02.0,M-2 
uss (P) = x (i,k) € Nor ee Ly 
k 
and 
Os 0 o _ 2 
AE = a re p 1,2, ‚m-2 


Fig. III.3 shows the relations between the various parameters 





om MTP (m). 
0 
a 
$ ——m-2) 
m | m-l po aes 0 + 
L D(t) 
ml ae in 
yo 
a 
7% 
(m) 
/ 
1 
A 
A 
A 
ES 
0 0 0 0 0 | 
0 tu tn tel Em-2 ta E 


Fig. 111.3. The m-th Minimal Time Problem 
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Using Thm. II.1 and the same reasonirgas in Section A 


we see that 


Theorem III.3 


The minimal value of the objective function of MTP(m) equals 


0 


ee 
m 


C] 


2. The m-th Minimal Rate Problem 
Berinjtion III.4 


The m-th minimal rate problem is given by 


MRP (m) : min ) r, (m) 
(i,K)eN, 


S.C, 


k(ži) + 1) 
(Ci ton) 
r$ (p) = a Dp =a; pie 
(1,k)<Np p 
IS k — = 
CoN er Bun, == A 
r,(p)» ye ee ee ee) No Ey 9’ P 
0 0 0 
non sol Ven See nt. €, r 
where 
eo = ) T ip - i Nor *(i,k) € Nor DD, = 12 en 
1 17 ql 
ALA j (ži) 
(ITI.54) 
DIN 0 0 = 2 
e i Ser p 1.2, re 





and e is any real number such that 0 < e < t0 m 


m m+1” Bau 


III.4 shows the relations between the various parameters of 


MRP (m). 
Dit) 
0 
“i 
0 nn 
me = T 
0 — l 
cies — E A m-1) 
o Am 
m 0 
Vs — Dy(t) 
/ 
J --- D (©) 
7 
/ (ml) 
/ 
/ 
/ 
/ 
/ 
/ 
0 0 0 0 E 0 
0 4 ml Em rl SJ = 


Fig. III.4. The m-th Minimal Rate Problem 


Again, using Thm. 11.1 to prove the existence of D(t) 


as shown in Fig. III.4 we conclude that 


Theorem III.4 
The minimal value of the objective function of MRP (m) equals 


0 


Am’ = 
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C. SUMMARY OF THE SOLUTION PROCEDURE 

In Sections A and B we showed how to formulate the MTP (m) 
guceuke({m) in LP form. The solution algorithm consists of 
an iterative procedure in which we solve that pair of problems 
Ber seach new corner point of the optimal delivery function. 
We refer to this kind of procedure as Sequential Linear Optimi- 
Zation (SLO). We will see in Chapters V, VI and VII that the 
SLO methodology is a very powerful tool for solving a variety 
of complex multicommodity network and certain other problems 


as well. 


i the Algorithm 


We now define conditions for algorithm termination. 


0 


Suppose we have just solved the MTP(M), and thus found En: 


The optimal flow schedule of order M, De) is actually the 


M 


optimal flow schedule we are looking for. Suppose that we solve 


now the MRP to obtain on. It is obvious that the new flow 


schedule Fu+1 (t) will generate a delivery function Dat) which 


D(t) 


Te 
a 0 
Ay — D(t) 
/ 
/ p, (M) --- D- (t) 
1 
1 BM = 


AS. Stopping Rule 
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nee 0 ! 
moe lidentical to Du), andi tche poinesS A and Beererjcg. II1,5 


will coincide. Thus the stopping rule can be written as 


RM = 0 = stop LETS 


M 


We can summarize the solution procedure "formally" in 


the following statement. 


Algorithm: 
m = 0 
Loop 
m toe 


solve MTP (m) and MRP (m) 


0 


(M) = Pu 


IE Py then stop 


Repede 


EJ 


2. Computational Complexity 


We now propose a conjecture which we believe to be 


correct although we have not been able to prove it rigorously. 


gemsjecture, IlI.l 


Let DYE) be an optimal delivery function. Then the number 


M 


of corner points M is bounded by |N,|. 


M < Ng! rer) 


where |N,| denotes cardinality of the set Np. 


o | 


[J 


We briefly discuss this conjecture in Appendix F. 
We do not intend to comment further on the issue of 


computational complexity. There is a vast literature which 


54 





deals with efficient solutions of medium to large linear pro- 
grams, and much of it is dedicated to multicommodity struc- 
tures (see [3],[4]). The objective of this thesis is to show 
how to apply these LP methods to the solution of certain dy- 
namic network problems. 

We conclude this chapter with a simple example of a 
computer solution to optimal delivery problem. An additional 
example is provided in Appendix G. 


3. Computer-Solution Example 


2 A 3 de 
q, (0) = 85, q; (0) = 30 


T = 1, 7i, j] €< Ly 


le oe TEAN o 
q3 (0) = 15, 33 (0) = 20 q, (0) = 10, g, (0) = 50 


Pecqemiiieo. Delivery Problem 


The delivery problem in Fig. III.6 was solved using the algor renn 
presented in this chapter. As a result the optimal delivery 
function and its generating optimal flow schedule were Found» 


The optimal delivery function is shown Po A 
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21 0, =2 
19 L 
0 = 
10 Sn 
0 3 
85 0, => 
60 o 
0 10 5 20 50 55 E 


Fig. III.7. Optimal Delivery Function 


The optimal flow schedule solution is shown in Figs. 


INT. /a to e. 


Fig. III.7a. Optimal Flow Schedule for t e ORO 
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) 


1 = = E= 
q,(10) = 95 (10) - q5(10) = 40 


U 


Ib. Optimal Flow Schedule for t < (10,15] 


q} (15) = a; (15) =70 


E. 


95 (15) =35 af (15) 


De. Optimal Flow Schedule for t <« (15,207 


so 





af (20) =65 


q3 (20) 2er 


Bigatti: 7d. Optimal Flow Schedule for t e (20 50] 


— 


= (50) =5 a (50) 


eo A 


ees ae. Optimal Flow Schedule for t e (50,57.5] 


It is somewhat surprising to find that even a simple 
delivery problem, such as the one considered here, should give 
rise to what might seem a complex optimal flow schedule. It 


is one of our aims in the following chapters, to show that what 
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— nn 


appears to be a "random" looking optimal flow schedule has a 
Hoe Of well defined structure to it. In particular, we shall 


return to discuss the above example at the end of Ch. V. 
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EI ZZZPROBERTILES OF THE OPTIMAL SOLUTION 


In Ch. III we have presented a sequential linear optimi- 
zation methodology for solving the optimal delivery problem. 


The availability of potent computational tools, like the sim- 


plex method (see [5]), makes the attractiveness of linear pro- 
gramming obvious. But it is not only the computational advan- 
tage that LP enjoys. For example, the linear dependence of 


optimal solution to changes, within certain bounds, in the 

right hand side (RHS) makes the LP formulation specially appeal- 
ing from the sensitivity analysis point of view. In this 
chapter, we use this and other known results of linear pro- 
gramming to determine some of the structural properties enclosed 


within a multicommodity delivery network problem. 


ADE E TRST CORNER POINT 


l. On the Minimum Time E 


a. The Primal Problem 
For convenience we restate the First Minimal Time 


Problem (see (III.1b)) in standard LP form. 


MEE LI): min ty 
Ent, 
` k k k 
E a y a DeO a EN 
o j(ži) 2° Š > 
k a 
2 — € IV: 
iso i k (Ši) s > Sr. 0, »li,jl € Lo (IV.1) 
k ae 
ty, s,4 (1), ay, fl) PO nee, 
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The slack variable Ser gives the amount of unused capacity 


of link [i,j] over the period or II 


Example: 





q3 (0) ; a5 (0) 


Fig. IV.l. Three Node Network Delivery Problem 


Mneematrix ormulation of (1V.1) for the problem in Fig. IV.l 


is given below. 


min ty 
SE. 
econo 1 0 0 0 o =i) |e, q7 (0) 
Doc oo oo) eu a, 
2» o 1 0.00 s,,(1) q3 (0) 
a 1 0 0% a 
oo oo 10 1 04) |u,,(1)] = |a,(o) (1v.2) 
o 0 0 0 o 0 1 uz, (1) q3 (0) 
E O 100000 Bau) 0 
A 0 ob 0 
moro 1 0 0 0 0 1 1 a 0 
ug (2) | 
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A 2 3 1 
x” (1) P7872 (1) S23 (1) 85) (1) a) (1), 875 (1) - 45, (1) 
3 
A a 


a 


b. The Dual Problem 


The dual linear program to (IV.l) can be written 


[6] as follows: 


DMTP (1) : max D 
(ee 
0 
See. 
- ) MSI al 
a 7 
Tj; O) O e (ES) 


> 


k K 
a; (1) a a, (1) a Mig) 


IA 


O SIK) e N 


Sa E 5 


Ea 


where 


i 
> 


and 
A ADO O) 


is a vector of dual variables. The vector )(1) has n(n-1) 


components corresponding to the nín-1) conservation constraints 
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AE SL). Vector I(l) has 2 components corresponding to the 
Æ pacity constraints in (IV.1). 

Applying the complementary slackness theorem of 
linear programming (see [6], p. 77) to (IV.1) and (IV.3) we 


have 


Lemma IV.1 
MER and A(1) be optimal solutions for the primal and 


dual problems, respectively. Then »[i,j] e L\ and »(i,k) < N 


0 0 


(1) af (2) > 0 + 7,0 = ee 011) 
ne > 05 (1) - 1) x ao = 0 
Gl) 5, )>0 + m, = 0 (IV.4) 
T5302) <0 > siz(1) = 0 
where 
of (2) eo 


E 
Tausch [1.8.2 wesconcluded that in our network 
model, the zero state is reachable from any initial state of 
the system, provided that q% (0) ZEOT R I E No: This is 
equivalent to saying that problem (IV.1) has a finite optimal 
solution. Using now the duality theorem of linear programming 


(see [6], p. 72) we conclude that the dual problem (IV.3) has 
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also a finite optimal solution, and the corresponding values 
of the objective functions are equal. An immediate consequence 


of this argument is the following lemma. 
Lemma IV.2 
0 oo 
Let E, be the first minimal time. Then 


0 k 
ao e (0) (IV.5) 


(i,k) «Ny 


[J 


co capado ty 
In this paragraph we briefly discuss the question 
of uniqueness of the optimal dual solution and its relation 
to sensitivity analysis. 


Consider a linear program in standard (matrix) form 


LP: mine ze) x 


Ea 


where T is the cost vector, A is the matrix of coefficients, 
b is the Right Hand Side (RHS) vector of requirements and X 


fh 
is the vector of primal variables.’ And its dual: 


"the interpretation of these quantities in the Minimal Time 
Problem context can be deduced by comparing (IV.6) to (IV.2). 
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DLP: max z = A-b 


N (Ivan 
where A is the vector of dual variables. 


Berinition IV.l 

Berzazesiven requirement vector b, an LP is said to be stable 
or to be at a stable point if the optimal dual solution A is 
unique. ] 


Let the LP (IV.6) be written 


EP: min zZz = Up’ Ap oh lat Xp 
S C 
BA, + DAD = 6 (TV oe 
Ans An a 0 es 
where Xa is a vector of basic variables. We can transform 
(IV.8) into an equivalent linear program 
LP: min z = T B th + (TL -T eo) x 
Í B D B D 
SE. 
-] -] 
= IV.9 
Ap + B DX, B `b ( ) 
> 
Xgr Xp 20 j 


TR X, is a basic optimal solution (X, = 0) it implies that: 
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(1) Xp = B =b > 0, i.e., it is primal feasible 
(IV.10) 


ea 2, =i 33 D > 0, i.e., it is dual feasible 


and 


A = PEB CIV L 


is the related optimal dual solution to DLP. It should be noted 
that the primal feasibility does not depend on the cost vector 
T, and that the dual feasibility does not depend on the require- 
ment vector b. 


Let us now consider a new requirement vector b 


such that 


where 


Thus, € denotes the magnitude of the perturbation vector Ab 
and es is a unit vector in its direction. In what follows 
we will be concerned only with perturbations resulting ina 
requirement vector b for which the LP has a feasible solution. 
We call such perturbations, feasible. If Ab is a feasible 
perturbation then, since the dual feasibility (IV.10(11)) does 


not depend on the requirement vector, the new solution Xp to 


the perturbed problem will be optimal if it is primal feasible, i.e. 
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X, = B(b-Ab) = X,-B tab > 0 (1v.13) 


An immediate consequence of condition (IV.12) is 


Wat if Xp > Oe. if the optimal"primal solution is not 


degenerate, then there exists some real and positive El such 


eee Ern (IV. I3) is satisfied for all feasible directions 
i aña all 0 < € < El: On the other hand, if Xp is degenerate 
it is always possible to find some direction of perturbation 
such that (IV.13) is not satisfied unless the magnitude, e, 
e menis perturbation is zero. 

Now, if primal feasibility is satisfied for the 


new (and therefore optimal) solution X then from (IV.9) we 


B’ 


have 


A =T Y 


minz = T,B"“(b-Ab) = minz- oe Ab (IV.14) 
The change 
az & min z - min z (Iv.15) 
in the optimal value of the cost function is 
na s ME Ab (IV.16a) 


or in terms of the related optimal dual solution (IV.11) 


zu = oe AD (IV.16b) 
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Persmition IV.2 


Consider an optimal primal solution to LP and let / be the 
related optimal dual solution. We say that a feasible pertur- 
bation Ab is acceptable if the change Az in the optimal value 


Seene COSt function due to that perturbation is 


Za = ae AD 


— 


— 


paa lle ar from our discussion that if an opcimal 
primal solution happens not to be degenerate then all feasible 
directions give rise to a nonzero acceptable perturbation. We 
now show that the same property holds for a weaker condition 
than primal solution non-degeneracy, namely stability. 

Let Ab be a feasible perturbation and assume that 
(IV.13) is not satisfied unless e = 0. This also implies that 
Xp is degenerate, i.e. there is at least one basic variable at 
zero level. If we force the magnitude of the perturbation to 
be somewhat larger than zero this causes one or more of the 
zero level basic variables to become negative. To satisfy 


AN 


primal feasibility, a new optimal basis B has to be found by 


exchanging the basic negative variable(s) with appropriate! 
non basic variable(s). The optimal value of the cost function 


for the perturbed problem can be written as 


1 


min z = BI Cea l (EY see) 


Ab? 


"this is usually performed by a Dual Simplex pivot step(s) 
(see [7]). 
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Mee let e go to zero, then in the limit we have 


az) = min z (LV 283) 
€>0 


from which we conclude that the optimal perturbed basis B 
is also an optimal basis for the original (unperturbed) LP. 
If we assume that the unperturbed optimal solution 


to LP was stable (unique optimal dual), then we have 


eae 2 (1V.19) 


where A is the related optimal dual solution to the original 


DLP. As a consequence of (IV.19) we state 


Theorem IV.l 
At a stable point, there exists a real and positive Er such 


that all the feasible perturbations Ab, where 


and 
I Fr Ex 
are acceptable. 


The application of Thm. IV.1 to the optimal delivery problem, 


aac particular to MTP(1) results in 


corolla IV.1 


Art a stable point, let Ed be the first minimal time and let 


AQ = {age} be a small change in the data queues sizes. Then 
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De a 
the change At, in the first minimal time which is caused by 


that perturbation of queues is 


MA 9; (1)4q; (IVO 
EN 
0 2 
Necessary and sufficient condition for optimal 
dual solution being unique is not obvious but partial results 
exist. The next two results (see [5], p. 144) indicate condi- 


tions for the existence of optimal stable solutions. 


Result iv. i 
If LP has at least one not degenerate optimal basic solution 


then the optimal dual solution is unique. 


C] 


Result IV.2 
If LP has a degenerate optimal solution then the optimal 


dual solution is never unique if it is not degenerate. 


paee 


— 


To make the digression from our main exposition as 
short as possible, we defer an illustration of stability in 
Optimal delivery problems to Appendix D. For the rest of our 
study we will assume, without loss of generality, the existence 
of stability whenever we consider perturbation problems. This 
allows us to evaluate the effect of a perturbation on the opti- 
mal value of a cost function without considering the perturbed 
optimal basis, which is very convenient from the mathematical 
point of view. 

Our interest in the stability concept is strongly 


motivated by yet another property that we find very useful for 
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analysis purposes. Consider the LP in (IV.6), and in particu- 


lar let 
AS (US) (IV S2) 


where S is a vector of slack variables. We say that a con- 
Aann 13 Critical if its slack variable is zero in al 
optimal primal solutions. In the next section we show that 
1f an optimal solution to LP is stable, then all the critical 
constraints can be uniquely identified with the help of the 
related optimal dual solution. It is again worth noting, how- 
ever, that there is no loss in generality associated with the 
stability assumption, since in Ch. V we present an algorithm 
to identify critical constraints even when the optimal primal 
solution is unstable. 
am Critical Sets 

In [8] Ros showed that the optimal solution to the 
static flow routing problem imposes a certain partition of links 
and commodities into hierarchically related sets. Our study 
of the optimal delivery problem (= dynamic flow routing) re- 
veals the existence of a similar structure. We obtain a some- 
what more general result which involves critical sets of links 
and commodities as well as critical flow rates. In hindsight, 
the existence of a more general structure is not surprising 


in view of the higher complexity of dynamic versus static flow 


routing problems. Here we present results related to the criti- 
eal sets of links and commodities. The critical flow structure 
is considered in Section A.2. Extension of these results to 


al 





subsequent corners of the optimal flow schedule will be 
Seuared im Section B. 


PemStareswitm definition of the critical set L,. 


BeeinleLlon LV. 3 
Piece Staa le pont, a link [i,j] € Lo belongs to the set Ly 


iets Mz (Ed) eo 


Mn 
kos 


Mie links in E, have the property that they must be saturated 


for the whole period og in any flow schedule P(t), 


0. < t < tj such that F,(t): Q(0) + Q(t,) = 0. To see that this 
is true we need to show that the links in L, are saturated in 


1 
all optimal primal solutions to MTP(1), and that no other links 
are Saturated in at least one optimal solution. The first 
fact follows from (IV.4(11)). Now suppose there is an optimal 
solution for which there is some link [i,j] e Lo such that 
Sij‘) = 0 and wa a = 0. Applying the theorem of strong 
complementary slackness (see [9], p. 54) to (1V.1) and (IV.3) 


we have 


Lemma IV.3 
Given a pair of primal and dual programs with feasible solu- 
tions, there exists at least one pair of optimal solutions 


ESAS and A(1) for which: 


k ey eee 
(2) De za, a = 9, (2) ol) 
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(i) $,,0) O Mig) = 0 (TV2) 


DI O e 5 (1) = 0 


27) 1J 


where 


Iœ 


of (1) 
Since by the assumption of stability the optimal dual solution 
1s unique, Lemma IV.3(ii) implies an existence of a different 
optimal primal solution for which the link [i,j] is not saturated 
any more. 

In passing we note that the optimal solutions re- 
ferred to in Lemma IV.3 are not necessarily basic solutions. 
Also, Lemma IV.3 differs from Lemma IV.1 in that the properties 
(1) and (11) are implied here in both directions. 

To see that the links in E, are not only saturated 
in all optimal primal solutions but in any flow schedule that 
empties all data queues by time eS it suffices to refer to the 


a 


Pseof technique of Thm. II.1l. There it is shown how a flow 


schedule F,(t), 0 < t < t, such that F.(t): Q(0) + Q(t}) = 0 


can be replaced by a constant flow schedule, Fy (t) = F, 


CS En such that F(t): Or Oo) Q(ty) = 0. The construction 
used there assures us that a link [i,j] <€ Ly is saturated in 


the constant flow schedule iff it is saturated in F(t) tor 


the whole period Ei TPhistoroves our point. 


"See Appendix A. 
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Consider some commodity (i,k) = No: which uses 
links in the set L,. Beeren) be a collection of all directed 
link chains connecting node i to node k. Suppose there is 
some chain Ps oi ae such that ps N Ez = Y, i.e. the set 
E, and the chain pi have no links in common. In this case it 
would be possible to divert more of commodity (i,k) flow into 
this chain. Recall from Lemma IV.3 that there exists at least 
one optimal primal solution for which all the links [i,j] / E, 
are not saturated. This makes that flow diversion possible. 
But if a flow is diverted from the chains cutting through L,, 
some links in L, will become unsaturated which contradicts 


1 


the definition of the set L.. Thus, any commodity (i,k) which 


is using links of L, enjoys the property 


1 


a A Z, “pe < (i,k) (1V.24) 


If we denote by Ny the set of all commodities that 


use links of the set E, then (IV.24) leads to the following 


theoren. 


Theorem IV.2 
The set Ly is a disconnecting’ set for commodities in N, . 
m 
We say that a chain ps eek) ets -dGttvVe wre. che 
tlovweOrrcommodity k is non-zero on each of its links. Let us 


paCcmrEoumany weommoatty (i,k) which is using links of the set 


une follow here the terminology used by Ford & Fulkerson 
AO 
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L,, one of its active chains. From Lemma IV.1(i) we have 


erce or each link of that chain 


K K 


k 
OHD = -Tag P 4 0,12), lo) ep: (US Za) 


Cano echan substitution" we find that 


K Lal + k wo b k 
aj (1) = Tia (P) - oH = mT, , CL) 7 43 (4) + Og (1) 
=... - } Tels) (EV 26) 
[a, 8] ep% ae 
ll 
Since at a stable point, 
II LE] e Ly 
T yg (1) 

= 0, otherwise 

We conclude that commodities (1,K) < No which use the set L, 


have their optimal dual variable 95 (1) positive. We summarize 


Pils Observation in a formal definition of the set N]. 


Definition IV.4 
Aras table point, a commodity (i,k) <€ No belongs to the set 


LEE oe (1) ei, 


iL 
O 


Again, as in the case of the set Li, it can be shown 


that the set Ny is unique, i.e. the same for all flow schedules 


0 


” such that F,(t): Q(0) + Q(t,) = 0. 


Ez ([t), 0 <t <t 


15 





Let us consider again the optimal delivery func- 


Eron Dy (t) and its generating optimal flow schedule Fy (t) « 
More precisely, let us look at the first part of F(t) , namely 


Ey) + = ut E The delivery rate in this interval aie is 


the minimal delivery rate over all flow schedules PL IE) y 
UNS E < En Sueca EL At): OO > Q(t?) = Ue Seon Chee rest kts 


wA NiS section we know that the critical set of links Ly has 


to be saturated at all times t < Ones | and consequently for 


h 


al te eee The only commodities that can saturate Ly 


are those in the critical set N}. This being the case, we 


may expect that a is the minimal value of flow rate with 


which the commodities in Ny Can saturate the link set L}. We 


address this proposition in the next subsection. 


2. On the Minimal Rate os 


From the discussion of stability property in Paragraph 
(c) of the last subsection, we conclude that there exists a 
real e, e€ > 0 such that all feasible perturbations AQ = {aay} 


that satisfy the equation 


) os (1) Aas = € (oye) 
(i,k) eN) 


where Aq, Zeh) € No are acceptable. This means that 
e describes the change in the optimal value of the minimal 
time En which is caused by that perturbation. 

Memo: in Fig. IV.2 depicts the optimal delivery 
Fumetion Dr Det E (the solution to MTP(1)) and its 


: 0 
perturbed version, say e USES ty = Er 
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Fig. IV.2. The First Perturbation Problem 


The perturbed flow schedule FPP (t), ES tie delivers all 


l 
put sae, ne (Cal) e No of each of the backlogged data queues 
qs (0) to its destination. 


ne 


Suppose we can find now another flow schedule Py (t), 


tyne Ss = such that Py (t): AQ + 0. This flow schedule 


noe 


will generate a delivery function D, (t), tye Sie. E which 


1s shown by the broken line in Fig. IV.2. The value of the 


delivery function Oe Ce at time ty 


tion AQ that was chosen, and the locus of its possible values 


-€ depends on the perturba- 


iS Shown by the vertical line in Fig. IV.2. Our objective 1s 
Oaa perturbation AQ, for which the value DP (t,-e) will 
attain its maximum. This obviously is equivalent to finding 


a delivery function D(t), tyne SES E, with minimal value 
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~r 


of the delivery rate oy (1). We can formalize now the above 


construction in the following Perturbation Problem: 


PEKI): min ) Ag 
(NS EN = 
0 
Sach 
(i,K)eN 
IL 
Aq’ = { ) El; = ) a = OF ¥ (i,k) = No 
j (#i) j(ži) J 
~% = 
) € Cjjr “*li,jl e Ly (10.282 


“[i,j] € Lo 


where the positive value e is chosen small enough so that all 


AQ that satisfy (IV.27) are acceptable. 
L 


Due to stability assumption such an ¢ exists and we 


can rewrite (IV.28a) as 


Rene: Kr rf 
(i,k) eN, 
s.t. 
Ze Oo. = 
(i,k) «Ny 
Mea) ez, =, = (0, Vik) e Ny © (AV. 288) 
zn > 
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) E = a A - Lo 
k(#i) + J 
eR u: 
rir fi; me O07 (1,k) € Nn aa Te Lo, 
where 
k 
BS Aq. 
k å l , 
jee N, 


Moreover, for each acceptable perturbation and in 
particular for the optimal solution to PP(1) there exists a 
feasible flow schedule er Ve tyne such that 
DP (t-e) = ) (15 (0)=4q%). We stress tnis observation 

TT i I 1 
(i,k) Na 


here, since as a result we can state the next theorem. 


Theorem IV.3 
At a stable point, the minimal value of the objective func- 
tion of PP(1) is equal to F 
C 
mediate implication of Thm. IV.3 is that the 
minimal flow rate os Can be found in a much simpler way than 
by solving MRP(1). The number of variables and constraints 
in PP(1) is considerably smaller (we are not concerned with 
flow variables corresponding to the period [0,ty-el) than in 
MRP(I). The only Limitation is that the formulation of PP(1l) 
is valid for stable points. 
A more important consequence of Thm. IV.3 is that it 


enables us to obtain a new insight into the problem by a 
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Sorerlleseuay Or (TV.28b). First, we want to indicate that 
problem (1V.26b) is feasible. As an example we can pick any 


ERO g < = and a perturbation 


k A k ; 
Aq. = er; (1), do bee Nos (v2) 
where E. (1) is the net delivery rate for the flow schedule 
0 ! 
Py (t), wil,K) € No: Then 
o in (IV.30) 
a ij , , rd E 0 e 
is a solution to PP(1). The corner point of the resulting de- 


livery function a t) is schematically indicated as point A 
in Fig. IV.2. We say that the flow pattern of the optimal 
solution to MTP(1) is a feasible solution to PP(1). 

Another simple observation results from the structure 


of PP(l). 


Lemma IV.4 


mernesopcıimar solution to PP(1) 


Ps Gt tt) O GV o 
T A — 
mnis, together with the fact that o% (1) = 0, (i,k) ¢ Ny sup- 


ports our previous observation that only commodities in Na 


play a role in the first segment of the optimal delivery function. 


> . e . 

'By flow pattern we mean the flow composition, without any 
particular reference to time, e.g. an m-piece flow schedule is 
Fonseruceted of m flow patterns. 
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The perturbation equation 


od ee 1 
(1,k)eN, > = 


has a very important interpretation: it expresses the condi- 


tion that the flows of commodities in Ny need to satisfy in 


order to saturate the set Ly + This interpretation can be de- 
duced by the following argument. The flow shcedule Py (t), 


0 0 
Bee 


of PP(1) must saturate Ly with commodities from Nie 


is a part of two segment flow schedule which delivers all the 


whose flow pattern is obtained from the solution 
Since it 


0 i 
queues by t Now, if we look carefully at the constraints 


1* 
of (IV.28b) there is nothing there, beside the perturbation 


Feuers stchat Can account for this saturation property. 


Theorem IV.4 


At a stable point, let ) (1) be the optimal dual solution to 


MPT(1). Then a feasible flow pattern F, of commodities in 
N} saturates the set E, aie 
AA 
(1, k)eN, a 


The next topic we want to consider is the lower bound 
on os. Suppose that we remove the feasibility constraints from 


(IV.28b) and end up with 


> min ) E 
o, | | 
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ae = 1l (DSZ 


U 


ro 
me 1S obvious that 97 < N and thus can serve as a lower 


bound. It also is easy to see that for this simple LP 


Oy = AE (LV an 
where 

o (1) = max a 

max 2 


(irk) N 


Theorem IV.5 


At a stable point, the minimal rate os satisfies 


os O ge (TV34) 
where 
g le max BR, 
Is (i,K)eN, 
# 


A related result states 


Lenman IVS 
At a stable point, if the optimal delivery function of 


order one, Dy (t), Oe eS = satisfies 
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(Vaa 


then it is the optimal solution to the delivery problem. 


Proof: 


SUppose that condition (1V.35) is satisfied. By Thm. IV.4 


we must have that (in this case N} = No? 
en a 1 (IV. 36) 
(ITK) eN 
0 
o e iad (IV.35) into (IV.36) results in 
a (1V.37) 
(i,kyeN, == 


Eussztenerlertet hand side of (IV.37) is equal, by definition, to 
the total flow rate p, (1). Comparing with (Iv.34) we must conclude 


that D) (t) is already the optimal delivery function and 


o7 (1) = 97: 


E) 


An interesting case occurs when 
o(1) = o __ (1), *(i,k) e ÑN (IV. 38) 
i max d A it 


and 


The delivery rate of the flow schedule solution to MTP(1), 


ae may be viewed as composed of two constituents. 
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p, (1) = O(N, ) > o (Ny Zn CL Vso) 


where o9(A) denotes the delivery rate due to commodities in the 
set A. Using Thm. IV.4 and substituting (IV.38) into the 


perturbation equation we obtain 


O(N, ) = lo (1) (IV.40) 


Since all the backlogs are delivered by the time =. we conclude 


that 


) cece) (IV.41) 


We now show that in this particular case there exists 


Seselmeien to PP(l) in which 07 attains its lower bound, namely 


1/0 max (1) : Fig. IV.3 will help the reader to follow the 
argument. 
D(C) 
i i U — 
TR — = 
ANAL > q; (0) 
F(t) 0 (0) meus ty (ik) aN, 


E Special Case of the First Perturbation Problem 
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Let us take the flow schedule solution to MTP(1), 


er DIE < = and increase the flows of all the commodi- 


ties in the set NaN] . This is possible since none of these 
commodities passes through Ly, and thus all use only unsaturated 


links. As a result all the queues of commodities not in N 


1 
will be fully delivered prior to ae say by some time el-e. 
Let us denote this new segment of a flow schedule by EE 
DS E < Se as shown in Fig. IV.3. In the remaining time 


interval er. only the commodities in N] will continue 
to flow with a total rate p, (1) = 1/0 „„(l). This is denoted 
by ur tyre es EN PG e IV. Comparing this value 
Or oy (1) with the lower bound on oF (IV.34) makes our point. 
Two observations are appropriate here. First, in this 
particular case it is not necessary at all to solve the PP(]l) 


since the minimal rate value is known to be 07 = LO ee L) 


a 
ahead of time. Second, there is a slight change in the notion 
of the next time problem. In MTP(2) we will be looking ex- 
Pigejtly for the minimal time = by which all the queues of 
commodities not in N, can be fully delivered to their destina- 
tions, given that queues of commodities in N, are delivered by 
the time oe These two observations have a considerable im- 
pact in the case of single destination networks where the 


Special case we consider here turns out to be the general case. 


ERZEREEZSEESND CORNER POINT 
We now begin a study parallel to that carried out in Section 


A, and define and describe the same type of concepts and 
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results that were derived for the first corner of the optimal 
delivery function. Since most results will be totally analo- 
gous we will not dwell on their proofs, unless the reasoning 
behind a proof is very much different from that in Section A. 


1. On the Minimal Time = 


The Second Minimal Time Problem can be written (see 


Ber. Ill.3) 1n standard LP form as 


MTP (2): min t, 
Seo. 
; ( ) ne (p) - ) i (p)) = OD “(i k) e N 
p=1 3(%i1) Y 3(%1) 3* : - 
k k k 
Fu. (2) - or rn ele ES 
3(#i) +? 3 (fi) J? 3 } 
ur. (p) = ) us. (p)- as (p) = 0, v(1,x) e Ny, 
a a 
peal, 2 
k = | 
A Be (ll) =. 0, ¥{i, 3) SL (IVa?) 
E Co, 1) j 
k a 
-t.c,., + ) A e A O L 
lij k:(i,k)en, 7? 1J ; 
-t,0,, + wa 2) = 0, *li,j] « Ly 
1J k:li,k)eN, 1) 
k oy / 
EG A ..(2 A L 
2) A Sea ) Lin ) El: 


86 








Be Ta e aC SO 
LE (i,kyen, j(#i) © ji) d 
0 
O), Ue. (p) (1) 20, e E 
etme 5 Eye +2 2 Oe vital © Le 


and ae 


for given t? 1 


1 


a: 
— 


In the formulation of (IV.42) we have incorporated some of the 
results we derived earlier in this section. All the flows in 
the first time interval caves are due to commodities in the 


set N}. This is so because any commodity (i,k) Z N} ceases 


to flow (its queue is completely delivered), in an optimal 


seluclen, prior to time EN Also, there is no need for slacx 


TE 
variables im Capacity constraints for links in Li, since we 


know that they must be saturated for all t e Cr: We note 


that the third constraint in (IV.42) accounts for no intermedi- 


ate data storage (as (p) denotes the surplus variable for this 


constraint, (i,k) e Nyy p=1,2). In Fig. IV.4 we indicate 
the basic relationship between the various parameters of 


iZ Js 


2 


2 O) R2 LY T2) 0 I)o, (2) ) 


C 
be the vector of dual variables for the dual problem to (IV.42). 


The vector )(2) has |N,| components corresponding to the first 


o 


ee ey Constraints of (1V.42). The vectors S(1) and 


(2), have |N components each and they correspond to the 


1| 
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Fig. IV.4. The Second Minimal Time Problem 


"no intermediate storage" constraints for both epochs (p = 1,2 
Pirie stead constraint of (IV.42)). Vectors (1) and I(2) 
have | Ly | components each and corresvond to the capacity con- 
straints in both time periods. Finally, g (1) and g (2) 
represent the dual variables corresponding to the rate and the 


time constraints, respectively. 


Lemma IV.6 a 


Peotmae meand (2) be optimal solutions for the primal and 


dual problems, respectively. Then ¥[i,j] e Ly and *(i,k) e No 


"We will indicate in the brackets the analgous result in 
Section A. 
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E ZO Rm. (2) = 0 if [i,j] Z Ly (IV) 


15) 1) 
Tl) < 0 > s;;(2) = Orr a L, 
ve egnconserained .ın sign if [i,j] < L. 
where 
o% (2) 20 
O 


Application of the duality theorm of linear programming results in 
Tenma IV. 7 (1V.2) 


Let E be the second minimal time. Then 


N O 


B k k 0 E 
= 5 0 (219,0 + 9, (2)t7 (IV.44) 
, E 0 


E 0 l 
It is worth noting that t, depends on the queue sizes and on 


2 
t9, 
Keersenosc difficult to see that at a stable point 
9, (2) < 0: Suppose that we resolve MTP (2) while we let 
t? > Eee. Now it is possible to deliver some small amount 
of each commodity in this additional interval and we expect 


> to decrease as a result by some Ben. Recalling that at a 


stable point, (2) relates the perturbation to the change in 


SE 


es tive tunetion, we conclude that o,(2) < 0. 


€ 
The discussion of stability in Section A.l.c applies 


wıicheur change to MTP(2) and we obtain: 


89 





Mente (Corollary 1V.1) 

At a stable point let > be the second minimal time and let 
(AQ, At y) be a small change in the data backlogs sizes and in 
the minimal total delivery time =. The corresponding change 


in the second minimal time At) is given by 


At, = y oe (2) Ags + 9, (2) At, (IV.45) 
(i, K) No 


N O 


m 
— 


Following the discussion in Section A.l.d we define 


the critical set Lo. 


Berinmelon IV.5 (IV.3) 
ea swaple point, link [i,j] e Lo belongs to the set L, 
Eicher if 11,3] e L, or E Ly and Dadas e Us 


=j 
ud 


The links in L_ have the property that they are saturated for 


2 
the whole period [0,t5] in any flow schedule Fy(t), 0 < t < es 


such that 


(1) F(t): OO) Y 
(IV. 46) 


TE) gli) = 

We note that the necessity to break the definition of 
L, into two exclusive cases results from the fact that 1,5(2) 
is in general not restricted in sign (see Lemma IV.6(1i)), 


Bar links [1,3] « L}. Thus, a simple statement like 
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Ra] € L, 2,12) < 0, will not be correct. With the same 


precaution we define the set N, as follows. 


Desınition IV.6 (IV.4) 


At a stable point, a commodity (i,k) e Ny belongs to the 
set N, either if (i,k) « N, or if (i,k) £ N, and 0%(2) > 0. 


Ü 
The commodities in N, have the property that they must 
flow through the set L, and saturate it during the interval 


5 | Or all Ey (t) such that (IV.46) is satisfied. Using 


the same arguments as those leading to (IV.24) we have 


Theorem IV.6 (IV.2) 
The set L, is a disconnecting set for commođities in L.. 
Ü 
A false impression may result from our discussion, 


namely that each corner point in the optimal delivery function 


implies a new pair of critical sets. This is true only for 
the first corner point. It is possible that, for example, 
L, = L. and correspondingly N, = N}. This would be the case 


when the corner point occurs because the commodities in Ny 


can not maintain the minimal flow rate os any longer (backwards 


in time). At this point a new rate 09 will be computed without 
change in either L, or N}. This brings us into the discussion 
of 05. 


2. On Minimal Rate 05 


E : À 0 
Our objective here is to study the properties of 9, 


and its interpretation. As before, we will rely on the stability 
assumption to simplify analysis. 
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D(t) 0 


tL ae 
Də (t) a 
o 
/ 
(b 
max D5 (t) 
DR D, (t) 
0 
E DS (t5-<) 
e € = + € ze 
0 
= = E 


Fig. IV. 5a, The Second Perturbation Problem 


Fig. IV.5adepicts the optimal delivery function er 


wee tS En (solution to MTP(2)) and its perturbed version 


1 
Be), Ve tyne. The perturbation equation derived in 
Section A (IV.27) is not applicable here. We want to find a 


perturbation AQ such that De will be maximal but at the 


22 
same time we must preserve previous results, i.e. the minimal 
rate o and its duration aoe The perturbed delivery function 
De) in Fig. IV.5a has the required properties. It is described 


by a generalized perturbation equation 


c = Fo (2) Age + en, (2) (IV. 47a) 
(MEN 
0 
which may be also written as 
I y k k an 
— San rt (IV.47b) 
7. 12) (i,k)en, 7 = 
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As in PP(1), we require from the flow schedule P (t), 
eo e =, to satisfy F, (t) s Bos 40. Tne broken line in 
Fig. IV.4 denotes the delivery function D, (t) which is gener- 
ated by F) (t). Now, if we permute the segments denoted by 
(5) and (a), by permuting the corresponding flow schedules in 


time (this is always possible), we obtain the desired struc- 


Prize as Shown in Fig. IV.5b. 





0 t3- € 40 4° E 


Fig. IV.5b. The Second Perturbation Problem (Permuted) 


Since maximizing Tea is equivalent to minimizing 


the delivery rate 07 (1) we can finally state the PP(2). 


I: 





PP (2): 


min A 
Cal 0 
Sis 
il ç K “k 
Se Cre) = 1 
male) (i,k)eN, + 2 
(IV. 48) 
) a = >) Er - E A AS 
j (Xi) 3 (Hi) ? 
` er 
R Ca’ Yin] L 
IA I EN 
gd EK a: 
ee gg O Lo: ¥(i,k) € No 


E 
Peele estabilicy assumption, all perturbations AQ 
that satisfy the perturbation equation (IV.47b) are acceptable, 
and for each perturbation there exists a perturbed delivery 
unction DAE) of the form discussed. Due to this fact we 


can state the following. 


maconem IV. 7 (IV.3) 


At a stable point, the minimal value of the objective function 


tOr PE(2) is equal to 05. 
O 


The analysis of PP(2) leads to exactly the same results 
as for PP(1), aside from the slight modifications introduced 
by the factor 1/1-0, (2) in the generalized perturbation equa- 


tion. We state those results here for completeness: 
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From the structure of PP(1) we have 


Lemma IV.9 (IV.4) 


Meade OOELmal solution to PP(2) 


Ber oo (IV.49) 


This, together with the fact that 0512) = 0, *(i,k) £ N, 
supports our previous observation that only commodities in 
N, play a role in the second segment of the optimal delivery 
eimetion, 

The interpretation of the generalized perturbation equa- 
tion is analogous to that in PP(1). Observing that this is the 
St condition in the statement of PP(2) that can account for 


Mercer enat commodities in N, must saturate L we have 


2 2! 


Mectem Iv. (LV. 4) 
t a stable point, let ) (2) be the optimal dual solution to 
MTP (2). Then a feasible flcw pattern F, of commodities in No, 
Saturates the set L» SEE 


i 


R k k 
=o tr a (2), = 1 


aSK 
O 


The lower bound on 05 has basically the same form as 


the lower bound on oe 
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moesrem IV.10 (1V.5) 


At a stable point, the minimal rate os satisfies 


1-0, (2) 
0 t 
a EA (IV.50) 
max 
where 
a = max lo, (2)) 
(1I KIEN 


2 


we can be certain that the flow schedule solution 
0 l 
D,(t), IAS E, generates the optimal delivery function if 


“er. Lemma IV.5) 


(IV. 3D) 


a 
N 
Il 
a 


i el, “(IS K) € N, 

The special case that we discussed at the end of the 
last section applies here as well. As we indicated there we 
will discuss it in detail in the single destination networks 


case in Chapter V. 


C. SAMPLE PROBLEM 

From previous sections it is apparent how everything studied 
so far generalizes and applies to any corner point. From a 
conceptual point of view only two corner points have to be 
studied, since corner three, or any subsequent corner presents 
no significant difference with respect to the second corner 


point. We believe that it will be more enlightening to present 
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a detailed study of a sample problem than repeat the same 


ideas again. 


Example: 


q (0)=4q 





Fig. IV.6. Sample Delivery Problem 


As the first step in solving the delivery problem in Fig. 
IV.6 we wish to find the composition of the first critical set 


Ni: Obviously, Ny can be only one of the following: 


A a 1(2,3)+, (iii) 1(3,2)}, Kiv) {(1,4),0372)}, 

ea, (vi) —1(1,2),12,3) 1, (vii) 1(1,1) (2,3),(3, 205% 
In what follows we consider explicitly each one of the above 
possibilities. 

(1) -- For this case commodity (1,4) must have all its chains 


(here, there is only one) going through Ly, which implies that 


Ac pone or the links [1,3], [3,2] or [2,4] must be in L}. 
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But this would imply that at least one of the commodities (2,3) 
era ,2) uses links in Ly and thus belongs to N}. This contras 


BMets (1). 


(ii) -- For this case commodity (2,3) must have both of its 


available chains going through L which implies that at least 


A 
one link in each of the following pairs, {[2,1],[1,3]} and 


me 4 |, (4,3); must be in L fers NOt Gitriculte to see thar 


1* 
Ia is saturated so is [1,3]. Similarly, if link 
[4,3] is saturated so is [2,4]. This implies that commodity 
(1,4) uses links in Ly and thus belongs to N. nie cons 
erets (ii). 

(iii) -- In principle the same type of argument applies 
here; we write in shortened notation: 


i 2 


(iv) - In this case L, = {{3,2]}. Suppose that commodity 
eeepeiseusing link [3,2] with some flow rate a, and conse- 
quently commodity (3,2) is using that link (remaining capacity) 


with flow rate l-a. Then it is true (remember that at this 


Eregerwe are Solving for a constant flow schedule) that 








4 2 
¿0 o q4 (0) _ q, (0) 
ine; OL l-a 
or equivalently 
T q l-a 
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N = eS 2 > SPA e L > (1,4) € Ni == Ny = Be 





We can eliminate the flow variable a by using the law of 


proportions, namely = = = = == Such that 
0 2a +q _ 
Se a + l-a — = 
Hew, we scan dlso solve for a, to obtain a = 2/3. The remain- 


mie apacity On links [1,3] and [2,4] is thus 1/3 (actually, 
mess enam 1/3, say e, E UE e Commodi cy (2,3) Z Ny and 
che solution to MTP(1) it must flow only through unsaturated 


links). We find that the time required to deliver the queue 


Stecommoaity (2,3) is 


Emrch ceoneradceits our assumption in (iv). 


(v) -- We have already seen in (iii) that if commodity 
D2) e Ny then also commodity (1,4) e N), which contradicts 
nv). 

Exercising our advantage over the reader in knowing the 


oe let us study possibility (vii) prior to (vi). 


(vii) -- This case is equivalent to (iv) if we let e = 0 
mow commodity (2,3) e Ny? in our discussion there. Then it 


Must be true (for an optimal solution to MPP(1)) that 


37 (0) q5 (0) q3 (0) 


E a Bea 








or equivalently, 
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Ae ac ee ese 
q 1-0 2 (l-a) 


HO 
NO 


where the last equality can never be satisfied for non-zero a. 


We conclude that (vii) cannot be accepted as correct. 


(vi) -- For this case Ly = {({1,3],(2,4]} (the remaining 
commodity (3,2) uses the unsaturated link [3,2] only, as 
required). The chain flow decomposition of the constant flow 


Sebecdule solutien to MTP(]l) is shown in Fig. IV.7. 











2310) 
8 < l-a 
Brest. 7. Chain Flow Decomposition 
It must be true then that 
4 3 2 
0 q7 (0) q, (0) G3 (0) A 
A ae 
where 
Bear en 
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If we consider in (IV.52) only those elements that have the 


variable a associated with them, we have (using the law of 


proportions) 
4 3 
BEN, ,;, 1. 
il T 2a + 2 - 2a >, qı os 2 q>10) LEVE 33) 


For the queues sizes we have selected, 


t) = 4q (IV.54) 
Using (IV.52) and (IV.54) we can evaluate au, 

ne = = >, (IV.55) 
and 8 

gs = g@ = Fr (IV.56) 
which is less than l-a = 5 (as desired). 


Now that we have established a Ny and Ly we turn to 


calculate ae the minimal delivery rate in the first interval 


of the optimal delivery function. Recalling that (see (IV.5)) 


0 k k 
Ey = ) o. (1)q,(0)» 


; Bu 
(i,k) «N, 


Ecco n par lnea to (IV.53) results in 


4 a E, all 2 
g (1) = 1, o-(l) = 5 (IV.37) 


JEON 





rO Obtain os we have to solve PP(1) (IV.28b), 


PP (1): min = 
(1,k)eN, 
SE. 
p 2 o, (IV.58) 
KIN. 
L 
and 


tri}, (LK) € Ny 1s a feasible set of delivery rates. 


L 


n ag TuS); problem (IV.58) can be solved by 


inspection yielding 


SN. ~g 


3 = 
(14,1) i (170); (IV -39 


Deelwiszchewn in Fig. IV.8. 





DRDS.) Solution to the First Minimal Time Problem 
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Thus, we conclude that 


The fact that the solution to PP(1) is unique, considerably 
simplifies the formulation of MTP(2). Let Q(t.) denote the 
E From our study we know 


system state at time t ESA 


2 2 
Bares 1ets a Constant flow schedule F(t), 0 < t < E, 
SUN ENaC F(t): O(0) > Q (to). This fact can be expressed by 
the chain flow decomposition in Fig. IV.7 if we substitute 
Q(0)-Q(t,) for Q(0). Also, we cannot be sure any more (nor is 
it required) that B < l-a, since it is possible that N, = Ny U (37 20i 


Because of the uniqueness property of the optimal solution to 


PP(1), we have that the components of Q(t.) are 
0 0 i 2 
(tı-t,)07, (eK) = (1,4) 
qe (t5) T (IV.60) 


0% otherwise. 


Since qe (t5) £ TO *(i,k) e N, then in particular 


0 
4 
qs (t,) = (es -t,)0) = q, (9), (LUCO 
or eauivalently 
t = q7 (0) ug (TV 76 ly 
2 — È it 


Equation (1V.61b) presents a lower bound on the minimal value 


0 
oe t5, 1.e. t > 2q. 
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Let us check whether it is possible that E) = 2q. This 


is equivalent to asking whether the queues of commodities (2,3) 
and (3,2) can be delivered within the time interval [0,2g}. 
Since commodity (2,3) and (3,2) have no common links in their 
respective chain flow decomposition we can assign a flow rate 
DO commodity (2,3) and a flow rate of 1 to commodity (3,2). 
It is easy to see that it will take = 2q and + = q units 


of time to deliver the respective queues to their destinations. 


The complete optimal flow schedule is now shown in Fig. IV.9. 








AO. Optimal Flow Schedule 
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The optimal delivery function is shown in ELO 0 


D(t) 

7q 
0 
Al =1 

A 0 

q 

os 
om 0% or 
0 t,=q t, =2q t, =%q E 


Fig. IV.10. Optimal Delivery Function. 


E 


A remaining problem, which does not arise in the foregoing 
example, concerns the possibility of loops existing in the 
£low solution. Evidently, such loops cannot affect the opti- 
mality of the delivery function, but nonetheless their exis- 
tence is not aesthetically pleasing. Two comments are in 
pude MED ESst, such loops cannot appear if the input traffic 
between all pairs of nodes is non-zero. Second, given the 
time and rate parameters of the optimum delivery function, all 
loops can be eliminated by solving the last flow problem 
again, but this time with the objective of maximizing the sum 


of the link slack variables. 
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V. SINGLE DESTINATION NETWORKS 


By a Single Destination Network (SDN, for short) we mean 
Eıeeza,lzehe data traffic in a network is destined to a single 
node. Without loss of generality we will assume that node to 


be n, n£ V. Accordingly, for present purposes we redefine 


No SE Or all nodes such that node i can communicate 


to node n, 1 X*X n. 


In general, we will sımplify the notation by dropping out the 
Aestanetion indication, which is implicitly understood to be n. 
Naturally, all the results that were derived for multi- 
commodity case apply to SDN. Their generality though, tends 
to hide some of the unique properties of SDN which we study 


here. 


Bess THPFSETRSIESCORNER POINT 


The single destination variant of MTP(1) is given by 


MTP (1): min Cy 
Since. 
al 
a = 2 U = Ol ZN 
j(i) Íl 3(%i) 31 i 9 
Se + DE + E O lle Lo (Val) 


O e Lo 
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where we have used the previously defined change of variables 


e aas] € Lp - 
The dual linear programming form to MTP(1) can be formu- 
lated exactly as in the multicommodity case. The same state- 
ment applies to the discussion of stability and to the defini- 
tions and properties of the sets E, and N. Let us assume, 
ENS pOint, that the optimal solution to MTP(1) is stable. 
Hence, the sets Ly and Ny are uniquely determined by MTP(1). 
Later on in this chapter, when we discuss the solution algorithm 
for SDN, we will relax this assumption. 
By Thm. IV.2, the set L, is a disconnecting set for nodes 


EN Nie i.e. every chain that connects any of the nodes I = Ny 


Heras sei nacton n has at least one of its links in the set 


L We shall see in the next few paragraphs that the maximal 


1 
flow rate IN) with which data can be delivered from the 
set N, to the destination node n is given by the Max-flow 


Min-cut theorem (see [10], p. 11] 


p ON = max ) Ese ES IN wu) 
De {F} ien, = 
where CS (Ny) is the value of a minimal cut-set, separating 
the set N, from node n. 
MOS important to realize that in the SDN case the mult#- 


commodity delivery problem turns into single commodity problem 


“Discussion of the differences between single commodity 
and multicommodity can be found in [9]. 
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which makes the notion of minimal cut-set meaning un TO 


prove (V.2), consider an optimal flow solution to MTP(1). A 
typical source node i, i « Ny Bw heown in Pig. Volta). The 
initial queue q, (0) is diminished with a rate r, (1) (net 
delivery rate of data from node i) such that 3, (67) = 06) Baro 


(b) in describes an equivalent setup made up of a 
virtual node v connected to node i by an infinite capacity 
link [v,i]. There is a constant data input with rate r, (1) 


to node v. 


r= r, (1) 





ALA Source Node in an Optimal Solution to MTP(1) 


If we extend the model to all nodes in the set N the result 


lee 
is as shown in Fig. V.2. 


Recall now that the critical set Ny consists of all nodes 


ied eE No which determine the minimal time oe Am implication 
of this characterization is that it is impossible to increase 


any of the initial queues (while not changing the others) 
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Fig. V.2. The Flow Pattern of F(t), (ee E 
without causing a to increase. Equivalently no delivery rate 
r¡ (1), 1 € Ny can be increased, even momentarily, without 


decreasing some of the other delivery rates. 

To show that the flow pattern originating in node v and 
zerminating In node n, as shown in Fig. V. 2, is maximal it 
suffices to demonstrate that there is no flow augmenting path 
(see [10], p. 12] from node v to node n. The existence of 
such a flow augmenting path would manifest itself in an allowa- 
Ba Eee As wOeerlow rate on exactly one ol the links [v,i], 
for some i < N.. But this would be equivalent to an increase 


I 


in exactly one delivery rate r; (1), for some i e N, (without 


i 
changing any of the others) which is just what we have shown 


to be impossible. This brings us to the conclusion: 
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Theorem V.l 


0 | 
Let F(t), USE $ ES be an optimal flow schedule solution 


EOmMNTP (1). Then 


) (13) 


As all the queues are reduced to zero by the time ES we 


can write 


t = eng (V.4) 


But it is also true that (see Lemma IV.2) 


= ER o,(1)q, (0) (US) 
a: 


m © 


Comparison of (V.4) and (V.5) raises a question about the 
functional relation between the set of optimal dual variables 
and the value of the minimal cut-set CS(N,). The next lemma 


answers this question. 


Lemma V.1l 
At a stable point, let INT) be an optimal dual solution to 


MTP (1). Then 


G, (1) = 1/C5(N,)» wi e€ Ni UE 
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Proof: 


Suppose that not all o,(1), i < N, are equal, and let a and 
b be a pair of sources in Ny Such that cu) > „m. Define 
a new delivery problem for which 


q, (0) -Aa, ar i =a 


q; (0) = qp (0) +Ab, if i=b (Ver) 
3,100), otherwise 
and 
“9, (1)Aa + 0, (1)Ab = 0 (V.8) 
where 


D > Aan 0 


At a stable point, we can always find a perturbation (Aa,Ab) 


for which (V.8) is satisfied and the perturbation is acceptable. 


Beerssmetdirrieult to see that condition (V.8) implies the 


fact (see Corollary IV.1) that 


ze (V.9) 


where E is the first minimal time for the new problem. The 
total delivery rate of sources in N, is, for the new flow 


schedule, 


TEL 





p (Ny) = o ___—_—_ CU 20a) 


Dien Can be also written as 


P (N) = P(N.) + € (V.10b) 
where 
2 a AB Aa 
¿O 
I 
we already know that o (Ny) = EN a and from (V.10b) and 


the assumption that g (1) > oa) we have that e > 0. Thus 


O ame = est) (v.11) 


een of course is impossible. This contradicts our initial 
assumption about an existence of unequal dual variables. As 


a consequence we may rewrite (V.5) as 


0 
A Ol (V.12) 
L<eN, 


and comparison to (V.4) completes the proof. 


C] 


It is interesting to observe in consequence that in SDN, 
Ebezfirse minimal time En is always equally sensitive to changes 


gn any of the queues in the set Ny > 


VEZ 





o 


u 


One of the significant properties of the set L, is that it 


1 
must be saturated by flows originating in the set N}, through- 
out the interval BR The total rate of the saturating 


flows was found to be maximal and thus equal to CSÍN Pe 


1) 
Mea te obvious that this flow rate cannot drop from its 
maximal value, even momentarily, since this would cause some 
of the data coming from sources in Ny to be delivered later 
than = This observation holds for all flow schedules, as 
long as they terminate by time ty. We conclude that no seg- 


ment of an optimal delivery function may have a slope less than 


Lemma V.2 


Let er Verne Er be an optimal delivery function in SDN. 


Then 


= CS(N,), meel ML (Ve) 
We will show now that (V.13) is satisfied with the equality 
for the first segment of an optimal delivery function. Suppose 
that in the optimal solution to MTP(1) we find that Ny = No: 
The total delivery rate of ON Ve = is maximal and 
eaual to CS (No? ° It can be easily seen with the help of the 
sketch in Fig. V.3 that the existence of a delivery function 
(shown by the broken line) that dominates oe would imply a 
delivery rate o" (No) > CS (No) which is of course impossible. 
We must conclude that Da OR a cS = is the optimal delivery 


minceron. This proves our claim for the case of N, = No: 
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2 O (Np) CS(N,) 
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Fig. V.3. Delivery Function with Maximal Flow Rate 


In general we may expect in a solution to MTP (1) that 

Ny = No: In this case it is possible to increase all the de- 
livery rates of sources in NoN] since they use (by definition) 
only unsaturated links. Hence, all the data backlogged in these 
nodes will be completely delivered prior to time e, say at 
some time e In the remaining interval only data from the 
set Ny will continue to flow in the network, with the rate of 
CS(N,). This simple construction proves the existence of a 


two part (e > 0) flow schedule with a delivery rate of CS (Ny) 


in the interval (EI-e,t1]. 


Corollary Ma2 


Let Dy (t) be the optimal delivery function in SDN. Then 


DZ 
pP] = CS (N, ) (V.14) 
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The proof leading to Corollary V.2 deserves some addi- 
tional discussion. Let us consider an optimization problem 
in which the objective is to minimize the total delivery time 
of data backlogged in NoN]? while keeping the delivery time 
of data queued in N, at time ae This is a slightly more 
formal statement of the construction method we used to prove 
Berellaryz’.2, Since the solution to this problem will auto- 
matically satisfy (V.14) we may use this new formulation as 
a substitute for our previous formulation, which required both 
MRP (1) and MTP (2) in order to obtain the value of ES. The 


substitute optimization problem can be written as follows. 


2 
Sch 
2 
Be - unap) = DN = l,2 
DR ij P? or sai q; 1, P 
Y) u..(p) - ea Cp) AA Nie = eee 
te) ee AG S i 
zur) i = 0), *i < N -N (ALIEN 
a ai | did oe 
0 
ESC ¡y - u < t1S;3 
E - rn el, ve Lo 
ir t, = 0, ¥ [1,3] € Los p = I 


1.355 





where we have used the notation 


iœ 


a) Det) 


0 
ij (Ey - E 


E 


ED 


E 


GE) ae) = t, E 


The formulation (V.15) depends on the knowledge of the set 


: 0 
Ny (and the time ti). 


we indicated before, we will later show that the stability 


and hence on stability of MTP(1). As 


requirement is not necessary. 

We now show nct only that MRP(1) is not needed in SDN, but 
that the optimization problem in (V.15) can be formulated in 
a much more efficient way (with regard to the number of varia- 
bles and the number of constraints). The basic idea is that 
the flow pattern of data delivered from the set Ny in the 
period eed can be made identical to its pattern in the 
interval Kr In other words, let Ps A Den active 


link chain used by some source i, i e Ny with a rate Eher 


es: (0,65). Then the same chain can be used to forward data 
from source i with tne same rate ze Ir *t € el, In 


order to see it we need the following result. 


Theorem V.2 
re etan Optimal flow schedule Coe Dre En 


for which 





u T 0 i 
r, (t) = , ¥t [0,t3) and “vi e Ny (V.16) 


Oy 


116 





Proo®: 

Before we start with the proof, we remind the reader that 
r,(t) 1s the net delivery rate of data from node i at time t 
(see (TL.2)). 

Let i and j be any two nodes such that ie N, and ee NoN}. 
a a E T e the sets of all directed chains con- 
necting node i and node j, respectively, to destination node 
n. Select any active chain pP; eo Ti and let x De the head 


node of the first link in that chain that belongs to L, (there 


1 
is at least one such link Since i e Nj). Let PG E Ps denote 
mic upasetaleenain connecting node i to node x. Our claim is 
Enac 


aa Ps 1S active then Ds N Ps = Y, “py TEN) (Ve) 


This claim can be easily verified by referring to Fig. V.4. 





3 ENJ=N, 


CU A ELrustration for Theorem V.2 


27 





Suppose that there is some active chain J Sa) such 


endet 


N Pi = Ra | (Velo) 
Since In, Ny + the flow on this chain rj [pily cannot go 
ehreugn Ly, so that the chain must have the form (see Fig. 


V.4) 


Ps A al (v.19) 


But this would imply the existence of a chain (not necessarily 


active) 


re ar... , O (72,20) 


which violates the fact that Ly is a disconnecting set for all 
nodes in N}. Wemmust® conclude that (V.17) is true. 

Now, consider the optimal solution to MTP(1) and let Oe 
be an active chain in that solution, for some ie Nj - Our 
proof of claim (V.17) indicates that no data flow from any 
of the sources in NoN} 


the partial chain De and this is true for all partial chains 


may ever (in any flow schedule) use 


Geeenisetyee for all 1 < Ny + In view of this observation we 
may require, without loss of generality, that an optimal flow 
schedule will have the same chain flow structure ene ler 


vieN,) Euer ro, Solution to MTP(l). This in turn implies 


NA cc nMI? () the delivery rates ot all the sources, 


TES 





and in particular those in Ny, have the form 


q,; (9) 0 
EACC) = pe le Noe ee [0,t,]. 





(V.2]) 


This concludes our proof. 


Thm. V.2 does not mean that the solution to MTP (2) has 


no effect on the chain flow structure of sources in N Due 


1% 

only that any effect must occur beyond the critical set L,- 

Tores is schematically indicated by point d in Fig. V.4. 
Suppose next that we have solved MTP (2), so that the chain 


flows of sources in N, incorporate interactions with the chain 


1 
mows originating in NaN]. Since in the period r only 
the sources in N, are active we may require, without loss of 
generality, that their chain structure in that interval, will 


remain the same as in the interval Dre) This result leads 


us to a new formulation of MTP(2), which is described next. 


Br ÜBSEOUENT CORNER POINTS 


1’ En and armed with the results of the last section 


we can formulate the Second Minimal Time Problem as 


Given N 


MPT (2): min t 





) __° ) u. 
En j(ži) 77 j(ži) J1 
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N 
Q 
F 
E 

LA 


rt 
to 

G 
|v 


07 E € L 


where we have used the transformation of variables 


> 


IE er u a) ic Ly - 
It should be noted that the number of variables as well as 
Constraints is exactly the same here as in MTP(1) (cf. V.1). 
At a stable point, the optimal dual solution to MTP(2) may 
be used to identify the sets L, and N, (see Def. IV.4). More- 
over, using similar arguments to those used in the proof of 


Thm. V.l, an analogous result can be obtained for sources in 


N,. Wefstate without proof: 


Theorem V.3 (Thm. V. 1) 
Let e) ES = be an optimal flow schedule solution 


COMTE (210% Then 


(i) ) r,(2) =CS(N,) 
1eN, 

(ii) O A y LID ES) 
LeN, ieN, 


[—' 
ha 


Proposition (ii) above was already proved in Thm. V.2, but we 
include it here for completeness. 
Similarly, it can be shown (by analogy to Lemma V.1) 


that at a stable point, the optimal dual variables corresponding 
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Co Sources in N -N} are positive and equal, i.e. 


o; (2) Z C eea e NN (200 


The functional relation of One) to the value of the minimal 


cutset CS (N, ) can be determined as follows. Using Thm. V.3, 


we have 
t2Cs(N,) + ) (0) 
2 1 di 
0 o 
ee Va 


or equivalently 


0 JU 
A ae NT ) a; (0) (TV. 2245) 
2 CS(N,) CS(N,) ieN,-N, i 


But we also have (from Corollary IV.1 and the form of RHS of 


22) that 


So 0 2 ) 120). (ve25) 
2 max ieN5-N, al 


Comparing (V.25) and (V.24b) we conclude that 


Lemma V.3 


At a stable point, let )(2) be an optimal dual solution to 


MP2). Then 


i 


ole) = CS(N,)-CS(N 


22, 
il 


A E No-Ny 


pel 





Next, using reasoning completely parallel to that in the 
proof of Lemma V.2 and its corollary we arrive at a similar 


result, which we state formally as: 


Lemma V.4 
Let Dy (t) be an optimal delivery function in SDN. Then 
05 = CS (N,) (22) 
g 

Tewm acon Of this result and its constructive proof 
is that we can proceed directly to solve a new version of 
MTP(3), without bothering to solve MRP(2). This new version 
of MTP(3) may be phrased as follows:-- "Minimize the total 
delivery time of queues in the set NoNo? such that all queues 
in N.-N 


Dae: 
are delivered by the time Et In order to show that this 


are delivered by the time E and all queues in Na 


MTP (3), like MTP(2), can be formally stated in an efficient 


way, we need to prove a result parallelto Thm. V.2. 


Theorem V.4 


There exists an optimal flow schedule Mor Oe = E for 








which 
| E 0 | 
(1) r,(t) = Ta et e [0, t] and vi < Ny 
El 
| jaa ee [0,t>] 
IR E : 
(1i) r(t) = 2 yi e€ NoN} (V. 238) 


O , otherwise 
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ESO OLE 
Proposition (i) is a restatement of Thm. V.2 and is included 
for completness. 


Enesgzoescrris in essence identical to that of Thm. V.2. 





1 
a 
$3 ; 
J eN,-N, 
KENN, 
Big. V.5. Illustration for Thm. V.4 
Dumas Claim as that for all sources jj, j € No-Ny, any active 


subchain = remains free from any interaction with flows 


Originating in the set N.-N 


g7No- Suppose that a flow coming from 


some node k, ke NoNo interferes with a flow coming from some 
node j, j € N,-Ny as lar o). Since the flow from k cannot 
use the set L, (and in particular L,-L,) Joy? (6 Gage) g\hie esol es he 


must use the bypassing chain which goes through node c. But 
this would imply that the node j has a chain outside the set 
L, which is impossible since L, is a disconnecting set for 


nodes in N3. This basically completes the major arguments of 


the proof. The remaining details are as in the proof of Thm. V.2. 


ro 
i 


Caw 


T29 





By now the reader has no doubt surmised (correctly) that 
all the results derived up to here can be extended to subsequent 
corners of the optimal delivery function. We therefore conclude 
that the general character of the optimal delivery function 


is as illusted in Fig. V.6, where 





A oprimal Delivery Function in Single Destination 
Networks 
0 = CS(N_) 1,2... 
m mas 
aD 
m~Aa A es, n= i, 
E CS(N,) ten. 2 
e = (V.29a) 
m 
Į i = 
u ze. m = 2,37. eee 
mn no 
Me mel 
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and 


mee eneral orm of MTP (m), m = 1,2,...,M is given by 





MTP (m): min tn 
Se 
q; (0) r S Rem 
-t a L uj; T u S ‚wieso 
oo a Ga) PRA 
k “k=l" 
KS 2737 ‚M 
lies, “= r E 
A ea) 2” a 0 m-l 
ne + De OA A qe nn (Veco 
= es a E == j 3 € = 
ee eis a ee oy 51 
En’ i G = ur w[i,j} = Lo: 
[ 
N 0 | 
where rt ea are given, and 


2 
'There is no need for slack variables since the links 


ea EL 2 0 
ae L are saturated for all t e LOrE E 
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SEE STEOBSZL OPTIMALITY 

We stated previously that one of the distinctions between 
a multicommodity delivery problem and a delivery problem in 
SDN is exhibited in the fact that every optimal delivery 
aseo DNS also globally optimal. This is not true 
for multicommodity case, as the counter example in Appendix 


B indicates. 


Theorem V.5 
Let Dy, (t) Dero pti mal delivery function in SDN. Then it 


is also globally optimal. 


Proof: 
Assume to the contrary that there exists some other delivery 


function D(t), Ta as. En os. atole! 


A 


1 0 f 1 
D(t > Da (€ PALOL Some t" e 10,t 


0 


1) - (SON 


With the help of the optimal delivery function let us find m, 


0 
m+] 


all the data stored in nodes of the set N 


D Mlzsweh that t' e (t Se Let us mark now 


-N_ so it will be 
Omen 
distinguishable from the data stored in the nodes of the set 


N Before applying the flow schedule F(t), OF -< Teme En which 


An 


generates the delivery function D(t), let us place an observer 
demno e e duty 1s to count how much marked Qn and unmarked 


on data is delivered to node n, up to time t'. 


Sd 


M6 





feos clear that 


oe) q, (0). ve 
1EN Nm 


O) 
A 


The observer cannot count more marked data than there was 


iaa iy. in the network. Also, 
Qa < CS(N_) +t (Vos) 


Since no more than CS(N_) of unmarked data per unit time can 


reach the node n at any given time. Thus 


D(t") 


IA 


ES NE (V. 34) 


] -N 
1eNo De 


But the right hand side of (V.34) is exactly the value of Dr, (4) 
An ten contradicts (V.30) and completes the proof. 
O 

The proof of Thm. V.5 is not dependent anywhere on the fact 
that the flow schedule P(t) is feasible in the narrow sense 
(see Appendix A), i.e. does not allow for intermediate storage 
of data in the network. Combining this observation with the 
result of Thm. II.2 We know that the optimal delivery function 
and its generating flow schedule also solve the minimal total 
delay problem over the class of flow schedules that allow 
intermediate data storage. As such, we obtain a much simpler 


solution algorithm to that problem than the one described in [11] 


"The work of Shats and Segall seems to be the only known 
open loop solution to the minimal total delay problem in SDN. 
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D. SOLUTION ALGORITHM FOR SINGLE DESTINATION NETWORKS 

Up to now we have established that the optimal delivery 
function and its generating flow schedule can be obtained by 
solving a sequence of specialized MTP's. In addition, the 
size of each of the problems is limited to 22+1 variables and 
n+£-1 constraints independently of which corner point we are 
solving for. The only condition that we assumed, for the above 
Rome Crue, 1S Ehat a solution to MTP(m), m = 1,2,...,M iden- 
tifies uniquely the critical set N a’ or equivalently, is stable. 
When we say the set N n’ we really have in mind the nodes in 


since the nodes in N were supposedly identified at 


m-l 


a nl” 


prevíous corners. 


Consider an optimal dual solution )(1) for MTP(1). Define 
Ny = set {i} of all nodes i « N, such that o, (1) > 0, 

and 
Ly = set {[i,j]} of all links [i,j] e L, such that a 


A useful interpretation of the sets Ni and de can be obtained 


il AL 
with the help of the following lemma. 
Lemma V.4 
Let )(1) be an optimal dual solution to MTP(1). Then 
o, (1) >0 + ie N, Cas) 
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Proof: 
Let i, 9i € No be some node for which a, (1) > ONE FOR (EV. 10) 


we have (for SDN, k = n) 


> = ) n Tag fh)» (V. 36) 
la ,8]ep; 


no é' : 
where p; is any active chain connecting source i to the destina- 


tion n. Since Tag ae. Ly we conclude that there exists 


at least one link in the chain DER forynleh TaB SS From che 


slackness theorem (see Lemma IV.l(ii)) we know that this particu- 


lar link will be saturated (zero slack variable) in all optimal 


primal solutions, and hence belongs to the critical set L}. Any 
source node using this link must belong to the set N; - We con- 
clude that ie N}. 

0 


It is clear now that the set Ny consists of all the members 
er Ny that were uniquely identified by the optimal solution to 


MPA TI). 


N 2 Ny, (VES) 
and similarly 
i = 
ee ty 


where the equality holds if the solution to MTP(1) is stable. 


It should be observed also that the set Ny 


Saneoeveastly deduced from application of (IV.5) to MTP (1): 


cannot be empty. This 


129 





O (V. 38) 


V 
© 
4 
= 
m 


Since Es Orand q; (0) No there must be at least one 


mode 1, ie No Rei wien o; (1) > 0. This also implies that 


Re tar least one link (a,8] < Lo such that Tag P => 0% 


Up to this point we have no way of identifying the remain- 
ing members of Ni: But let us try to proceed with the solution 


algorithm in spite of this fact. We will use the set Ny in- 


stead of N as a result we obtain a slightly different formu- 


E 
Sons for MTP(2). 





MnP(2): min E, 
Si Es 
Gono) 
-t, = at l üa T u teen Dee Ny 
en a A 
5 | 1 
) A ae nec 0) aN (V.39) 
IR, 1 
Sd E m 
A = 0, »[li,j] e 4 
2 7 iJ , , 1 
sE GC oo Fun +H Ss Olia] £ T 
27 1) ia 19) Á í 1 
it to, “mij = 0, *[1,3) E Lo 
O 


T « s 
‘In this equation we don't need slack variables since we 
know that all the links in L} have always to be saturated. 
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Suppose now that actually Ny c N}. Then it is impossible 


to deliver all the queues in the set Ny -NÌ prior to time En 


while keeping at the same time the delivery rates of all the 


1l 


nodes in N” at the value 4,10) /t). We therefore conclude that 


Lemma V.5 
0 l ; 
Let t> be the optimal value of the cost function for problem 


wo). > Then 


(V.40) 
This does not necessarily seem like progress until we con- 


sider the optimal dual solution to (V.39). Using again (V.5) 


for our problem here, we obtain 


0 0 : 
) h o.(2)q,(0) (v.4l) 


1fNy 
which implies (similarly to (V.38)) that there must be at least 


one node i, i / Ny 


that this node belongs to N 


such that g, (2) > 0. We want to show now 
1: Let us pick some active chain 
for that node. If the chain passes through nee we are done, 

since this implies that i < N}. If the chain does not go through 


it 


1 then we can repeat the argument used to prove Lemma V.4. 


Defining 
u i 
Ne = Ny u {iljo,(2) > 0 and i ¢ Ny} (V-42) 


dl E 


and 
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} (AD 


E: e un 1 
E ({i,jlin,, < 0 and [i,j] # Ly 


we can solve (V.39) again, this time using NT and = instead 


of Ny and L 


cussion we are assured that after a number of iterations k 


7 respectively. In view of the precedding dis- 
1’ 
k) $ IN] | the whole set N, (and L,) will be identified and we 
may proceed to solve the original MTP(3). 

The general idea behind the solution procedure and how it 
applies to subsequent corner points should now be clear. The 
only consequence of instability in any of the corner points 
is to increase the number of iterations needed to reach the 
next corner. If we let K denote the total number of iterations, 


where 
M 
K = ) X, (IV.44) 
then we obviously have 


Moe Koc | N (V.45) 


o! 
where M is the number of corners in the optimal delivery 


muneecLon in SDN. 


E. REMARK ON MULTICOMMODITY FLOW SCHEDULES 
In this subsection we briefly discuss a class of multi- 
commodity problems for which the optimal dual variables satisfy, 


at every corner, the following relation: 
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A a (V.46) 


In particular, it can be shown that the data flow rate with 


respect to the set N. is maximal in the period Gr i.e. 


k E 0 
A en oor), e e (V.47) 
4 -Zm 


The delivery of all data queues corresponding to the set 
No "Nm can be accomplished prior to ar say by tyne, Since by 


definition the flows of commodities in N Bi only use unsaturated 


0 


ea 
links. Thus the generalized perturbation equation must be 


Satisfied by flows of commodities in the set Nn (see Thm. IV.9) 


max y eos = 1, vt e (tome, to] (V.48) 
(i,k)en > 


Using (V.47) in (V.48) we have 


wane ee C=“ (v.49) 


As a conseguence, the optimal delivery function for this class 


of multicommodity problems is characterized by: 


| 0 ee 
(1) Pm a ~ 07 m) > 
max 


at 
IA more detailed argument can be found in Chapter IV.A.2. 


1 > 





l 
(1,k)<N_-N 
En. 0 f m -1 
lat) ie = —  z- ==  M= 2,3,...,M 
= ee eee meat: 
and 
k 
i q; (0) 
O (MS EN 
E = l 
l P max (N) 


Moreover, it is easy to see that the proof of global opti- 
mality for SDN applies without change to the multicommodity 


case considered here if we use instead of CS (Na 


P nax im? 
since the two are not equal in general for the multicommodity 
case. 

It turns out that computer solution example which we con- 
srderec in Chapter III,3 falls into this category of "SDN 
like" multicommodity problems, i.e. multicommodity problems 


Aca ica ne "optimal" solution is also globally optimal. 


For convenience we restate that delivery problem. 


2 u 3 u 
qj (0) = 85, q,(0) = 30 


ea wurd, et 


ij 0 
| 
q3 (0) = 15, q3 (0) = 10, 
2 (0) = 20 >(0) = 50 
G3 22 


A Delivery Problem cr Chapter III.C.3. 
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The critical sets for this problem are (from the computer 


SO Iie er) ©: 
I De 2 emer; 3) } oe ee 8 Ny 
= E E No, E) y u No, (Va S 
and N, = mr) o Ny. 
NOE dr Eicuit o see from Fig. V.7 that the maximal 
flow rates with respect to the critical sets are: 
ze. Sec ne) a. Pmax (N3) Ze 
(Vas2ca) 
men =>, and nase) E 
Using (V.50) we obtain the optimal corner times 
Bes NA D0 se 
ism a, t4 = 3757 >0 
(25525) 
O20. DES 
SS S 75770 
Ino. 
and ts = 535 = 10 


Comparison with Fig. III.7 will immediately reveal that (V.52) 

is a correct description of the optimal delivery function. 
This completes our promise at the end of Chapter III, to 

show that what looks like a random flow schedule has indeed a 


Mes ate ure, and thus simplicity, to it. 
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F. SAMPLE PROBLEM 

We conclude this chapter on single destination networks 
with a short study of a sample delivery problem. This gives 
us an Opportunity to illustrate some of the special proper- 
ties that are characteristic of SDN. 

In se Shats and Segall presented an interesting algorithm 
for solving minimal total delay problems in SDN. We have shown 
(see Thm. V.5) that the optimal delivery function is also 
globally optimal in SDN and thus (cf. Thm. II.2) serves as an 
optimal solution to the minimal total delay problem. Unlike 
in [li], the sizes of linear programs that are required in our 
solution procedure are independent of the number of corners of 
the optimal delivery function. This makes our algorithm ade- 
quate to solve large network problems, a task which could not 
be handled by the authors there. We believe that there is also 
another advantage to our methodology, namely the additional in- 
SIGNE T it provides., 

We have adopted one of the computer solution examples from 
(11, p. 73] as our sample problem. We intend to show that by 
using concepts introduced here it can be solved, with very little 
Sort, actually by inspection only. 

As the first step in solving the delivery problem in Fig. 
ea sho find the composition of the first critical set 


Nae ebyaGgusiy, N 


1 can be only one of the following: 


1 


"The study in this reference aroused our interest in SDN, 
and inspired many of the results obtained in this chapter. 
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capacities are indicated 
on links 


Fig. V.8. Single Destination Delivery Problem 


Der), (iii) 13), (iv) {1,2}, 


A A o5y A 272 


Let us now consider all the possibilities for which node 
Ale Se Ny (i secs em a vi) and (yii)). As a result of our 
study we know that the chain flows originating in N, must satur= 
ate CS (N, ) amsc-heropeimal solution for all t e (oe. An 
immediate consequence of this statement is that node 2 has no 
available chains, during that period, to send its data to the 
destination n. This leaves us with possibilities (ii), (iii) 


and (v). By using exactly the same reasoning we can discard 


possibility (111). 


17 





Now, suppose that N, = {2,3}. Since the flow originating 


I 


in Ny must saturate CS(N,), one of the following chain flow 


decompositions must be part of the optimal flow schedule: 





Fig. V.8a. Chain Flow Decomposition for Ny ZN 23) 


Or 


q, (0) => 





Fig. V.8b. Alternate Chain Flow Decomposition for N, = Ma 


Then, it also must ke true that 


| 


05 
ae) - Ita 


> 
| 
e 


Ox. 


Ib Sys 





DI ES SE 4 
E T R 
Bersrmersfirstet Case we have a solution for a, a = l À and for 
the second case, B = -15. Both of course are infeasible and 
we conclude that Ny ar 


We are ready to construct the optimal flow schedule. The 
flows from node 2 must saturate CS(2), which is equal to 2. 
The chain flow decomposition which achieves this is unique and 


iseonewn in Fig. V.9. 





Been o tamal Plow Schedule for Source (2). 


We are left essentially with a delivery problem shown in 
Pige- VLO. 

The final solution should be obvious to the eye at this 
point. For completeness though, we continue with our formal 
exposition. 


Suppose now that 


fe, 3) ne NS => ON 








avallable capacities are 
indicated on links 


Fig. V.10. Delivery Problem for Sources in N,-N 


Orr 


this would imply the chain flow decomposition which is shown 


MELO. Vall, 


q, (0) =2 
g,(0) =4 
2 
3-al 
Omer 
Hee ieee Chain Flow Decomposition for NN = low 
Also, 


qe. 
2 2+0 


ger 
3-0 


140 





O dS to a = - 1/5, which is unacceptable. The only 


doming possibility is shown in Fig. V.12. 


SS 





Uo] 


Biever Gpeimal Flow Schedule for Sources (1) and (3) 


The solution presented in Fig. V.12 is equivalent to the 


following statements: 


(i) In er Ny 
(33) N, = {1} vu N, 
where 
N, = (1) 


The resulting optimal delivery function is shown in Fig. 


V. 13, where 
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8 E > 
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een Optimal Delivery Function 
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Be a I 
Pee mes? 3)=-CS (2) 3 
3 CS5(2,3,1)-CS(2,3) 
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VI, APPLICATION TO STOCHASTIC DELIVERY PROBLEMS 


We have studied in depth so far the optimal delivery prob- 
lem. The mode that we have used in our discussion is based on 
the assumption that during the period of interest the data 
input rate is identically Zero. We find it convenient, espec- 
lally in view of the forthcoming discussion, to refer to that 
class of delivery problems as "deterministic." 

We now focus on stochastic delivery problems. Here the 
data input rates are assumed to be governed by some stochastic 
process. In this new framework, the time necessary to empty 
a data queue (deliver its contents to their destinations) is 
no longer a deterministic value but a random variable. In [12] 
Meemsuggested Use of the expected delivery time as a performance 
measure for dynamic routing. We demonstrate that the Sequential 
Linear Optimization (SLO) methodology, which we used to solve 
the deterministic case, can be applied to the stochastic 
case with Yee's performance measure. Before we do so though, 
we briefly summarize the most common stochastic routing model 
(13] and indicate why the new performance measure seems to be 


advantageous. 


i 
A. BACKGROUND 
Beer onnmtes or departure is the original (cf. Chapter II) 


2-link, n-node model for a communication network. Data entering 


Tin this section we closely follow the discussion in [13]. 
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the network from external sources forms a Poisson process with 
a rate of a (messages per second) for those messages entering 
the network at node 1 and destined for node k. All messages 
are assumed to have lengths that are drawn independently from 
an exponential distribution with mean l/u (bits). The combined 
effect of finite link capacities and random fluctuations in the 
actual arrival rate of messages to the network causes queueing 
delays. In order to accommodate these queues we assume that 
all nodes in the network have unlimited storace capacity. At 
any given time, the state of the network (or its congestion) 

is described by the set of data queues. 

With each link [i,j] e€ Ly we associate, in addition to its 
capacity, ae a queue one of all messages waiting to be trans- 
mitted over that link. The routing of messages (flow pattern) 
in the network is accomplished by determining for each node 


Merde rion Of the incoming traffic, for each commodity (i.e. 


destination), is to be directed to which link queue. 


other commodities 
pr Nerea, 


\ 

x 
\ 
| 


conmodity K, 
to other links 


ee 
C 


E to ne j 
1J 





Fig. VI.l. Schematic Representation of Node-Link 
Queueing Model 
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The data rate conservation equations for any node i, i< V 


= 


eiampemwartten (with the help of Fig. VI.1) as 


Q Er) - J 2 = a, *(i,k) eN (V1.1) 
ine) j (ži) 


A 
o> 
Lt 


“{i,jl e L “k (VE 


f.. +a. 0 
i 


=, sic ic 

We are now faced with analysis of a network of queues. A 
Similar problem was studied by Jackson [14], and he was able 
to establish that an imbedded queueing and serving facility 
offered a solution identical to the same facility acting inde- 
pendently from the network, but with Poisson arrivals at a 
rate offered by the network. In order to apply this remarkable 
result here it is necessary to assume that every message, once 
it arrives at its intended queue aE has its length randomly 
selected anew from an exponential distribution with mean 1/⁄/u. 
This destroys the dependence between interarrival and service 
(transmission) times. This assumption was studied extensively 
Dotes], with the conclusion that the so-called 
"independence assumption", albeit rigorously unjustified, 
leads in practice to useful results. 

With the independence assumption, we see that any link 
[i,j] « L, is now representable as an MIM|ı" queue with Poisson 


0 


> 
‘Detailed study of M|M|l queues can be found in any basic 
text on queueing theory. 
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grrıvalıs Of rate fij =) 2: and exponential service rate 
k 


with mean l/ue, For notational simplicity we will assume 


j . 
tnat alí capacities in the network include the factor u, and 
thus the mean service (transmission) time of a message over 


mas (apa) is ee In what follows we assume steady-state 


M 








operation of the M l system, and thus require that 


fij k Cij (VIL) 


For the queueing model described above the average delay 


T of a message passing through the network is given by [13] 


a j C.-f.,. 
[i jle Lg 3) 1) 
where 
A k 
a z ) a; 
(i, K) <No 
and 
= y K 
1J k(fi) + 


The most common statement of the routing problem involves 
minimization of the average delay T as the objective function, 
eRbseer tor constraint (VI.1). Various solution methods for 
this non-linear problem have been presented (e.g. [16] through 
[201) in past years. In [8] a different approach leading to 
a linear programming formulation was taken, namely a satura- 


eion ratio i; is defined for each link [i,j] « Lor and 
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the worst of them is minimized. This procedure is then 
iterated until all saturation ratios have been minimized. 
When implementing any of the above methods in actual routing 
control it is necessary to update the estimates of input data 
rates more or less often (depending on the nature of external 
data sources), and to recalculate the routing variables in 
order to adapt to possible changes. It should be noted that 
in (VI.4) the flow rates (ff), and respectively the routing 
variables (aja) depend only on the input rates {as} and not 
on the actual congestion (9,3! in the network at the update 
time. It is reasonable that inclusion of global congestion 
information in determination of the routing variables should 
improve the adaptivity properties of any routing methodology. 
In particular we follow [12] in suggesting the expected time 
needed to empty a gueueing system as a practical performance 


measure which uses congestion information. 


Er EL HEZEXREEETED TIME TO EMPTY A QUEUEING SYSTEM 
In this section we will derive (following Yee) the formula 


Por ce the expected time needed to empty, for the first time 


1 
an M|M|1 system. 

The first question in this respect that we wish to answer 
is: If a message arrives to an empty system, how long will 
it take, on average, before the system becomes empty again? 


Some thought will show that this is exactly the expected length 


of a "busy period" t,, which is known to be for MIM|1l 


A ee CVT 25) 
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where a and 1/c are the arrival rate and the expected service 
time, respectively. 

Now, suppose we look at a queueing system and find (q-1) 
messages awaiting transmission, and one message being trans- 
mitted. Let us agree to put any message that arrives from now 
on in buffer (b) (see Fig. VI.2). Also, we will always empty 


BRO) PDuUtrer prior to servicing the (a) buffer. 


new arrivals 


Fe 
|) ul (a) 
2 


Fig. VI.2. Queueing System 


It is clear that the next message in buffer (a) must wait, 
before beginning transmission, the length of a busy period 
(VI.5). (We note that the change in queueing discipline that 
we introduced does not affect the distribution of busy/idle 


Periods for MimM|1.) Following the same argument it is clear 
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that the system becomes empty for the first time after 


waiting, on average, 


t = gt, = $ (VO) 
We summarize this part of our discussion in the following 
lemma. 


Lemma VI.l 
The expected time ty needed to empty, for the first time, 


an M|M|1 system with q messages is 


where a and 1/c are the arrival rate and expected service time, 


respectively. 
E 
Consider a commodity (i,k) e No which is characterized at 
a given moment, say t = 0, by its queue qs (0) and its arrival 


rate a$, Suppose we dedicate to this commodity a part of the 
capacity resources of the network in such a way that they can 
support a constant flow rate = or zeommoadrey (i,k) from node 1 
to node k. From our previous discussions we know that it is 
impossible to assign to commodity (i,k) more capacity than 
CS(i,k) the value of a minimal cutset separating node i from 


node k, i.e. 


ff < CS(i,k) (VI.7) 


MSs situation is depicted schematically in Fig. VI.3. 
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i £* < CS(i,k) k 


Era. 2 0ueueing System tor Commodity (i,k) 


K 
At this point we can write that the expected time t, needed to 
1 


empty, for the first time, a queue of commodity (i,k) is 


k 
A q4 (0) 
E; = S K (VIB) 
tes a 
L jE 


We would like to consider the same construction simultan- 
eously for all commodities in the network. Since the capacity 


resources are limited the following constraints must be satisfied: 


(i) Dinkscapaceity constraint 
Dre se (VI.9) 
el) Ze ee 
where g5 denotes the part of link capacity c.. that is dedi- 


I. 1) 
cated for commodity (i,x) use. 


(ii) capacity assignment continuity 


k | 
ES E, IR K) e Ny (Vilas 
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The constraint preserves a constant total assignment, say 


a Semenetwomn Capacity to commodity (i,k) along all paths 


from node i to node k. Substituting (VI.10) into (VI.8) we 


obtain 
k 
k EOS 
te, 1(i,k) ce N (VEAN 
a MS E. as : 


ea) Ze 


From all feasible capacity assignments that satisfy (VI.9) 
and (VI.10) we ask for one which has the minimal largest ex- 
pected delivery time. This is the same Min-Max criterion we 
used in the formulation of the First Minimal (deterministic) 
Time Problem. The formulation of the corresponding First 
feminale expected Time Problem (METP(1)) and subsequent optimi- 


zation problems are the subject of the next section. 


C. SEQUENTIAL LINEAR OPTIMIZATION FORMULATION 


i 
With every feasible capacity assignment, (E) we asso- 


e e e ptor vector T = (t,,t,,...,t,), M< EN of 


M 


distinct expected times to empty the queues. We assume, with- 


out loss of generality that the components of T are ordered 


such that ti Pen iC i > 3. Now, given two feasible capacity 
l A B 
assignments pa and E we say that ES dominates D IEE E, < Es 


dol for some 3, 3 < mia (1, Mp) 


A The definition of optimal capacity assign- 


ment now follows directly. 
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Beminition VEI.l 
We say that a capacity assignment is optimal for a given 


network state la; (0) ) and arrival rates {at}, i Re N 


O 
BoE © 
OF | 0,0 os a 
eS so ft |tyrta,----t) 4}, m=1,2,...,m (VI 12 
ee 


Our objective is to show that the optimal capacity assignment 
can be obtained by solving an appropriate deterministic optimal 
delivery problem with constant rate data inputs. 

We now present the statement of the First Minimal Expected 
Time Problem (METP(1), for short), in which the largest expected 


queue delivery time is minimized. 


METP (1): min ty 
a (0) 
Zw ee O E AO 
ME A Jee 
es rigen (VI.13) 
ke 1J — 1) 0 
E ae Seo * (17) 1 se & * (i,k) £ N 
ID ij A í 0’ A 


Qq 


mee eemal olution to problem (VI.1l13) is a function of the 
initial network congestion {qe (0) } and the expected arrival 
rates Car This open-loop solution can be implemented, at 


least in principle, as closed-loop routing control by continuously 
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Beeteeubating it in time, with the current network state (con- 
gestion) as initial condition for each problem. It should be 
pointed out that we do neither imply nor believe that this kind 
of implementation is possible unless the frequency with which 
a new solution is recomputed can be adjusted to match the 
actual transmission and computation capabilities of a network. 
In spite of this fact, we will assume in the sequel that 
it is possible to recompute the capacity assignment with every 
new message arrival to the network. As a result we obtain a 
theoretical model which provides some new insight into what 
Howenererttrect Of including congestion information, in addition 
to that of arrival rates, on dynamic routing strategies. 
Following the assumption above we specify the set Noe in 
this chapter, to include those commodities (i,k) for which 
a: (0) zer doing so we do not allocate in (V1.13) capacity 
resources to empty queues (as long as they remain empty). 
If one "deletes" from problem (VI.13) those commodities for 
which qs (0) = 0 (the expected time needed to empty those queues 
is not well defined), then its interpretation as a "minimal ex- 
pected time to empty" is strictly valid, and the network could 
Operate for some finite time with links saturated while also 
obeying the capacity assignments. 


Using the transformation 


E A 4 ef, A ES 


OF 0 


and introducing slack variables we obtain the LP formulation 


OR PIEI) in standard form. 





METP (1): min ty 


| 
5 
il 


k 
q; (0), MC Ko) = No 


=C- CM H DAS = 0 (VEAN 


C 
|v 
© 
+ 
H- 
LI. 
Ty 
E? 


The reader will recognize problem (VI.14) as a statement of 
the rst Minimal (deterministic) Time Problem (cf. III.lb) 
with constant rate data inputs {ay} to the network. Although 
both formulations are mathematically identical, here we inter- 
pret ta ore netexpected value of a delivery time. 

It is easy to see (cf. VI.14) that as ty reduces to its 
minimal value E a subset of links, say Ly, 1s bound to become 
erraleal, For this set of links the second constraint in 
(VI.13) holds with equality or, equivalently, the respective 
Slack variables in (VI.14) must be identically zero. 

It can be expected that many feasible capacity assignments 
will solve problem (VI.14). Because of this circumstance, we 
may ask for the "best" among many solutions. One way to ap- 
proach this question will be to apply the same min-max criterion 
A ETE (tte those Links which are not critical, i.e. 
comia So such that [1,3] Z L,. We therefore seek next to 
minimize expected delivery time to while retaining = tor 


all commodities that were assigned capacity in the set a 
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we denote the set of these commodities by N The Second 


1* 
Minimal Expected Time Problem can then be stated as 


Miele (2 Js min t 


a; (0) 0 | 
+k _ TK = tiy ae ee E N, 
ia at + 


k 
a ( 15) 
A 2 6 AK). EEN EN VER 
Be LE: -aj Es m 
7 (#i) Zi 
k ora 
y Be = Moca, Alte L 
k(X%i) i] 27 J 1 
mm £0) FO en eee 
A iy == aby O “1 
Es ES A EN 
7! 5 Z ’ J = 0’ ’ = 0 
Again, using the transformation 
z a n Dale 


and rearranging results in an LP formulation 
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METP (2): min to 


SS E. 
k 
q3 10) k k k 
meters t 2.) + i} un ee e E AO 
3 (41) 3 (41) 
k c _k k 
red + IE RE De = i ) = N- 
== ae j(i Ji A Sen 
+ Yun = 0, *[i,jl e L (vd) 
Be) en 37 f Á il ` 
SENC + ) u“ es = ¿Ola ce L.=1 
y 33 ij > Gee = OA 
Be > 0, *[i,j] « L 
27 ij’ Si; = , , = o 
x (i,k) <€ No 
0 
for given E; - 


may also happen that the solution to METP (2) is not 
unique. The procedure could then be repeated by identifving 
adenrsemalseritieal links [i,j] and defining L, to be a set 
of links including the new critical links as well as those in 
L, . Similarly we define N, to be the set of all commodities 
that have capacity assignments in the set Lo. Goneinuang to 
iterate in this way until we have exhausted all commodities, 
we ultimately generate a capacity assignment for which (VI.12) 
tS satisfied and thas is optimal. 

For completeness, we present the linear programming formu- 


lation of the m-th Minimal Expected Time Problem. 
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METP (m): man E 


Sol. 
qs (0) D: l 
an - u = 0, (i,k) e N -N 
0 , GS i , m , 
5 to A A AO pa 
DI-L ‚m-1 
K k k : 
-t_a. + Y u, = u. = qa: (0), »(i,k) € N¿=N 
m j (ži) q) 5 (Fi) Ja I Oia 
mec: - ) a == ei, ee as 1 (VI. 18) 
I a 2 
Bee + ) ue eG = Oye (aan) Ly-L 
J k(%i) 1) £) Ma 
te ur. uber 6. 
m’ ij = , 1] = 0’ 
* (i,k) <€ No 
: 0,0 0 
BOE Given Error ro 
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The SLO methodology establishes a hierarchical order among 
the critical sets, as indicated in Fig. VI.4. The arrows 
there indicate, for a given set of links, which are the 
commodities that may (must--for a horizontal arrow) use it. 
The question of identification of critical sets (or their 
partial composition) was considered in detail in Chapter V.5 


and we need not repeat it here. 
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Dicey isos Hierarchical Structure of Critical Sets 


The following example of a stochastic delivery problem 


illustrates some of the concepts that we have introduced in 


this section. 






capacities are indicated 
on links 


Mer ee Stochastic Delivery Problem 


ore 





Following the argumentation we used in discussion of the net- 


Work shown 1n Fig. VL.5 (cf. Chapter V.F) it can be shown that 


en Sr rst critical set N} is composed of commodity (3,4) only. 
In this case the assignment of capacity to commodity (3,4) 


Beomeneque and shown in Fig. VI.5a. 





ac oa. Optimal Capacity Assignment for Commodity (3,4) 


The minimal expected time to empty the queue q5 is given by 


Commodities (1,4) and (2,4) belong to N27Ny and the correspond- 
ing Capacity assignment is shown in Fig. VI.5b. The second 


minimal expected time is computed (using the values in Fig. 


Vie Sb) tO be 


Al 
O A 37° 
E 1 


CJ 
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EA SD. Optimal Capacity Assignment for Commodities 
in N_-N 


2 1 

we conclude this section with an observation regarding the 
optimality criterion in Definition VI.l. Since the mathemati- 
cal formulation of MTP(1) and METP (1) are identical we can 
interpret the stochastic delivery problem as a deterministic 
delivery problem with constant rate inputs. We could ask then 
for an optimal capacity assignment (that may change with time) 
which consists of two consistent parts, one that accommodates 
the constant input rates in such a way that the other enables 
optimal delivery of the backlogged data (in the optimal de- 
livery function sense). Since the modifications necessary 
to make the time (MTP(m)) and the rate (MRP(m)) problems handle 
constant inputs are trivial, this approach would result in an 
optimal capacity assignment schedule (utilizing a sequence of 


corresponding METP's and MERP's). 
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Although our previous results indicate that this sequen- 
tial optimization approach of interleaved time/rate problems 
is more powerful than that of time problems only, we do not 
use it because of conceptual difficulty that arises between 
firm scheduling of events in time and the potential necessity 
to recompute a new capacity schedule, For example, it is possi- 
ble that in an optimal capacity assignment schedule, a queue 
of some commodity (i,k) € No will be assigned capacity only 
following a point in time which is beyond the recomputation 
instant. This can lead, at least in theory, to tremendous 
delays in the delivery of that queue. For this reason we find 
1t conceptually more satisfactory to consider constant capacity 


assignment (as derived in this section) in our model. 


Er DISCUSSION 
in the last section we have shown that an optimal capacity 
assignment can be obtained as a result of solving a sequence 
of linear programming programs (METP's). We also have indi- 
cated that the mathematical formulation of a stochastic delivery 
problem is identical to that of a corresponding deterministic 
problem with constant rate data inputs. Due to this similarity 
we can apply the theory developed in previous chapters to derive 
and understand the structural properties of an optimal solution. 
The derivation of routing variables has been studied mainly 
in two extreme situations. In the first case only the steady- 
state arrival rates of messages are taken into consideration. 


The performance measure objective usually is to minimize expected 
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delay experienced by a message traversing the network or, as 
M RO minimize a sequence of saturation ratios. In any 
e eeehe underlying network model is that of Kleinrock [15], 
as described in Section A. In the other case, only the con- 
gestion state of a network is utilized to compute a routing 
flow schedule that attains optimal delivery function. This 
approach is studied in detail in this thesis. The dynamic 
delivery problem introduced in this chapter, provides a theo- 
retical model which will combine both types of information, 
i.e. the expected arrival rates of messages and the existing 
congestion in a network. We combine the objective of empty- 
ing the set of initial queues with the probabilistic informa- 
tion about future expected arrival rates of messages. We 
believe that the resulting optimal capacity assignment pro- 
vides a fairly accurate analytic model for a desired dynamic 
(short term) routing strategy (routing variables). How to 
use such routing variables (whatever their origin) in the 
implementation of an actual network control system is a com- 
plex matter and only initial" results are available. We 
will not consider this issue here. 

We have indicated before that it is unrealistic to expect 
that the new dynamic model can actually be used to compute 
routing variables in real network environments. We do suggest, 


however, that it may prove useful in simulation studies as 


ola 
'For relation to other results, see Chapter I. 


Mee cample 121], 122] anawi23]. 
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a reference against which actual routing strategies can be 
compared. 

DADS connection, 1t is worth noting that it is also 
possible to solve the dynamic routing problem which results 
when all of the as (not Just those for which qs (0) is non-zero) 
are involved in allocation of the ee In practice this 


J 
requires interpretation of the set N, in (VI.14) and in subse- 


0 
quent time problems as the set of all commodities. This does 
not change anything in the mathematical procedure used to solve 
those problems. The conceptual difference, however, arises 
due to the assignment of capacity resources to commodities 
for which a: (0) = 0. Since the message arrivals to the net- 
work are stochastic in nature, it now may happen that for a 
certain period of time the assigned capacity E (for the com- 
modities in question) will not be fully utilized. This model 
variant, although in theory inferior to the one discussed 
earlier, is probably an acceptable compromise between the static 
and the dynamic network models. Its obvious advantage results 
from the fact that it is not necessary to recompute the 
capacity assignment with every new message arrival to the net- 
work but rather as frequently as the actual computational 
facilities allow it. Finally, we mention that the new stochas- 
tic delivery model seems to be free from the "independence 
assumption." Moreover, observe that 1f all the = are set 
equal to one, the resulting (static) solution optimizes the 
objective function which (sequentially) seeks to minimize the 


maximum expected length. of the queueing system busy periods. 
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Vil. APPLICATION TO NETWORKS WITH TRAVERSAL DELAYS 


A. INTRODUCTION 
Up to this point we have considered multicommodity delivery 
problems for which the delay in delivery was caused by finite 
Srplemey sor network links. A natural extension of this model 
is to consider networks with traversal delays. Association 
of traversal delay with each one of the links complicates the 
mathematics of the optimal delivery problem, but provides a 
considerably improved model for transportation applications. 
The classical transportation problem (see, for example, 
[7]) refers to the shipment of assets’ from a set of sources 
to a set of destinations, to satisfy given demand at minimal 
Sect.) Imporcant class of extensions of this problem recog- 
nizes the existence of queueing and traversal delays, and conse- 
quently looks into the question of minimal time demand satis- 
fiability. We prefer to view this issue in a more general 
framework of Minimal Time Redistribution Problem (MTRP); given 
initial and desired distributions of assets, in some geographi- 
cal locations, the objective is to redistribute the assets 
deeerdimeiy in minimal time. Thus from our point of view 
there is no inherent distinction anymore between "source" and 


destination" nodes. 


"We use "asset" instead of the more common term "commodity" 
ocio aa sat from Our definition of commodity in Chapter 
II. There commodity was identified with destination node, 
where here one type of asset may be demanded in many locations, 
and a location may have demand for many types of assets. 
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meeps and important application of this class of 
Peeoremo is CO military logistics planning. In particular, 
assume that an outbreak of hostilities in a number of locations 
Bequires redistribution of assets (troops, tanks, supplies, 
etc.) in minimal time. The same model may be used for supply 
of aid to disaster-struck areas or for transportation of 
perishable supplies, and many others. 

The usefulness and applicability of minimal time problems 
has attracted great attention (see [24] for survey and exhaus- 
tive list of references). Most of the research, however, has 
been done in the area of bi-partite transportation networks. 
Hammer [251° provided an algorithm to solve a single asset, 
uncapacitated minimal time problem. Tapiero and Soliman [27] 
have treated the multi-asset version of the capacitated minimal 
time problem as an optimal control problem, using a maximum 
principle and a continuous state-space framework. Their paper 
does not contain proof of their algorithm. Bookbinder and 
Sethi [24] use basically the same approach but elaborate more 
on the mathematical programming aspect of their algorithm, 
for which only convergence to a local minimum is assured. In 
both cases it is unclear whether or not the algorithms are 
computationally manageable for problems of practical size. 

Ineiz4! an important observation is made, namely that at 
least for bi-partite transportation networks, capacity linking 


constraints (to be discussed in the next section) cause most 


'For similar results see [26]. 
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Ie coOmplexity in MPRP. The authors predict there that it 
may not be easy or even possible to include a capacity linking 
constraint and still provide a linear programming formulation 
of the problem. In this chapter we study this question and 
come to a conclusion that Sequential Linear Optimization 
methodology can be used to solve MTRP. Also, we analyze the 
Special nature of capacity linking constraints and their influ- 
ence on solution procedure complexity. Then we consider 
possible extension (by discrete time mođelling) of the results 
Porgeneral networks. Application to military decision problems 
1s provided in the formulation of the Maximally Delayed Decision 


Problem (MDDP). 


B. TRANSPORTATION NETWORK MODEL 
IE cios cal Representation 

A useful way to view the redistribution of assets over 
a set of locations is in terms of a network model composed of 
nodes and links. The links represent unidirectional means of 
asset transportation and the nodes represent physical locations. 
A typical transportation model is shown in Fig. VII.l. 

With each link we associate a traversal time, i.e. 
the time required by the corresponding transportation mode 
to traverse the distance between the locations represented by 
the head and tail nodes of that link, respectively. 

In general, we allow for a variety of capacity con- 


Sans eater most common of which concern loading, unloading 
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Base vıı.ı. Transportation Network 


and link capacities.' The loading and unloading constraints 
(sometimes referred to as "linking constraints") provide an 
upper bound on the volume of assets that can be loaded or un- 
loaded (in general, with respect to particular transport mode) 
Te en location per unit time. The link capacity constraint 
represents the upper bound on the volume of transportation mode 
and thus on the total amount of assets per unit time that can 

be sent over that link. 

With each node we associate, at every interval of time, 
the amount of assets stored at the corresponding location. The 
collection of these descriptors for all the nodes in the net- 
work and for assets in transit constitutes the state of the 


system. 


iene also may consider an upper bound on the amount of 
assets allowed at any given time at a particular node. 
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Since none of our results is dependent on the number 
of asset types nor on the number of different transportation 
modes we will, for the sake of simplicity, limit the notation 
to represent a single type of asset and a single mode of trans- 
portation. Generalization to more than one type of each is 
straightforward. Furthermore, since most of the notation and 
the meaning of network parameters are the same as in Chapter 


II, we will keep our discussion brief whenever possible. 


Consider a transportation network G(V,Lo) +» where 
Ve= T1172,...,n) is a set of n nodes and Ly =) is a 
set of links. We also define: 
q; (t) = amount of asset stored at node i at time t, 
Vera: 
f..(t) = rate of asset flow leaving node i on link 
= [i,j] at time t, w[i,J] < L, 
NO ae o lia aa Or Tof the transpor- 
1J tation mode represented by that linx), 
»[i,j] E Lo 
a, = loading capacity at node 1, *1 e V 
>: UN adiias Capacity at node jJ; ~] e V 
t.. = traversal delay from node i to node 3 along 
*J Laie er, ze 


We reserve the use of respective capital letters for set and 
vector notation, interchangeably. For example, the quantities 


of assets stored in network nodes at time t are given by 


MO a (E) a le) ra, (€) - 
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2. Dynamic System Equations and Constraints 


The quantities just defined must satisfy three basic 
constraints: non-negativity, conservation and capacity. The 


non-negativity constraint states that 


ee O “[1,3] € Lọ; *t, 
(Vile) 
qj it) = 20 Tie) E 
The conservation constraint may be written as 
> ) 
SER q. T1t.) - j ( Eada 
ne 5 ty Hal 
t, (VIL: 2 
- IO da eV 
a De 
and 
t > ti: 


Somoteommemyil.2) accounts for the fact that flows arriving 
at node i, i« V over link [j,i] at time t, left node jJj at time 


R where "a is the traversal delay associated with link 


ipl). ewermally, the Capacity constraints are 


o las Loose 


ij ga 


eee eee cog Ve V (VIDEOS) 
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OS IM ~] 


= 


Deramıceıion VIL.L 
A set of flows F(t) is a feasible flow schedule if it 
Sermrori1es (VIL.L)-(VILI.3) for all t. 

Let Qo and en denote the initial and the desired 
(terminal) distributions of assets (network states), respec- 
tively. We say that en is reachable from Qo 1f there exists 
a feasible flow schedule F(t), E, SS Es such that F(t): 
Qo (to?) > Q,(t,) and 0 < t, < t] < œ. Consequently, for any 


Such. pair (QoQ? there is some minimal value of E; - We 


define the minimal redistribution time ES as 


ct 


I 


HE EE O AO) OE) (VITA) 
ne e 0 ie =: 


Derinition VII.2 
We say that a feasible flow schedule F(t): Qa (0) > Q, (t,) 1s 


wee emal solution to the minimal time redistribution problem 
0 

LF ti = Er: 

C. BI-PARTITE NETWORKS 


de, Problem Statement 


In this section we study the minimal time redistribution 


problem on bi-partite networks. A bi-partite network 1s one 


"we of course exclude the trivial case where Qa = Q.. Also, 
we assume without loss of generality that the initial time Eo 
is identicallv zero. 
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whose node set can be partitioned into two subsets S and D, 


so that each link has its head node in S and its tail node in 


D. A typical network of this nature is shown in Fig. VII.2. 
ə ; 
S D 
Fig. VII.2. Bi-partite Network 
It is customary to associate the set S with "supply" 
nodes, and the set D with "demand" nodes. Various amounts of 


assets are stored initially at each of the I £ |S| supply nodes 
and there is a specified requirement for assets at each one 


of the J £ ¡D| demand nodes. Following our notation we write: 


A 
Qo io (S¿1S7r--«15710,0,...10) 
(VETS) 
A 
Q u ( ır 1 1 101 Ayres era) 
Teo eea on (~-) indicates "don't care" situation. -It is not 


important how many assets remain in the supply set S (as long 


as these are non-negative quantities). 
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we can restate now the Minimal Time Redistribution 


Problem as follows: 


For any given pair of system states (07-21) and net- 


work parameters T = { ae = oT A = fa, }, B= be, 


J 
Miami ly, 2,...,l and j = 1,2,...,J, find a minimal time 


La 
+) 


flow schedule. 
2. Structure of the Minimal Time Flow Schedule 
In this subsection we analyze the structural properties 
of the minimal time flow schedule. We derive a result which 
parallels that of Thm. II.1l and enables us later on to formu- 
late the minimal time redistribution problem in LP form. We 


Serctewhen the trivial but important observation that 


Lemma VII.Ll 
LES A AS ti be a feasible flow schedule such that 


E(t): Qg (0) = Qı (t1). Then we may always take 


WCA = 


fx 
17) | 
SEE 


PEOOF: 


Eger (VIOD 


Consider a bi-partite network (like that shown in Fig. VII.2) 
and assume that there is some feasible flow solution F(t), 


eee SUC that Be je: 29,10) > Q (tı). It is quite oọobviġus 


I 


that this flow schedule cannot use any link [i,j] for which 


west since all the flows on this lirk will arrive at their 


13 ke 
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destination j later than SE and thus be of no use in the con- 
text of the minimal time problem. Also, there is no poamtsin 
mung clows over a useful link [i,j] ae < ty) beyond time 
m ij because they will not reach their destination in time. 


ra 


== 


Combining (VII.2), (VII.5) and (VII.6) we can express 


the desired system state with respect to the demand nodes as 


con | eee (tetera un (VII.7) 
1) 
O 

Similar iy, the non-negativity constraint (VII.l) with respect 


to the supply nodes can be written as 
E a ea e Ss (UI) 
J 


The minimal time redistribution problem reduces to 
finding a flow schedule which satisfies (VII.7) and (VIT.8), 
subject to capacity constraints, in minimal time. 

We now show that the search for a minimal time flow 
schedule may be confined similarly to the optimal delivery 
problem, to the class of piecewise constant flow schedules. 
Furthermore, it is possible even to narrow this class to flow 
schedules which we call linking flow schedules. In order to 
present this subclass we need to introduce some additional 
notation. 

For a given bi-partite transportation network we define 


mortech Supply node i, ie S a vector T; E (1, (1),7,(2),...,7,(n,)) 
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of all distinct traversal times associated with outgoing links 
of that node. The components of T are assumedto be ordered 
so that T; (K) < T,(r) Is nes and n; is the number of 
vector components. By n.. we denote the ordinal position of 


1) 


Be mene or, such that T.[(n..) = 1... For the network 
Y E] T) 


in Fig. VII.3 we have 


0 Traversal delays are 
indicated on links. 


© ad 
5 D 


Fig. Vil.3. Bi-partite Transportation Network 


ere f= (2,10), T, = (8) 

ny = on ny = a n, = L 

A 2 el ya a a 
n = 1 
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we now use this example to demonstrate the structure 


of a linking flow schedule. 


To f 4) E, 1 (2) 
I ' 
` pa 
£15 (t) 
Fe) 
27 
En | 
| DM: a 
2 SE) 4 
u 22 ae 
23 ' 
EET) 
Ze £,.(1) 
33 
a; 1 £33(t) 
N —— 
E, -10 t,-8 t 77 E,73 t73 ty 
Fig. VII.4. Linking Flow Schedule (ty > max Ty 4) 
L J 
0 
The length of a horizontal line in Fig. VII.4 expresses the 
time duration of effective flows on that link (Lemma VII.1). 
The labels above a line identify the constant segments of a 
flow. The crosses (or corner points) indicate the time instances 


at which the flow may change its value. We do not consider 


Flow termination points as corner points in this context. 


"We consider loading and link capacities. Unloading 
capacity constraints are studied later on in this chapter. 
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The flow scheäule in Fig. VII.4 was obtained by using 


the following construction rule. 


Construction Rule: 
A a on a link [a,b] has a corner point at time t iff 
link [a,b] is linked by a loading capacity constraint to some 


and T Sr 


Aena [c,d] such that t = t.-tT od 


i ea ab’ 


Now we are in position to formulate a general descrip- 
tion of a linking flow schedule for the case when ty > max Tis! 
a 
i.e. all the links are usable. 
behumLe tone VII. 3 


A feasible linking flow schedule F(t), 0 < t < t, is described 


by 
E = fij, ty-T,4 (m1) AS E (VI 
ur E E o TON ESSE Se 
A _ we 
where Ta = a is å vector all distinct 
names Mdeleays tot links linked to [1,3] by a loading con- 
| Toe | A 
straine, lia) . Ly: Also, we define an = t]. 


Tof each of the asset-types, in multiasset case. 
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For notational convenience, we define ©, (m) to denote the 


duration of the m-th flow segment f; (m), 


ID 


9, (m) a - Tij w: m = I (VILE TOJ 
Since our interest lies, as indicated before, within 
the class of linking flow schedules, it is useful to rewrite 
the constraints (VII.l)-(VII.3) (we also use (VII.7) and 


(VII.8)) as they apply to linking flow schedules. 


(i)--non-negativity 


g A e Lo’ ne ensure 
(VET TI 
Be 
1) | 
Sr ) ) £55 (mo, (m), wies 
J m=l 
(ii)--conservation 
n.. 
a (VII: 12) 
d = O A e ED Ver: 
J o 1) > 7 
(iii)--capacity 
Ze) = ES A € Lo: m= 41,2, mas 
(VII L33 


Finally we state and prove the main result of this 


section. 
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Theorem VII.l 
Let Qo and Q. be any two states of a given bi-partite trans- 
portation network such that 9, (44) is reachable from 9,10) by 


some flow schedule F(t), O Secs E, - Then there exists a 





feasible linking flow schedule F(t), Oe tos ty such that 
E(t): Qg (0) ES Q. (t). 
Proof: 
Define (for useful links) 
E A 1 aa ¡340 
EY ~ 5. m) ENE ae (VIT. 119 
a t_-T..(m-1) 1) 


lien) € Lo y m = Loto 
If we start now in state Q, (0) and apply the new flow schedule 
P(t) as defined in (VII.14), the system will be transfered by 


time Es into a new state, say Q. (t1) rA VIT 2); suchithat 


11 


a 
a IAS E “j e D (Vid sess 


Il oa.: 


oa) 
A O 


Spa par. 14) into (Vil.15) results in 


>) f,.(t)dt, j < D (VIL. ep 
1 e 


) = d (ti), wj e D (VIL 
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Equation (VII.17) is essentially the desired result 
but we still must show that F(t), as defined in (VII.14), 
ees Constraints (VET.11)=<(VII.13). The first part of 
(VII.11) is trivial due to the requirement fij t) > 0, 
way) € Lo’ vt. The second part of this constraint is shown 
Mmemoe Satistied by substituting (VII.14) into it and comparing 
PAresresult with (VII.8). 


To show that Zy a los Los A A 


u i 
merecoemelude from (VII.14) that 
fij (m) < max tf, (80), (VIT. TD 
[t] TT; : ana) Se (m) ] 
ea eee ar i xa. 
1] 


But since by assumption the right hand side of (VII.18) must 
zer scherlınk Capacity constraint, so must DEE 
Desshew that the loading Constraint holds for the new 


flow schedule, we write (by assumption) 


) fis(t) < ay, vie S, ve (VII.19) 


Mee beth sides Of (VII.19) over the interval E 


£ . 
ty Denen fornfany mn such that 1<m< nij results in 


(JE; 5 (t)) at a (VIT-20p 





Pee vit.20) can be written as 


pa 
Aq 
H- 
t. 
> 


(m) < a; (VII 2238 


QU. 


1] 


use = ds aa y Mo 
1] 


which completes the proof. 


An obvious consequence of Thm. VII.1, with respect to minimal 


time flow schedules, is given in the following corollary. 


sereolarv Vil. 

If there exists a feasible minimal time flow schedule then 
there also exists a feasible minimal flow schedule of the 
sino type. 

U 
It is worth noting the difference between the optimal solution 
to MTP (1) for a network without traversal delays, and the 
minimal time flow schedule here. In the first case it was 
sufficient to consider a constant flow schedule, whereas here 
we need a piecewise constant solution, which serves as a first 
indication of higher complexity of delivery problems on net- 
works with traversal delays. 
3. Solution Algorithm 
Armed with the results of the last subsection we may 


state the Minimal Time Redistribution Problem as fcllows. 
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ME REE bl 


I 
Sre. 
DE 
1] 
f mom = d., *¥ 
O Om) AO 
nee 
. Ñ 
2 A O a) SS MOL © S (VARIOS 22) 
en a a 1 
Es YE e Cis a[i jE Lo my as 
) E 3 110) < ass AS ML idos an 
J 
th E ¿Un > 0, > la € Los nee, 1 
g 


We have assumed that = > maxtt;.t, ena nas are 


usable. 


Spec Ion Of the first two constraints of (VIL.22) 


{T 


SONS that in both the coefficient of zen A € Ly 


hħas the form (VII.10): 


ô. (1) = t, = 7,0), “ies 


tie tee OF this coefficient 1S unknown since it is a func- 


tion of the variable tje Define 


a 5 ~ c x 7 
ujj H = £,¿0)8, (1), e Lo (VITAE 
Using (VII.23) and rearranging (VII.22) results in 
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MTRP: min t 


I 
Such 
De 
iD 
wel) + A IS a, , “ie D 
So] i m=2 + o J ; 
13 
A IL E O << s., vies (VII.24) 
Fe s A ij T = 
i j a 
-t è + -T, l s : € 
ne „eo = ‚lle,., *[1,)] Lo 
Ly (m) = ij’ vi, = Lo’ 
m = 2, E 
17 
-t,4, + jp E le “wies 
cag HU = aj; “ie S, m= T 


En am f..,(m) 


| v 


O Ll E Lo’ m = 2 ON 


which is a linear programming problem. 

Mas noularbe noted that the formulation in (VII.24) 
is valid only when t)-t; (1) > 0, *1 e S which makes the inverse 
transformation of (VII.23) meaningful. But this is taken care 


of by the assumption that 20 


1 
1 > max LT ig Let 


j 
Do 


T 2 (Weert) be the vector of all distinct traversal 


delays of the bi-partite network. The components of T are 
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geenmedwteesbe In descending order of their size, and n is 


their number. Thus for example the statement: = armaz eee 
L 


is equivalent to 0 


Ele (VIT. 29) 


Suppose now, as may be the case, that a Er, aa 


1 
particular that 
0 

T(k+1) < Es r (K) (VII. 26) 
for some k, 1 S, It should be clear that the solution 
to MTRP will yield (due to the nonnegativity of ut), 
yl e Ly) 

0 
ti =, (Vil 22% 


HoWaeton (Vil.27) suggests itself as an indicator of the fact 


Pat = er. We may think Of solving MTRP again, but this 


time we "remove" all the links for which ze = Til) It may 
happen now that ty = 1(2), in which case we repeat the proce- 
dure with respect to links for which Ti =a 2 oer LE 

0 


T(2) < t} < t(1) then we have found an optimal solution. 


Pofinition VII .4 
fae macnn > 1, Minimal Time Redistribution Problem (MTRP(m) ) 
is equivalent to MTPR for which all the links [i,j] = Ly, such 


chat Ty. > t(m) are considered to be not usable. 
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Actually, we are free to search for = in any order, 


and in particular we may implement a "binary search" technique 


[28]. The resulting algorithm will have the following form: 
Algorithm: 


=- fal 


LOOD 
Solve MTRP (m) (VII 225) 


er EN = T(m+1) then m=- m+ E 


repeat loop. 
: 0 m 
sibe ty > t(m) then m + E repeat loop. 
else stop, 


where n is the number of distinct traversal delays in the 
network and [x] denotes the ceiling (i.e. the smallest 
ese 1 Of the number X. 

L 
Since the complexity of a "binary search" is logarithmic with 
respect to the number of elements being searched we conclude 


that 


Lemma VII.l 

The MTRP can be solved in 0 (1g, n) number of steps, each 
being a solution of an LP, where n is the number of distinct 
traversal delays in the network. 


In the next subsection we give a simple example of MTRP. 
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4. Example 


a, =4 
Os 
Sz Aa) er 
S, =50 
10 
5 
3 
a Ens E Lo 
GH & 
2 
s, =30 _ 
2 a, =4 dz =40 


Fig. VII.5. Minimal Time Redistribution Problem on 
Bi-partite Network 


The MTRP shown in Fig. VII.5 was solved using the methodology 

of the last subsection. The resulting minimal time flow schedule 
is shown in Fig. VII.6a. We find it useful to present the 

same flow schedule from a demand node point of view, i.e. 

with respect to the arrival time. This is done in Fig. VII.6b. 
We conclude the example with the observation that the existence 
of unsaturated links and/or loading constraints (cf. F511), 

f 


De fna (2) ) in the optimal solution suggests the optimal 


Do 
delivery function as the final goal of optimization. We leave 
this point open for a future study. We will touch again upon 


this issue when we consider discrete time networks. 
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£. (t) ma £, 4) =2 
17 
ze iE (2) =2 
a, £ (t) 12 
f} 3 (t) f3 (1) =2 
P fn, (1) =0.67 
4 | £, (t) £n, (1) 087 
E = = = 
£43 (2) £5 (1) 2 £3 (2) 4 
0 2 0 
a iS 1475 t,=17.5 E 


oi ca. Minimal Tame Flow Schedule (flow departure) 





fo, (t) £ (1)=0.67 
fn, (t) £5 (1) =0.27 
f3 (t) £, 301) =2 

fn (t) fn (1) =2 £,,(2) =1.14 





0 3 5 SEO EERS E, =17.5 8 


Fig. VIl.6b. Minimal Time Flow Schedule (flow arrival) 
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EA ETORRI Ting Unloading Constraints 


Up to this point we have considered the MTRP (on bi- 
partite networks) with link and loading capacity constraints 
only. In this subsection we introduce the unloading constraints. 
We will use the network in Fig. VII.3 for illustration purposes. 

The unloading constraint introduces linking of the 
flows with respect to their arrival time at demand node set D. 
The loading constraints had the same effect with respect to 
eirescdeparture time from supply node set S. Fig. VII.7 shows 
the unloading linkage of flows. (Note that from a demand node 
point of view, a flow on link [i,j] commences at time Cas and 


terminates at time Ey +) 


£ \ A 

N | 114) £,10) EZ 

l A === O 

f, (t) £,] (1) 
fi> (£) £, es 
fo, (t) £,5(1) £,,(2) 
ao a 

uy a eg 
f(t) £33 (1) 23302) 


A a Linking by Unloading Constraints 
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when both loading and unloading constraints are in 
effect, problem formulation is complicated by the requirement 
that the two grids of corner points, namely that in Fig. VII.4 
En Aa in Fig. Vil.7 must be consistent; if any flow variable 
Changes value at a particular point of time (corner point), 
Pruentatt other flows linked to it (now by both constraints) 
must also be allowed to change. Consider for example the 


Paeieue lal network in Fig. VII.8. 


Fig VII.8. Network with Linking Constraints 


The sequence of plots below illustrates the procedure of find- 
ing a consistent flow schedule. In each figure, the crosses 
identify already established corners, while the circles denote 
new ones. 

After some thought the reader can convince himself that 
if the traversal times (and/or the solution time es) are 
incommensurate, the number of corner points would be infinite, 
in which case the solution is not piecewise constant. This 
difficulty may be overcome by approximating the continuous time 


axis by a finite sequence of discrete time values with some 
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£, (t) 
0 E, To E, Ea E 
PE — o 
£, (t) 
N seem 
0 Ta T, E 719FT7 ty 
f (t) 
£, (€) 
0 Tat] t1 To t 


Fig. VII.9. Linking Flow Schedule 


resolution At, At > 0. We note that the same reasoning applies 
to general networks with traversal delays anc hence defer further 
discussion of discrete time networks to the next section, where 


we no longer restrict consideration to the bi-partite network. 


DS DISCRETE TIME APPROXIMATION 
In this section we briefly discuss a discrete time approxi- 


mation of a minimal time redistribution problem. We have shown 
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already that this approximation arises naturally in the con- 
text of bi-partite networks with both loading and unloading 
constraints. The same reasoning applies to general networks. 
Hence we must assume that the corner points of a minimal time 
flow schedule may occur at any integer multiple of some basic 
ieee OF time At, At > 0. 

It 1s worth mentioning that for general networks the number 
of corner points in a minimal time flow schedule can be very 
large, even if only link capacity constraints are present. 
This is so since a flow over any link in the network has to 
reflect all the possible time patterns of flow arrivals (in- 
tended for that link) to the head node of that link. Although 
this phenomenon is unrelated to the notion of commensurability 
its effect for large networks may be in practice the same. 

We start by breaking a given interval of time [0,T], where 
Bene ee < hee intowk = i eepabes.weA feasible flow schedule 


IN SE SE As assumed to be constant throughout each 


1 
segment of duration At. Also, the traversal delays are approxi- 
mated to the next nearest integer multiple of At. Thus from 
here on we will interpret TE as the integer specifying the 
number of basic time units At which make up the traversal time 
from node i to node j. We use the following notation for the 


eeann desired states (i.e. distribution of assets) of a 


network: 


p> 


Q = (S¿1S71+».15,) 


A 
= (dyrdyı...,d,) 
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enesalso the convention 


2 


Ke 


A . 
eg = the amount of assets shipped out of node i, 
SVemmaw linge, |) at ElMeeitistant rat, 
and arriving at node j at time instant 
E E as 
q; (r) = the amount of assets stored at node i 
enr roughout ENeNI=tn interval, er = 0,1... 
We also observe that a link [i,j] can carry useful flows 
.)-th inter- 


that were initiated up to and including the (K-T; 


val. Any flow launched beyond this time will not reach node 


¡Moby tıme K-At. 


We formulate now the minimal time redistribution problem 


as 
MTRP: Find a minimal value of a time segment index k, say ae 
for which the optimal value of the cost function in the 
Co lowing linear program is zero. 
LP: min a 
Si. Le 
a q; (k) d., *l e V 
e A 2020) s., *i< V 
z ay ae : 
E O] EA) E T pn 
= ny j(zi) J% : 
ines ape 
"For Tasse ronal simplicity, here and in the sequel, we use 
A ee tehougn f.. (r-t...) = 0, unless 0 < ET < K-T,, 
ren) e L 32 J z J 
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) air) Ra ren SOL y K 
MA e 
(Vis 
Terme) a Re ee er: 
wo gt eL Ì 
a ] E = 
Po) Euer *[1,3] e Lg, r= 0,1, y K 
Ur q; (r), Ejz (E) A OF “l < V, ALG j] E Los 
SOF E 
where Qg and en are given. 
0 


It is not difficult to see that k At approximates the 
optimal value of the minimal time =o within one basic time 


interval. More precisely 


0 


(k°-1)at < to < ae. (VII. 30) 


0 
i 
This is achieved of course at the expense of solving the LP in 
(VII.29) a number of times, for different values of the time 


segment index k, until the smallest value of k, say Ko 1s 


discovered for which min a = 0. Using a binary search tech- 
nique (similar to (VII.28)) results in needing 0 (lg, i) steps, 


where AT is the size of the potential range of values of oe 


The additional computational complexity that is introduced 
by this approximation scheme is counterbalanced,at least 
Barea lily our ability to implement thenotion of an optimal 


delivery function. We remind the reader that in the case of 
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exact solution (of a bi-partite network without unloading 
constraints) it was in general impossible to implement any 
requirements with respect to delivery time (the optimal delivery 
function and in particular the minimal rate problem fall within 
this category) while preserving a piecewise constant flow 
schedule. 

The discrete time version of a delivery function (compare 


weh (11.6)) is given by 


De ee, x = 0,2,..+,k° (VII.31) 


Ma 
where the summation is over all nodes for which the demand 
is defined (demand set). The delivery function in (VII.3l) 
was defined for flow schedules which do not allow intermediate 
storage or, in general, the amount of assets at any node i and 
any time t cannot exceed the demand A. at that noce.: In Ene 
redistribution problem this assumption is not true anymore 
and a generalized formulation of delivery function is necessary. 


We define 


a a (VII.32) 
ep 
where 
q; (r), q; (r) < d; 
w; (r) = 
diy otherwise 
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ae ormulatıon in (VII.32) accounts for the fact that a 
surplus delivery of assets (beyond the demand d,) is not con- 
sidered helpful and thus must be disregarded. The reader will 
Botıce that if q; (r) = diy *r, as the case is in the delivery 
meoolem Of Chapter IL, Definition (VII.32) reduces to (VII.31). 
Since the corner points of the minimal time schedule are 
fixed by the approximation method and the minimal time index 
solution k5 the optimal (generalized) delivery function is 


defined as follows. 


Derınitien VII.S 
A delivery function DY (r), = Glee is said to be 


optimal if it satisfies 


Dr) = max {D(r)|p°(k°),D°(k°-1),..-.,D° (r+1)}, ass) 
Bas) 
r = io 
where 
Boa) > ede 
: a 
leD 


(23) 


The sequential linear optimization procedure now consists 
of solving the minimal time problem (VII.29) followed by a 
sequence of x9-1 Macao lens lie formulation of the m-th, 
0 


m = k P | problem in that sequence is derived next. 


We start by observing that 


max D(r) = miní ) d, - D(r)) 
ieD 
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Semeauivalencly 


ieee ae = min Beir) | (VII. 34a) 


where 


| d.-q, (rx), Gave. q; (x) 


ab 
d; > w, (r) = | 


h, (r) 


O, otherwise 


We can state now the m-th optimization problem as follows: 


min ) h.(m) 
i 
le 2) 





Sit. 
0 E 
q; (K ) = dis “i < D 
A 1 EOS) ISLE Y 
i eO u 
Pr) - ga.(e-))4 ) E£..(r) men) = 0), vice V, 
> Š A eee SO Aa 
mos Men Me 
Bega (=) > d., vie D 
a = eee at 
ai 


See (VII.29). 
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0 
mD 1 


y fi Ea a avs meer Ver = 0,1 ont 
jie) 
0 
E (r-T ) MDEE RNV r = 0,1, ar VII.34b 
es eos ni | | 
Sita ales) < eee |) e E r= 071 = 
lj Dr PA t o i f 7 
NS A E AA lil € Eg, 


where Qg” SE nn and Den r = ee are given. 


m 
& 


Basis zimrortant to notice that the definition of h, (1) in 
(VII.34a) is satisfied in the linear problem formulation 


(VII.34b). The underlined constraint and the form of the cost 


á 


AA eeoa n anb) provide for this fact. If q; (r) > d, 


then the cost function will force the corresponding A, (x) 
to be equal to zero and when q; (r) < d; then h, (r) is exactly 
the difference between the two values. 

We use the following example to illustrate some of 


the ideas presented in this section. 
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Example: 


= 3 
= 
ST E Lo 
dy =d, =s, =S} =0 
=]3 





Fig. VII.l0. Discrete Time Redistribution Problem 


For the problem in Fig. VII.10 the traversal times (in units 
of At) are indicated on the corresponding links. The relevant 
loading and unloading constraints are shown per unit of time 
At. This problem was solved using the methodology presented 
earlier in this section. The resulting minimal time flow 


snc dule iS schematically illustrated in Fig. VII.1l. 
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Flow 
£ 
£3 (t) 
a,=2 
Ey (E 
£53 (E) 
a,=2 
iog 


Bion VII. ll. 


The same flow schedule is shown with respect 


meld.  VII.i2, 


Flow 


Ben E 


Shipment instance 
0 I 2 


Arrival instance 


I 2 3 
1 0 

il 

2 
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Minimal Time Flow Schedule 


4 





(loading) 


to arrival time 


0 
K 
6 7 8 
2 2 
2 2 
L 
> ~ 


Minimal Time Flow Schedule (unloading) 





Besoprrmalsgelivery function is shown in Fig. VII.13. 


Bie) 





Pacey tte Se Optimal Delivery Function” 


An observation with respect to Fig. VII.13 is appropriate here. 
In contradiction to the optimal delivery function studied in 
the main body of this thesis, here the objective function is 
no longer convex. This is still one more distinction between 
the delivery problem on networks without and with traversal 
delays. 

We have mentioned in the introduction to this chapter that 
the minimal time redistribution problem has important appli- 
cation in military logistics planning. The SLO methodology 


and the results obtained thus far allow us to extend this 
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Se 2 ment beyond the classical transportation preblem, into a 
domain of complex decision making processes. We study one such 


application in the next section. 


E. MAXIMALLY DELAYED DECISION PROBLEM (MDDP) 

In this section we are concerned with a particular aspect 
of the decision making process in a military environment. 
Assume initially that a possible strategy to follow has been 
determined and that the problem of interest is the redistribu- 
tıon of available assets to implement this strategy. Since 
the various possible distributions of assets over a given set 
of locations constitute a state space of our model, any strategy 
implies one or more corresponding state trajectories (from initial 
to desired distribution) in that space. If there is a set of 
possible strategies under consideration, the ability to analyze 
their corresponding state trajectories can help identify which 
of these strategies is most attractive. 

Consider a military environment characterized by a network 
G = (V, Lo? and an initial distribution of assets Qo = (Sjre..,8,), 
and suppose there is a set P = (Py Pores +s Py) of M strategies. 


Each strategy is described in turn by its corresponding terminal 


distribution of assets Qn! nN =e eee ee We assume, without 
loss of generality, that all the states On’ oa | cee, SM ame 
reachable from Qg’ 1.e. there exists rd: Qg (0) > Qan (E? Tör 


some K < ©, which is called the horizon time. The decision 
maker is faced with the problem of selecting the most desirable 


eaa ary context) strategy from the set P. An inherent 
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FOPerty Of most military situations is the necessity to make 
decisions subject to a high degree of uncertainty (as to the 
enemy's state, for example). It seems natural that the deci- 
sion maker should seek to delay the decision process as much 
as possible, while still ensuring that all the terminal dis- 
Bras utions E) are attainable within the prescribed time limit 
K. By delaying the decision instant, a number of advantages 
are achieved, for example: 

(a } minimize unnecessary, sometimes irreversible, 
commitments of assets. 
(11) maximize enemy's uncertainty as to which strategy 
is selected. 
(iii) gain time to acquire additional information which 
may influence the decision process. 
Aces Step in delaying the decision instant by, say, 
Ky units of time, we consider finding a common flow schedule 
Pix), UNS Ky ntc easter Ene system from its 
initial state Qo (0) to some intermediate state O(n) Ob- 
viously, the flow schedule Flik), the state 0 (k,) and the re- 


sulting decision delay of k, units make up an acceptable solution 


y 
to the delayed decision problem iff all the terminal states 

Q , m= 1,2,...,M are reachable from 07 (k,) within the remain- 
ing time K-Kk,. This leads us to the statement of the First 
Maximally Delayed Decision Problem. 
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MMDP(l): Find the largest time index K], Say Ke, for which 
the optimal value of the cost functiön in the 
following LP is zero. 


Be: min en 
Sec. 
I Er 1 
EAS) 0,0) Q (k,) 
(VETAS) 
- 3 — 

F A) Q (k,) a Or (K) - 4), m= IN 

where a, = (Uy 1 q ro. 107) is an n component vector and a, > 0. 


a 


To simplify the notation we have used vector formulation, 
where each of the constraints in (VII.34) represents a short- 
WwemeOneall the constraints in LP(VII.29) with the distinction 
That O” (k,) here (which corresponds to 91 = (dq dar...) 
there) is unknown. 

It is true that as Ky increases to its maximal value i, 
some subset of trajectories connecting QO" (ke) 


tive terminal states is bound to become critical, so that the 


to their respec- 


corresponding terminal states become unreachable from ot (x) 
for any k > ky We denote by Pr the subset of strategies that 
become critical at k7. The decision to select any one of these 
strategies has to be made prior to or at time n 

er Ps c P, then we may follow the same line of reasoning 
as before to formulate a second optimization problem, which will 


let us identify the second critical set of strategies and 


enci corresponding critical time. 
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MDDP (2): Find the largest time index k., say oe Torsushnich 
the optimal value of the cost” functioñ in the 
following LP is zero. 


P: in 
L mın 5 


F (k): @g(0) - otırd) 


0 
Fk): OKI) > 0Q,(K), eme P, 
Be Or B of (k3) (VII.35) 


2 
Baik)? Q (ko) > Q.(K) - a,, ¥m É Pr 


= 0 = 
for given Kis where a, = (Agr Agresa 


> and A, > 0. 


>) 


El 


Continuing to iterate in this way until we have exhausted all 
strategies, we will generate a finite number Mor Mo Á M-1 of 


pairs (Kup), m = 1,2,...,M, with the following properties 


0 0 
(i) 0 < Ky < Ks < < Mo < K. 
(11) P. n Pe = 9, i,j = 1,2,...,Mp+ A (VII.36) 
-ku 
In 
(211) UP a= P 
isn 


We see, as before, that the SLO methodology gives rise to a 


wrae hical structure (Fig. VII.14) in which the solution of 
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the individual optimization problems is carried out in a 
sequential manner, and that each such solution constrains the 


solution space of problems lower in the hierarchy. 






0 Qo 
yo Q (Ky) 
T 
| 
{ 
| 25 0 
0 ! Q (kJ) 
K 
2 
i 
f ! 
: 
i 
¢ | | ae 
i | 0 
E i ! un 
¿ 0 | 
| i i 
| | | 
! | 
| 
| f 
| ! 
| | 
| J | 
K = ° 4 
P P E 
L 2 Mo 


Fig. VII.14. Hierarchical Structure of the Strategy 
Sets w/r to Decision Instance 


We conclude this section with an example of MDDP. 
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Example: 


n 
t? 


Q 
H- 
U. 
| 


EL 


a 
En 


A 


Ku 
iI 
O 
il 
8 
t 
Pp 
m 
< 


Fig. VII. 15. General Transportation Network 


Eere NENNE WOK in Fig. VIL.15 thé initial distribution 


of assets is 
Qo = (2070704070700) 


and there are two possible strategies characterized by their 
respective terminal distributions Q1 and Q5 at time KAt, where 
K = 10. 


O = ORo 00 0,20), 2, = 1H 0. 0,0715,5.,0) 


The decision maker is faced with the problem of identifying 
a £low schedule which will enable him to delay as much as possible 


the final decision as to which of the strategies he is to follow. 
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The problem was solved using the methodology of this sec- 
tion. It was found that the latest instant of time by which 


a decision wnich strategy to follow has to be made is kyAt, 


where KY = 4, 





Beauv et boo Hierarchical. Structure of the Decision Process 


The first critical strategy set P, consists of strategy Po- 


1 
The fact that strategy Py 1s not critical (and the only one 
left) indicates that once it is decided to select py at time 
Ky = 4, it is possible to meet the desired distribution 91 
prior to time K = 10, say at time K* (see Fig. VII.16). The 
appropriate minimal time flow schedule can be found by solving 
a minimal time redistribution problem with Qo = ot (x?) and 


TE The complete flow schedule solution to the MDDP 


Inner TITS 1s given in Appendix H. 
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VIT r ONCLUSIONS 


A. SUMMARY OF IMPORTANT RESULTS 
In Chapter I we expressed the desire to obtain new insight 
into the problem of dynamic routing assignments. As an initial 
step toward achieving this goal a new performance measure for 
efficient delivery of backlogged data to their destinations is 
presented and called the optimal delivery function. The first 
important result is that the corresponding optimal flow schedule 
can be obtained by solving a hierarchical sequence of linear 
programming problems. Due to this fact the optimal delivery 
problem is computationally tractable even for moderately 
large networks. Furthermore, the well established results of 
linear programming are exploited to derive and understand the 
properties of an optimal flow schedule. Among the important 
results, at each hierachical level, the optimal flow schedule shows: 
(1) A critical set of commodities that share a common 
delivery time. 
Che A critical set of network links that must be saturated 
by those commodities throughout their delivery period. 
(iii) The minimal rate with which saturation can be achieved. 
Another important result relates, at each hierarchy level, 
the properties of an optimal flow schedule to the optimal cual 
variables associated with the corresponding linear program. 
Moret la Ene composition of the critical sets can be 


uniquely identified and the minimal rate of the saturation 
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flow can be calculated, whenever the optimal solution is 
stable. 

The optimal delivery problem is then studied in the case 
of single destination networks. Two important results are 
obtained. The first states that the optimal delivery function 
is globally optimal, i.e. no flow schedule can deliver more 
commodity to the destination at any time, than the optimal 
flow schedule. The second result is the algorithm for solution 
of the single destination delivery problem. It is composed 
of a sequence of linear programs whose size is independent of 
the total number of hierarchy levels, which makes it compu- 
Puemonally Crrtcient. In this context a method for identifying 
critical sets, when the optimal solution is not stable, is 
presented. 

An important result of this thesis is presented in the form 
of a new dynamic network analysis in which the optimal capacity 
assignment for routing purposes takes into account not only 
the backlogged messages but also the expected arrival rates 
of messages at the network. It is shown that the "minimal 
expected time to empty a queueing system" objective leads to 
a mathematical formulation which is identical to that of the 
deterministic delivery problem with constant flow inputs, and 
thus most of the previously derived results apply. It is also 
shown that a slight modification of the dynamic network analysis 
leads to a reasonable problem formulation which can be expected 


to provide a bound on the performance of more computationally 
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tractable routing procedures which might be implemented in 
a real network environment. 

In the last chapter a number of results concerning the 
delivery problem in capacitated networks with traversal delays 
are obtained. In particular we mention the solution algorithm 


for the bi-partite network case. 


B. AREAS OF FUTURE WORK 

Many additional areas of research appear to be ready for 
further investigation. 

One important point for future study concerns the efficiency 
ern spr0poOsed algorithm for solving the multicommodity 
delivery problem. In particular, decentralized computation capa- 
bility should be investigated for its efficient implementation 
in dynamic routing schemes. 

A further study of the delivery problem on general net- 
works with traversal delays is needed. Here the question of 
computational tractability seems to be crucial. Our results 
for the bi-partite and general networks should provide an 
go Opri ate starting point for this effort. 

Throughout the thesis we presented a number of conjectures. 
The questiors of the number of corner points in the optimal 
multicommodity delivery function and its possible global 
optimality,in the case of intermediate storage flow schedules, 


remain open. 
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APPENDIX 


INS ERCOE OF THEOREM II.1 

Thm. II.1l was stated in terms of flow schedules which satisfy 
constraints (II.1)-(I1I.3). We prove this theorem here for a 
wider class of flow schedules, namely those that allow inter- 
mediate data queuing. More precisely, we say that a multi- 
commodity flow schedule is feasible in a wide sense (w/s) if 


ME satisties the following conditions: 


ER 


(i) fij (t) Bee Ore Vl E Lo and *t. 
(ii) arte) > 0, w(i,k) e Ny and vt. (A.1) 
a A A eb, and +t 
k(#i) 7? y 
where 
= 
k k k 
re = LL) ESE) = 0). £3,(€))dt: 
u rn j(#i) 7° 


The basic difference from Chapter II is that here the delivery 


Rare or commodity k from node i at time t, 


k A k k i 
ao ) Seele (2932) 
> si) + j(ži) J7 


is not restricted any more to be a non-negative quantity. 
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Theorem II.l (Extended) 
Let Q (ta) and Q(t,) be any two states such that Q(t,) is reach- 
able from Q (ta) by some feasible (w/s) multicommodity flow 


schedule F(t), to SS t,- Then there exists a feasible (w/s) 


multicommodity constant flow schedule F,(t) = F, t 


1 Der 
such that P. (t): Q(to) > Q(t,). 
BROOÉ : 
Define 
E 
ak A Ji k T 
Es = SET u E (le al deal € Los *k (Ar) 
0 


Tf we start now in the state Q(t) and apply the new flow schedule 
D o a duration (t, "tp? we will transfer the system into some 


new state, say Q(t). 


k ex ex 
E (Ets Sd. Edo (A.4a) 


m 
wee Ucing A.) into (A.4a) gives: 
En . 
ak k k À 
O IA ESE) E. (tE) dt, (ab) 
= Zr er j (ži) J? 


and from (A.l) we conclude that 


A 


qe (t,) = TET T (A.5) 
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Equation (A.5) is the desired result but we must still show 


Baer as defined in (A.3), satisfies the constraints in (A.1): 


(i) obvious, integration of a non-negative (by assumption) 
function results in a non-negative value. 


(ii) From (A.4a) we have that for any t e [tot] 


^k k ^k i 

qit) = lt) m (t-tọ}rir *(i,k) < Ny (A.6) 
where 

rx g A le Tepe N (A.7) 


e > a 

A A 

Since re is not a function of time we see (using (A.5) and (A.6)) 
ER k 

that q (E) falls on the line segment Joining a and qj (t1). 

But q; (tp) and a) are non-negative (by assumption) which 


leads to 


q(t) a een Nor WE € [tort]. (A.8) 
(iii) Define 
ae (2.9) 
= k(#i) 2 
We need to show that Es, ES Enel “li € Lo: Using (A.3) we 


can write 





Hh > 
iI 
ES 
co 
BJ 
Hh 
es 
ct 
a 
or 
rt 


la] eL (A210) 
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By assumption 


k 
max Gere re.., »iaslen (A.11) 
ee ae ne Be) 0 
Ge =], 
which results in 
fij < ae ae en Lo (A212) 


and completes the proof. 


Deucerıve two corollaries of Thm. II.l. 


eonoularys 171.1 
Pa et. |Ll.i let Q(t) = 0. Then the constant flow schedule 


F(t) = F, to SENS ta is feasible in the narrow sense (n/s). 


Proof: 


Using equation (A.6) together with the requirement that 


Q(t) = 0, we have 
a E IA ke (A.13) 
re aC OA aries iene 0 ai 
Beck) EN. and thus cieseQ, (i,k) eM, and 
mp’ 2°: eo en K) O 


there is no buildup of queues. 


Somo La TL 


Pasac. let F(t), to su ti be a feasible (n/s) multi- 
commodity flow schedule. Then F (t) = F, to Bean ty is also 
feasible in a narrow sense. 
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Proof: 


From equation(A.3) and the definition of a net delivery rate, 


we have that 


K 
ct 
As 
ct 
= 
H 
En 
mM 
Z 


(A LON 


(A 


Obviously if = 63 > Ol, we 


K 
g + . 
ap 


ity, ty] and w(i,k) e Nor so is 


C 


we prove now a result which is frequently used in our study. 


Lemma II.l.l 


Let Q (Ea) and Q(t) be any two states such that Q(t) is 


reachable from Q(t). Let cht: eee 


all flow schedules such that F(t): Q(t) > Q(t). Also, let 


be the set of 


{P (8) >, Es < t < t, be the set of all constant flow schedules 


Jl 
such that Py (t): Q(ty) > Q(t,). Any Vink 11,3) =< Ly that is 


Saturatea in all flow schedules in the set {F, (t)} is saturated 
for the whole period [Eyr ts] in all flow schedules in the set 


MRE: 


Proof: 


Smpseserchatsa link [i,j] < La is saturated in all the flow 


0 
schedules in the set {F,(t)}. Let there be some schedule F(t), 


to eo ty for which this link is not saturated for the whole 


ee 


period [tort Men Oy Thm. II- l and in particular (A.10) it 


1l- 


is possible to construct a constant flow schedule Py (t), 


is = E < € 


0 suehsthat E 


res Q(t) > Q(t), puer Or which 


l l 


Zale 





the aggregate flow f.. < 


ij Cis! so that the link is not saturated. 


Dee ontradiets the definition of the set {F. (t). 


il 


DEC TOBAL OPTIMALITY AND MULTICOMMODITY FLOW SCHEDULES 

In this section we show by counterexample that the conjec- 
ture that every optimal delivery function is also globally 
optimal is false for the multicommodity case. The counterexam- 
ple applies only to flow schedules which are feasible in the 
narrow sense, i.e. do not allow intermediate data storage. 
It is still an open question whether the same conjecture holds 
for the more general class of flow schedules, namely those that 


allow for queues buildup. 


Consider the delivery problem in Fig. B.l. 


6 3 Er 
q, (0) =30 © & gg (0) =60 
(2) Cee = ly 
(5) Ca 


[i,j] € Lo 


An lo) =40 
g3(0) =70(-3 ) (6) a} (0) 


Fig. B.l. Delivery Problem 


oe definition and related issues see Section A and 
Section © Of the appendix. 
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Using the algorithm of Chapter III an optimal flow schedule 


Souto 1S obtained. It is described in Fig. B.2. 


O 28833 


OQ O 


(4,3) 


SOTS 
eo u WI 


©) 


(273) 


a E 





O 
(a) (») 





© 


‚6 (6,1) 
(4,3) (1,6) Sen 


ee 0 
Fig. B.2. Chain Flow Decomposition of F(t). 
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The optimal delivery function which is generated by this 


flow schedule is shown by the solid line in Fig. B.3. 


D(t) 
200 
a + 
eT ° - os 
Zr 
a 
138.33 ‚7 
A io 
) = 
85 2 E D, (t) 
A A 
4 ean 
Ws inn D, (t) 
2 0333 
4 
4 
h 
4 
28.33 35 44.13 55 110 t 


Optimal Delivery Function 
The delivery function D, (t) shown by the broken line in Fig. 


B.3 corresponds to the flow schedule P, (t), SNS. E which 


is depicted in Fig. B.4. 


zu 





0 <E no 





(3,4) 


33 <t 19 





(4,3) (1,6) (6,1) 


Fig. B.4. Chain Flow Decomposition of F, (t) 


We have that 


DE > e ed) (B.1) 


from which we conclude that not every optimal multicommodity 


flow schedule is also globally optimal. 
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Ge INTERMEDIATE QUEUEING OF DATA 

In this thesis most of the results concerning the deter- 
ministic delivery problem in communications networks are derived 
in terms of flow schedules that do not allow intermediate queue- 
ing of data (i.e. that are feasible in the narrow sense ). TAIS 
approach avoids unnecessary complexity which could obscure 
insight into the problem. 

For completeness, we now show that the solution algorithm 
in Chapter III can be easily modified to allow intermediate 
data storage. The implication of w/s feasibility (see (A.l)) 


is that the net delivery rate ere, 


E E ES O. (GA 
oa (Ai) 42 
is not limited any more to be a non-negative quantity. Asa 
consequence the conservation constraint (II.2) must be changed 


REO 


ae) £ o - ft J fístar- ] £í,(aJida > 0, (C.2) 
rt 


ora er to account for the non=negativity of all queues at 
all times. 

Now, it is not difficult to see that the formulation of 
MTP(m) (see (III.4)) in terms of a w/s feasible flow has the 


following form. 


+ . . . 
'This section should be read after Section A of this appendix. 


2.159 





MTP (m) (w/s): 


2 3555 


min € 
m 


(p)) Ss aa v(i,k) 


O A mM Lee 


0 
ee 


5 + ) ( u, 
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< £ No 
(2) < 310), “(i,k) e N 
E Si r 1 = 0’ 
K re Ly 
Soo, “[i,j] « Lo (C. 3a} 
art Gij” 
Zee e Los 
Dr il, a” 
lu = 0 
DO 
we = At r 
C R p Po 
n= Aa ‚m-2 
0 
= De 
a D aem ee Nos 
ee Bee 


DS Le al 





0750 0 0 
for given tyr Py ree ert 1 m1: 


Lat 


inspection of (C.3) reveals that it differs from the original 
statement of MTP(m) in the form of one constraint only (under- 
lined). A similar modification can be applied to MRP (m) 


(III.6) to allow for intermediate data storace. 


MRP (m) (w/s): min ) en 


na- l 


0 k k Ok er | f 
(t,7e)r, (mel) + er, (m) + I, Sa) = q, (0), ~(1,k) < No 
0 k k | 
(te) xr, (mtl) < q;(0), ¥(i,k) < Ny 
0 k | k K : 
= ; ; e N y 
(€, 7elry (mtl) + er, (m) < g,(0), (i,k) o © ) 
m-l 
0 K k O k k : 
— o à A , , = N i 
er E er, (m) + Èn E O < G70), ¥ (i,k) o 
n=m-1,m2,.: 2 
y £5. (p) ae) = Lee 
P =i a ‚m+l 
) ri (p) = fon A 2na 
(i,k)eN P 
0 
£$. (p), rip) A Kae Na 
2, / zer = 0 


¥[i,j] e L 


ie ee es oy ee 
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oa. 0 O 
for given O EE E 
where 

K k k 

ete) = ) £,(P) - ) £,,(p), v(isk) € Ny, p = 1,2,.-.,mtL 
j (ži) C 
O _ ODO a 3 
En = = De DIPS LL mel 

and 

e is any real number such that 0 < « < a = n 

m m+1 


The change in the form of the cost function is necessary to en- 
sure that we minimize delivery rate and not tne total flow rate 
in the network. 

We conclude this section with an example of a simple delivery 
problem for which we compare the delivery functions resulting 


from both types of flow schedules. 


Example: 
3 4 
= 0) =6 
dy (0) =18) AS ew 
oo © 
a] (0) =6 | 


capacities are indicated 
on links 


Doel Delivery Problem 
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An optimal flow schedule solution for the problem in Fig. 


Bee echown in Figs. C.2a and C.2D. 
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AE Zac nchaingElow Decomposition for the Perlod Kor 
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ay (1) E 





© > 
© Me 


EE ER Chain Flow Decomposition for the Period (1,4 ] 

The resulting optimal delivery function is shown in Fig. 

Suppose now that we allow for intermediate data cueueing 
and in particular consider the flow schedule in Figs. C.da and 


Observe that in the first period (0,11 commodity (1,7) 


flows with rate six from node 1 to node 6 and continues 
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D(t) 


42 
>, =10 
>> =12 
12 
0 1 2 3 
Bremse 3. Optimal Delivery Function 


ay (0) 


(intermediate 


ag (t) 
storage) 





Dee. a 
Queueing for the Period [0,1] 


225 
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Chain Flow Decomposition with Intermediate 








2 


t~ 


Fig. C.4b. Chain Flow Decomposition with Intermediate 
Queueing for the Period (1,4] 


with rate three from node 6 to node 7. As a result, three units 
of commodity (1,7) are stored in node 6 in an intermediate 

queue denoted by delt). The contents of this queue are delivered 
in the second period (1,4]. 

The delivery function which is generated by the new flow 
schedule is superior (dominates) the optimal delivery function. 
Both delivery functions are displayed in Fig. C.5. 

We conclude that flow schedules which allow intermediate 
data queueing may have advantages in certain instances over 


flow schedules which are feasible in the narrow sense. Never- 


theless, it seems clear that wide-sense (w/s) feasible optimal 
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ui 
= p, =10 
1 
15 ae en 
u f CE — (w/s) a 
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EE. Comparison ot Delivery Functions 


flow schedules (which allow for intermediate data queueing) 
Ano: produce delivery functions that differ substantially in 
their basic characteristic (piecewise linear, convex, etc.) 


from those produced by narrow sense (n/s) optimal flow schedules. 


Meee Mone ON STABILITY 
In this section we complete the discussion of stability 
(from Chapter IV.A.c) with examples of delivery problems that 


melustrate this concept. 


l. Unstable DELI Veny Froblem 


3 
a, (0) q; (0) 


ij 20 


Fig. D.1. Delivery Problem 
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3 
For the case when qı (0) = q3 (0) = q, an optimal con- 
ge CEER low schedule solution can be found by inspection. It 


Momsiown in Fig. D.2, 


Fig. D.2. Optimal Constant Flow Schedule 


Formally, the MTF(1l) for this problem can be written as: 


MTP (1)": 


0 Peo 1 0 0 oj |t | 
0 LL 0 0 0 0 Ur) 
-1 | 0 0 1 0 0 4,2 = 0 (D71) 
=l 0 10 0 1 0 | ne o! 
E o o uf . 0 
E25 
513] 


| Y 
[O 


(tiru 27u23 Y1 37 5127823513) 


el- 


'Since both commodities have a common destination (node 3), 
wWemwill mot use the upper index notation to indicate the destina-— 


Prom, 2.e. ze > De: ee Lo: 
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An optimal basic solution ze iS round to be 


T 
er (£1 75370) 378) 9/8) 3) ee ee (D2) 


and the related optimal dual solution is 


where 


I, = O O OOO 
Now, suppose that a, (0) ge ae Wa ES MOE dia 
cult to see that in the new optimal flow schedule, commodity 
meee 8 muse use the link chain ({1,2],{2,3]) in addition to 


Is] Therefore it has the form shown in Fig. D.3. Also, 


DeD Sciructure Of the Optimal Perturbed Flow Schedule 


aca di that the flow rate variables Ej are related to the 
A PS by the transformation 


ze 0 pr 
fij = uij/ tir ee Lo. 
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memset be true that 


Oe q | € 
a E = mo E Oey 








the change AZ in the Optimal value of the cost function 


that was caused by the perturbation is 


(D.6) 


NN 


AE th = da» (a = = 


If we use the optimal dual solution (D.3) to evaluate Az, we 


obtain 
fae oe = (0,1,0,0,-1)-(0, 0,00). = €, (D.7) 


which is incorrect. We conclude that the original optimal 
solution is not stable (optimal dual is not unique). Actually, 
ehesriew optimal basic solution is (cf. Figs. D.3 and D.5) 


£ s E E E 
X = (Ey 1077173107, 31512) = (q 3,3,973,975,9 E) (D.8) 


"Using the law of proportions, 


a c _ ate 


b d b+d° 
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and the new related optimal dual is 








Me A A (D.9) 
where 
2 = 
Bp Pa: 
Obviously, 
zen (D.10) 
2. Stable Delivery Problem 
4 4 
q (0) q, (0) 
ir ZA] < L 
Fig. D.4. Delivery Problem 
4 4 ) 
For the case when q7 (0) = 2q and q, (0) =g, an optimal 


constant flow schedule solution can be found by inspection. 


Ms shown 1n Fig. D.5. 


Zou 








Fig. D.5. Optimal Constant Flow Schedule 


Bormallv, ene MTP (1) for this problem can be written as: 


MEP{l): min ty 
So. Es 
O OO 00 00 0 0 Ol |e, | 2q 
O OO 0 oe u, q 
ro ro oo 0 0 0 0 lu. 0 
OO o 1 0 0 0 0 0) m, 0 | 
A O 0 0 0 1 0 0 0 Ola, = 0 | (D.11) 
Owes 6 6 GO 1 6 0 ol fu 0 
Peon o 8 O O LT 0 Ol faa | 0 
-ı 0 0.00 (0 0 0 0 0 Ot o ss, | o | 
0 o o oo Li sja o] 
S93, 
A 
za 
Is, 
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it a 34 5191513159315 1415241534) 2 


Pmsortimal basic solution X, is found to be 


B 


Xp ME Uy Uy 31 Ug 37 By gr Ugg 4 78) 578)3) 


nana, q, 2,9,0,0) , 


and tne related optimal dual solution is 


where 


Now, suppose that q7 (0) = 2qgte, £ > 


thought the reader will convince himself that 
else link chain ([1,3],13,4]) in addition 


([1,2],{[2,4]). Therefore an optimal solution 


Sewn in Fig. D.6. 


Also, it must be true that 











AE q E 
a 7 a 
We can now calculate a from (D.14) 
e 
DB ale g = 2 
ee 177 
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|o 


(D T2) 


(DS) 


with a little 
commodity (1,4) 
to ([1;4]) and 


has the form 


(D.14) 


wim 


(DLS) 








Fig. D.6. Structure of the Optimal Perturbed Flow Schedule 


tnewenange Az in the optimal value of the cost function 


that was caused by the perturbation is 
= bs = en 
t. = gq (q +3) 3 (Dike 


If we use the optimal dual solution (D.14) to evaluate Az, 


we obtain 


1 
bof 


— e 1 1 sl —z 
ze = 7 N° AD ‚3,3,0,0,0, 3, 
z ORT 


which is the correct value. Actually, the new basic solution 


Demeter Fig. D6 and -D.15 ) 


PA 


Bp = (€¡+077%] 3.13% 411941 U3415]3:823) 


= 


age E € € G 
(G43,943,3,9,94+3, 95,437, 7) (D12) 
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and the new related optimal dual is 


= 
il 
> 
uu > 
lI 
wm 
wol m 
~~ 
wj — 


‚ur.0,.0,- 


Lo] 
wm 
wm 


a he (DTO) 


where 


Thus 
Du =. A. (D. 20) 


Formally, we have not shown that the original delivery 
problem in this paragraph is stable (i.e. has a unique optimal 
dual) but rather have illustrated that although the perturbation 
causes a change in the optimal basis, the optimal dual solution 
is not changed, a behaviour which is characteristic of a stable 
Ferne. Le 1S worth noting that further study of this example 
would show that the only unstable point here results from a 
requirement vector b, such that q% (0) - 59 and q5 (0) = q. The 
corresponding flow schedule is shown in Fig. D.7. In general, 
for a given network we expect a randomly selected requirement 
vector b (queues sizes) to be stable though the optimal primal 
solution may be degenerate. This observation is backed up by 
the experience we have gained in solving a number of delivery 
problems. This serves as an additional motivation why we are 


interested more in stability than in non-degeneracy. 
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Eee, Unstable Optimal Flow Schedule 


ORTI MAL DELIVERY FUNCTION IS PIECEWISE LINEAR 

In this section we show that the assumption we made about 
the piecewise linearity of an optimal delivery function is jus- 
tified. Assume to the contrary, that there exists a continuous, 
non-linear optimal delivery function. This function must be 
convex, since otherwise we could improve on it by generating 
its convex hull by the method of constant flow substitution 
(see Appendix I.A). A typical delivery function of this kind 
is shown in Pig. E.l. The broken line in Fig. E.l represents 
the delivery function which corresponds to an optimal solution 
COSTEL). 

consider the optimal flow schedule F(t), 0 < t < = that 


K 
Ben zstves Dit), and in particular the net delivery rates r,(t), 
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D(t) 





Fig. E.l. Optimal Continuous Non-linear Delivery Function 
*(1,K) e No at time t = oe Define a perturbation vector AQ 
such that 


E (E.1) 


1% 
where 
O< 6 < $11 


and 6, is the maximal value of 6 for which the perturbation AQ 


I 
is acceptable (for simplicity we assume that the problem is 
at a stable point and thus $1 > 0: 


The optimal value of the cost function for the perturbed 


problem is tle where (IV.20) 
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e = o (Ag (E.2) 
(1,k) eN, = 


From Thm. 1V.4 we have that any flow pattern of commodities 


in N, that saturates the set Lı (and in particular the flows 


R a Be must satisfy the perturbation equation, i.e. 
1 = eee (E.3) 
(i,k) eN} 


Weamoa (E.1) and (E.2) with (E.3), we conclude that 
Be 6 (E.4) 


Equality (E.4) ensures us that it is possible to deliver 


the set of queues Ags - Be ¥(1,k) € N) within the period 


(&)-e,t)] by using the set of feasible flows (£p5 (tq)}, 


Pei) «© N of the optimal flow schedule Fit), 09 < = < SL The 


a 


perturbed (two segment) delivery function is shown in Fig. E.2. 


D(t) 





Fig. E.2. Perturbed Delivery Function 


DS 





Since D, (t) dominates D(t), the last cannot be an optimal 
delivery function. 

Following the same argument it can be shown that an optimal 
delivery function must be composed only of linear segments, 
i.e. it is not possible to have an optimal delivery function 


as shown in Fig. E.3. 


D(t) 






(2) = linear 
segment 
(1) = non-linear 


segment 


Fig. E.3. Mixed Type Optimal Delivery Function 


Basically, this completes our argument about piecewise 
linearity of optimal delivery functions. We note that in case 
the optimal solution to MIP (1) is not stable, the argument does 
echange significantly: 

imsteadmor sh. 2) weewrite 


k k 
e = ) o,(1,R) Ag,» (es) 


(i,k) eN} 
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where 


K0 


= f ] = 
R Ej (ty), (ES) e N 


1* 

mare dual variables SE A € Ny in (B.5) \correscend 

to the related optimal dual solution for the perturbed problem 
k k Tae 

(note, 9; (1,8) = a (1), *(1,k) € Ny ir, the Original solution 

is stable). All other properties remain as before. 


We summarize our discussion with the following theorem. 


Theorem E.l 


An optimal delivery function is piecewise linear. 


E Ohman NUMBER OF CORNER POINTS 

In Chapter III.B.2 we proposed a conjecture (Conjecture 
III.1) which upper-bounds the number of corner points of an 
optimal delivery function with Nol” the number of non-zero 
data queues. The conjecture is based on one aspect of our 
experience with delivery problems, the essence of which may be 


formalized as follows. 


Hao position F.l 
At a stable point, let K n denote the number of commodities 


for which the dual variables at the m-th corner are equal to 


their maximal value or are negative, i.e., 


| ko k | 
K, = |((1,k): 0,(m) =0,,,¿(m) or 0, (m) < e 
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then 


Conjecture III.l is a trivial consequence of Proposition 
F.l. Actually, we proved a related result (cf. Lemma IV.5) 
which says that if K = [Noi for some m, then the delivery func- 
tion cannot be improved any more and thus is optimal. There 
seems to be a slight difficulty in proving the opposite direc- 
Bean, ise, that an optimal delivery function implies (at a 
stable point) equal dual variables. We must leave the proof 
of this conjecture and of Proposition F.l as open topics for 


further research. 


COMPUTER SOLUTION EXAMPLE OF OPTIMAL DELIVERY PROBLEM 





4 1 3 Dein ae 
0) =70, 44 (0) =35, 45 (0) =50, 43(0) =25, 44 (0) =25, 45 (0) =35. 


Fig. G.1. Optimal Delivery Problem 
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The problem in Fig. B.l1 was solved using the algorithm of 
Chapter III. We will not present the partial solutions in de- 
tail, but only the optimal flow schedule and its corresponding 
delivery function. The optimal flow schedule is composed of 
four segments, the chain flow decomposition of each is shown 
in the following figures. We append each segment with the 
elevant information obtained from the solution of the corres- 
ponding minimal time problem (recall that in our notation the 
0 


"first" segment is G. Also, we demonstrate the pertur- 


bation equation for each one of the flow segments. 


MTEP (1): N = MS dr LS, L = o ley 






T l 
5 1.4 eee a 
2Sa see TD = (oc) (1) =7,0, (1) =7,903(1) =7, 
5 ei 
o,(1)=5,0,...,0) 
k k 1 
) r = 5-2 =1 
A = 
e 


Fig. G.la. Chain Flow Decomposition for t « (507515251 


242 





MTP (2): N, = N}; L, Ly 
B PRT ee. 
2) = (0) (2)=1,05(2)=1,93(2)=2,03(2)=2,0,...,0 


re) 





q(t) 


ga (t) 






ko. k 
— o 
OS 


z =(0.8-1 A) Doula oooh eee al 


0_ 3 


ice Globe chain Flow Decomposition for t « (31.25,50] 
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MTP(3): N ES 2 Ms U loser)! 





I 
JE 
3 
HE 
I 
q, (t) 
1 Be ini. 2... 2: 
Iss (2) la een o. (3)r (3) = 355 2+5 1 +>5-1 +5-1) L A 5 


0 


Sc aia Flow Decomposition for t e (25,31.25] 
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MTP (4): N. = N, A L. =L ao 


4 4 E 

_ 3 E SS 
) (4) = (97 (4)=1,0),(4)=1,03(4)=1 07 (4)=1,0(4)=1, 
07(4)=1,0,...,0,0,(4)=-5) 





L 
t 
= ¿1.7 +] +1 +0.15 +1.15 +1) = 1 
0 © q% (t) 
O 4 


cae na in Flow Decomposition for t e [0,25] 
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The corresponding optimal delivery function is shown in 


Fig. G.le. 
D(t) 
Je 
237 oo 
0) = 
180 05 =5 
150 o9 aa 
0 25 3l L5 DUO Zo 


Fig. G.le. Optimal Delivery Function 


Poe DE TALGE@ SOLUTION OF THE MDDP EXAMPLE 


For completeness we restate the MDDP. 
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tJ 0 
“hs = 1, #l1,j] € Ly 

cy E a. =b, =œ, wi eV 
a T 


Q = 2020707070,0)7, a 
Q) =207,.070,0,07 20), 
Q- =a OFO 155,0). 


Fig. H.l. Maximally Delayed Decision Problem 


«ba 


Optimal Flow Schedule: | 








aE (kK) 

TR period flow 

a] Il 2 
2 2 2 
úl 8 2 

Me E 0 
et 2 0 
4 3 0 


"We assume here that the Seamemg cima sls k = l, and- not 
k = 0 as before. 
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[5,4] 


[5,6] 


period 
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flow 








period 
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jene period flow 





[5,4] 4 0 
i 5 0 
: 6 0 
i 7 0 
g 8 0 
" 9 0 
(S55) 4 2 
A 5 0 
i 6 2 
dl 7 2 
e 8 2 
i 9 2 
F, C 
[1,2] 4 2 
Á 5 2 
i! 6 2 
‘t 7 2 
i 8 2 
Be 4 0 
i 5 0 
3 6 0 
‘ 7 2 
i 8 2 
al 4 0 
i 5 0 
" 6 2 





IR period flow 





i273) 7 0 
: 8 0 
ee 4 2 
4 5 2 
£ 6 2 
N I I 
" 8 0 
9 2 
(245) 4 0 
: 5 0 
i 6 0 
y 7 0 
i 8 ml 
á 9 2 
Sol 4 0 
n 5 0 
i 6 0 
i 7 2 
t 8 2 
1 9 2 
(Son 4 0 
f 5 0 
4 6 0 
: 7 0 
y 8 0 
i 9 0 
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